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ES 
Avant-Propos 


Cet ouvrage contient l’essentiel du matériel couvert dans le cours de Mécanique Clas- 
sique II (PHY-10492). Il est basé sur les notes de cours de P Amiot et prennent leur 
inspiration comme il est coutume de plusieurs livres de références. 

Les notes couvrent la mécanique classique avancée, soit le formalisme de Lagrange, le 
formalisme canonique, la théorie des perturbation et le mouvement d’un corps rigide. 
Les notions de mécanique sont rappelées dans le chapitre 1. Le formalisme de Lagrange 
est introduit au Chapitre 2. Suivent quelques applications et propriétés (Chapitre 3), le 
formalisme canonique (Chapitre 4), la théorie des perturbations (Chapitre 5) et finale- 
ment le mouvement d’un corps rigide (Chapitre 6). LC appendice contient un résumé des 
notations, un aide-mémoire et quelques références complémentaires. 


Québec Luc Marleau 
Mai 1997 Département de Physique 
Université Laval 
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1 RAPPEL 


1.1 Trajectoire et cinématique d’une particule ponctuelle 
| 


La particule ponctuelle est sans dimension. C’est une création de l’esprit, un modèle, 
représentant un objet physique qui n’est animé que d’un mouvement de translation (pas 
de rotation sur lui-même). On admet ici que notre espace physique est à trois dimensions 
auquel on adjoint le temps qui n’est pas ici une dimension mais un paramètre immuable 
et indépendant des objets physique et de leur évaluation dont il sert à mesurer le taux. 

Nous représentons l’espace physique par un espace à trois dimensions à l’échelle, 
doté d’une origine notée O et de trois axes orientés. La position instantanée de la parti- 
cule y est notée par un point P dont la position est entièrement définie par un triplet de 
nombres appelés coordonnées du point et qui mesurent généralement des longueurs ou 
des angles (voir figure 1.1). Ces coordonnées seront souvent notées x; ou q;. Il est sou- 
vent pratique de parler du vecteur position de la particule, noté x ou p qui va de l’origine 
O au point P. 


Figure 1.1 Trajet d’une particule 


L'évaluation du système physique sera décrite par une courbe ou trajectoire C, décri- 
vant le déplacement continu du point P dans notre espace de configuration. On conçoit 
cette évolution comme résultant d’un paramètre invariant qui augmente. On le choisit 
généralement et pour des raisons pratiques comme étant le temps, noté t, mais ce choix 
n’est pas unique. Le point P se déplaçant avec le temps sa position, r, variera dans le 
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temps et la trajectoire sera décrite par r = r(t) en terme des composantes par: 


m=m(t), i=1,2,3. (1.1) 


Qui dit mouvement pense intuitivement à une rapidité de mouvement. Cette notion, 
ce concept est quantifié par la définition de la vitesse V 


V(t) = Lx) = x(#). (1.2) 


Notons par la lettre p le paramètre (arbitraire) dont la variation génère la trajectoire (il 
peut être ou non le temps). Alors la longueur s de la trajectoire entre po et p1, est donnée 
par: 


(1.3) 
où p varie de façon monotone entre po et p1. Alors on peut écrire (voir figure 1.2): 
dx ds dx dx 
= = = V— 1.4 
dt dtds ‘ds oha 
— 
\ 
Figure 1.2 
On voit immédiatement que : 
d Pa 
Z? (1.5) 


ds 
un vecteur unitaire dans la direction du vecteur T qui donne la tangente à la trajectoire 
au point P. En effet 
T= lim — = —. 1.6 
, Ji P As>0 As ds { ) 
On obtient ainsi V =Tv ou 7 donne la direction et v la grandeur de la vitesse (vectorielle) 


V. Par abus de langage v s’appelle aussi la vitesse. Ce qu’il faut souligner, c’est que V 
est toujours tangent (c’est un vecteur) à la trajectoire. D’ailleurs, pourvu que le paramètre 
p varie de façon monotone (et continue) le vecteur = est tangent à la trajectoire, le cas 
V= de n’est qu’un cas particulier. 

Intuitivement la vitesse V peut varier le long de la trajectoire (voir figure 1.3). Pour 
quantifier cet effet nous définissons l’accélération a 
dV dx 


a 
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et clairement 


Pop dV  d(Fv) 
= dt dt 
dv. d? 
= — — 1.8 
d tou (1.8) 
Parce que 7 - 7 = 1 alors ae? T D =9F -2 = = (. Ainsi d? est perpendiculaire à 7 qui 


est tangent à la trajectoire. Donc %& T est normal à cette re Appelons ñ le vecteur 
unitaire normal à la trajectoire dns la direction de £ i.e. dans le plan instantané de la 
trajectoire). On calcule 


dr dî ds dF? a 


2 = 1.9 
a la Ga -F Th. a) 
On écrit par définition, p7! = pa . On a donc pour a 
E de 
= — — 1.10 
a : ñ + de F: ( ) 


Ainsi l accélération a une composante tangente à la trajectoire (F) de valeur 4 ce et 


Figure 1.3 


une composante normale à la trajectoire (f) de valeur + . On peut montrer que p est le 
rayon de courbure de la trajectoire. En effet, dans le voisinage immédiat du point P, la 
trajectoire peut être approximée par un arc de cercle, p serait alors le rayon de ce cercle. 
Plus la trajectoire est courbée autour de P, plus la vitesse ne rapidement selon ñ. 
De fait, plus p sera petit et plus la composante normale de a, a , sera grande. 


1.2 Plusieurs particules ponctuelles 


Pour représenter la position de N particules dans notre espace de configuration à 3 
dimensions nous avons besoin de N triplets de nombres (total 3N) 


r, = (tir, 2v2, Pa): W= 1,2, …, N. (1.11) 


L'évaluation d’un tel système sera représentée par N trajectoires (une par particule) dans 
cet espace. 
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Il est souvent utile d’imaginer un espace abstrait comptant 3N dimensions, 3N co- 
ordonnées y sont nécessaires pour décrire la position d’un point de cet espace qui donne 
à lui seul la position instantanée des N particules. Par un léger abus de notation on note 
les coordonnées de ce point {x;; i = 1, 2, ..., n = 3N }et on peut parler de la trajectoire 
du système dans cet espace. 

Ainsi, assez typiquement on écrira alors des expressions comme la force par exemple : 


F(x;,t) (°% composante); i,j —1,2,...,n. (1.12) 


1.3 Éléments de dynamique 


Depuis Newton on connaît l’équation fondamentale du mouvement : 
F = ma. (1.13) 
Elle prend plusieurs formes (pas nécessairement équivalentes) 


dr dv dp 
dt? 7. dt 7 dt 

La quantité F est la force. Elle détermine le système et est déterminée empiriquement, 

i.e. c’est l’expérience qui nous en donne l’expression. 


Cette expression qui est vraie pour 


(1.14) 


r = X = (£1, £2, £3) (1.15) 
le demeure pour un nombre n de degrés de liberté. Pour alléger, écrivons 
x= (£1, 22, L3, T4,- En) (1.16) 


un vecteur à n composantes. Intégrant m dans F (qui n’aura plus les dimensions d’une 
force mais celles d’une accélération) écrivons l’ opération de Newton : 


X = f(x,%å,t); n composantes, n = 3N (1.17) 
ži = filxj, t,t); n équations, i = 1, 2, ...N (1.18) 

ou encore 
X= fo (Xp žpt); v,u=1, 2, ...N particules. (1.19) 


L équation de Newton, en tant que loi physique se doit d’obéir à certaines symétries 
que nous fait découvrir l’observation de la nature. On dit alors que la mécanique classique 
doit être invariante sous les transformations de Galilée. Cette invariance est valable pour 
les systèmes physiques fermés. Il n’y a qu’un seul tel système, c’est l’ Univers mais en 
pratique les effets des corps éloignés sont souvent négligeables et on fait l’approximation 
que le système est fermé. Cela signifie que tous les corps qui jouent un rôle significatif 
sur le système sont inclus dans le système. Il n’y a pas de force extérieure. Cette dernière 
notion de force extérieure peut également être utile, mais nous y reviendrons. 

L'étude d’un système physique peut se faite entre to et t ou entre to + s et t + s (on 
peut refaire aujourd’hui une expérience faite hier et obtenir les mêmes résultats). Ainsi, 


T= filzj, t,t) = file, t +s) (1.20) 
où s est quelconque. On on conclut que f ne peut dépendre du temps et donc 
ți = filxj, t); n = 3N équations. (1.21) 
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Postulons que les résultats d’une expérience sont indépendants de l’endroit où elle 
est faite. Si je déplace d’une même distance orientée, l, chaque particule du système 
physique alors sa position passe de x, à x, + 1 (v compte les particules) alors que x, 
demeure x, puisque 1 = 0. La loi de Newton x, = f,(x à Š) doit être indépendante de 
1, ce qui impose que f, dépende de x, sous la forme x, — xx puisque 


Xa — X) > (x, +1) — (x + l) =x, — x) (1.22) 


donc 
X, = fp (Xp — X\,Xy). (1.23) 


On sait également par expérience que la physique est la même pour deux observateurs 
se déplaçant l’un par rapport à l’autre avec une vitesse constante (translation de vitesses). 
Cela impose soit 

fo (xa — Xa) 
f, = ou (1.24) 
fy (x, a XX, Xu due: xı). 

On admet également que la physique au Canada est la même qu’en Australie, même 
s’ils ont la tête en bas. Par conséquent les lois physiques, telles l’équation de Newton 


ne peut pas dépendre de l’orientation de notre système de référence. Un tel changement 
d’un angle @ autour d’un axe f se note, en coordonnées cartésiennes 


r> xr (1.25) 


ou, si on écrit r sous forme (matricielle) d’un vecteur où les éléments sont les compo- 
santes de r, 


r —Gr; où G = matrice 3 x 3 pour une particule (1.26) 
Clairement, si 
r —Gr (1.27) 
alors 
r—Gi et r—G* (1.28) 
donc l’invariance de 
F=f(ri) = r=f(Gr,Gi) (1.29) 
implique 
f(Gr,Gi) =Gf (r, ù). (1.30) 


Complétons tout cela avec les autres lois de Newton avant de revenir plus tard sur 
certaines conséquences des résultats ci-dessus. Dans un système fermé, la loi d’action- 
réaction stipule que si un corps, noté par l'indice v agit avec une force F, sur un corps 


u alors ce corps agit sur avec une force Fy = —F,,. Ainsi si nous n’avons que deux 
corps, avec ry = (£y1, L2, y3) 
Mi = Fu (1.31) 
mo) = EF: = —F;2 (1.32) 
ou de façon générale, pour N corps (sans somme sur v) 
N N 
mË, = 5 Fp =- 5 Fv. (1.33) 
u=1 u=1 


Cette loi a une conséquence immédiate et importante: la conservation du moment 
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total. Sommons ci-dessus sur v 


N 
> mn Nr =0 (1.34) 
2 LU 
donc 


2. = 7 D mir = = 0. (1.35) 


La quantité dérivée est donc une constante dans le temps, i.e. >. Myly =C. 
Il est habituel de définir le moment p, = mt,. Nous aurons donc 


SP —=C=P: le moment total. (1.36) 
Remarque 1 
En conclusion : le moment (linéaire) total d’un système fermé est une constante du mou- 
D vement. 
On définit le moment angulaire d’une particule par 
l =f; Xp, = mr, xX i, (1.37) 
donc 
I = moi xi, + mr, xË, =0+m r, xË, (1.38) 
N 
= r xF, =r,x D Es (1.39) 
H 


Définissant le moment angulaire total du système 
L=, (1.40) 


alors 


N N 
L=¢ `r, xò Fu). x Fpp: (1.41) 


H ; 
Avec F,, = 0 (la particule n’ agit pas sur elle-même). 

Or, le vecteur r, — r, est dans la direction relirant les particules v et u. Si la force 
entre ces particules est dans cette direction, comme sur la figure 1.4, alors le produit (x ) 
sera zéro et L = 0 donc L = constante. 


Ev 
Put 
v 


Figure 1.4 
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Remarque 2 
Par conséquent : si les particules constituant un système fermé n'agissent les unes sur 
les autres que selon la droite qui les relie, alors le moment angulaire total du système est 


une constante du mouvement. 


1.4 Travail et Énergie 


Lorsqu'une force F agit sur un système physique, disons une particule, on dit qu’elle 
fait un travail sur ce système. Ceci cause un changement de l’énergie de ce système. Soit 
une trajectoire entre les temps to et t. Calculons le long de cette trajectoire la quantité 


F-dx 
z(t) t , tp 
| Fdx = [ra m [ ee. 
s H a de dt 
m f* d {dx dx m dy 
= PJ See) 24 
ft IL | 50 
0 0 
l 0 l 2 
= 3y (t)— LM (to) =T — To. (1.42) 


Appelant T = mv? Ténergie cinétique, on voit que l’application de la force F se 
traduit par un changement de cette énergie cinétique. Notons cependant que l’ intégrale ci- 
dessous se fait le long d’une trajectoire. Le résultat peut donc dépendre de cette trajectoire 
(voir figure 1.5), i.e. de façon générale 


Î Fdx + F'4x. (1.43) 
Ci C2 


Dans certains cas cependant, et ils sont physiquement importants, l’intégrale ne 


C; 


Figure 1.5 


dépend pas de la trajectoire mais uniquement des points initial et final, on dit qu’elle est 
conservatrice (la force). Strictement parlant, il s’agit d’une propriété mathématique, i.e. 
qui résulte de la façon dont F dépend de x, v, t. Il se trouve que dans monde physique 
réel, plusieurs forces peuvent être décrites par de telles fonctions. Lorsque tel est le cas, 
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l'intégrale de F-dx sur un parcours fermé est évidemment nul. 
0= f Fax | VxFdS (1.44) 
C SEC 


où l’application du théorème de Stokes est responsable de la dernière branche de cette 
équation avec S' une surface dont la courbe fermée C marque la frontière. Comme cette 
surface est arbitraire mais que le résultat de l'intégrale doit toujours être nul alors la 
fonction à intégrer doit être nulle 


V xF=0: force conservatrice. (1.45) 


Dans ce cas il est toujours possible d’écrire F comme le gradient d’une fonction scalaire. 
On écrit 

F =-VV(x) (1.46) 
et on appelle V(x) l'énergie potentielle. Ainsi le travail fait par une telle force entre les 
points Xọ et x sera 


L F-dx = — n VV(x) = V (xo) — V(x) = V — V. (1.47) 


On avait vu que ce même travail était donné par T(x) — T (xo). Nous aurons donc 


T+V =D +V : Énergie conservée. (1.48) 
Lorsque la force qui agit sur une particule est conservatrice on peut définir une constante 
du mouvement (indépendante de t) qu’on appelle l’énergie E = T + V. Physiquement 
la force est donnée par — V V, on peut donc remplacer V par V + constante sans changer 
la force F. On change alors la valeur de E en E+ constante. L échelle d’énergie ne peut 
donc être fixée qu’à une constante additive près. En pratique on fixe la valeur de V(x) 
à une certaine valeur, Vo, pour une x = Xo, Xo et Vo étant arbitraires. 


1.5 Systèmes à N particules et forces extérieures 


Supposons un ensemble de N particules interagissant entre elles et sur lesquelles 
peuvent également agir des forces extérieures. Notons m; la masse de la à" particule, 
F; la force externe qui agit sur elle et F;; la force due à l'interaction de la j°* particule 


sur la iiè"%, Évidemment F;; = 0 et par la troisième loi de Newton F;; = —F;;. Pour 
la 47% particule, l’équation de mouvement est 
miki = F; + X Fij. (1.49) 


j 
Sommant sur toutes les particules 


5 MX; = 5 F; + 5 F;; = 5 F;=F: force externe totale (1.50) 
i i ij 


(3 


parce que Dn j Fij = 0. Avec M = X; m; : masse totale des N particules, 


1 . & |1 
F =M pri mx, = Me Z Drz (1.51) 


1.5 Systèmes à N particules et forces extérieures 9 


d’où 2 
N 1 


donne la position du centre de masse du système. Le mouvement du centre de masse se 
fait comme si toute la masse y était concentrée et que la force externe totale s’y appliquait, 
quelle que soit l'interaction entre les particules. Définissant le moment linéaire total où 
P = MX, on aura 


d ` 
F= ge ÀP Sh mx. (1.53) 
Si la force extérieure disparaît, alors P = constante. 


Après le moment linéaire total, étudions le moment angulaire total. Nous aurons évi- 
demment par rapport à l’origine Ox 


N 
L= X mixi x X; (1.54) 


mesuré à partir de l’origine du système de coordonnées utilisées. Il est utile d’utiliser les 
coordonnées relatives que nous noterons les y; (aucun rapport avec le y des coordonnées 
cartésiennes), et définis par y; = x; — X — x; = X + y,, 


N N 
L X mixi X Xi = X mi (X + y;) x (X+:) 


N 
= Dmfrixi+xxX+X x ÿi+y x X) 


(1.55) 


mais $; miy; = 0 donc $; m;ÿ; = 0 aussi, et alors 
N 
L = X miy: x ÿi+ MX x X = L, + Lom (1.56) 


où L, = D miyi x yiet Len = MX x X = X xP. 
Ainsi le moment angulaire total par rapport à l’origine d’un système inertiel est la 
somme vectorielle du mouvement angulaire relatif des particules par rapport au CM et 
d’un moment angulaire correspondant à la totalité de la masse centrée au CM par rapport 
à l’origine du système inertiel. 
On peut passer d’un ensemble de particules ponctuelles à un corps de volume fini 
en remplaçant de façon adéquate les sommes par des intégrales. Dans ce cas on voit 


apparaître des densité de masse p(x) telles que 


M = p(x)déx. (1.57) 


Volume 


Exemple 1.1 


Système simple unidimensionnel: 
Si la force F = F(x) et qu’en une dimension il existe une fonction V (x) telle que 
ð OV 

V 


(1.58) 
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Supposons V (x) comme sur la figure 1.6 et étudions une particule qui serait soumise à une telle 
force. Nous avons 


E= Ta +V(a) =T +V. (1.59) 
Évidemment T > 0 et donc E > V(x) toujours. Donc E > Eo. Ceci contraint le mouvement. 
Par exemple si E = F, alors le mouvement sera limité à la région entre x1 et £2 . Par contre 


si E = F», alors non seulement la région £o < x < xs est-elle possible mais aussi la région 
£ > Ta. 


Figure 1.6 


En une dimension il est simple d’obtenir la solution à partir de l’ équation pour l’é- 
nergie ci-dessus. En effet, isolant de =f; 


dx 2 1 
= ZE- vo) (1.60) 
m dx 
dt = | — ————. 1.61 
V 2 /E-V(x) Ra 


Intégrant, 


1.62 
TEG E-— VEO ue 


ou formellement t — to = f(x, E) — f (£o, E) où on isole z = z(t — to, E, xo) : solution 
unique si on connaît FE et zo = a. 


1.6 Degrés de liberté 


La notion de degré de liberté jouera un rôle important dans les chapitres qui vont 
suivre. Cette section est consacrée à la première étape de cette notion. 


1.6 Degrés de liberté 11 


La première caractéristique des degrés de liberté est qu’ils se comptent. Un système 
physique a un, deux, trois,..., N degrés de liberté. 

Le degré de liberté est généralisation du nombre de directions indépendantes selon 
lesquelles une particule peut se déplacer dans l’espace physique. Ainsi, une particule 
ponctuelle pouvant se déplacer dans une direction possède un degré de liberté; elle en 
possède deux si elle peut se déplacer dans un espace à deux dimensions , etc... . Des 
forces agissant selon une ou plusieurs de ces directions peuvent limiter le mouvement 
de la particule à un domaine fini selon ces directions sans faire disparaître le degré de 
liberté. Par exemple, si une particule est libre de se déplacer selon l’axe Ox seulement, 
elle a un degré de liberté. Si une force, disons harmonique, Fy = —kzx, agit sur la parti- 
cule, le domaine de variation de la particule sera réduit de —x0 à +x0 selon son énergie 


E = koi, et la particule a toujours un degré de liberté. Cependant si cette force est ca- 
ractérisée par une tige rigide qui empêche tout mouvement, alors le domaine de variation 
du mouvement est réduit à zéro et la particule perd son degré de liberté. Dans l’ exemple 
considéré ici (voir figure 1.7) la direction du mouvement est une droite (cartésienne). 
C’est un espace à une dimension géométrique correspondant à un degré de liberté phy- 
sique. La particule pourrait de ne pouvoir se déplacer que selon une courbe quelconque, 
disons la deuxième courbe de la figure 1.7. Encore une fois la particule n’a qu’un seul 
degré de liberté, une courbe étant un espace à une dimension, un seul nombre ou coor- 
donnée étant suffisant pour déterminer la position de tout point sur la courbe, par exemple 
la distance orientée (+ ou —) par rapport à une origine O quelconque. 


d x 
e . 
— 
x 


Figure 1.7 


On peut donc prendre pour règle que le nombre de degrés de liberté d’une particule 
est égal au nombre de coordonnées nécessaires et suffisantes pour déterminer la position 
de la particule. Compter le nombre nécessaire en général n’est pas difficile, un système 
physique comptant n particules pouvant toutes se déplacer dans un espace à D dimen- 
sions aura nD degrés de liberté même si ces particules sont en interaction à condition que 
ces interactions ne limitent pas à zéro les domaines de variation. Prenons par exemple 
deux particules ponctuelles, 1 et 2 dans un espace à deux dimensions (voir figure 1.8). 
Ce système compte 2 x 2 = 4 degrés de liberté. Pour décrire ces 4 degrés de liberté ou 
peut choisir les 4 coordonnées £1, y1, £2, Y2. On peut aussi choisir z1, y1, 0 et r, cette 
dernière coordonnée mesurant la distance entre les deux particules. À chaque fois, quatre 
coordonnées sont nécessaires et suffisantes pour décrire les directions selon lesquelles 
les composantes du système peuvent se déplacer, i.e. définir exactement la position des 
deux particules du système. Dans ce problème il existe des familles de solutions, cor- 
respondant à des conditions initiales spéciales, qui ont comme caractéristique, soit 0 = 
constante soit que r = constante et où il apparaît donc que le domaine de variation de 
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certaines coordonnées est réduit à zéro, semblant indiquer que le nombre de degrés de 
liberté est maintenant de moins de quatre. Il n’en est rien, le système continue d’avoir 
quatre degrés de liberté, un simple changement des conditions initiales demandera quatre 
coordonnées encore une fois pour décrire le mouvement. Le nombre de degrés de liberté 
ne se compte pas dans la solution mais est une propriété intrinsèque du système physique. 


Figure 1.8 


Supposons maintenant que le ressort soit remplacé par une tige rigide sans masse 
de longueur l (voir figure 1.9). Le domaine de variation de la distance entre les deux 
particules est réduit à zéro. Un degré de liberté vient de disparaître. En effet on peut 
écrire soit 


r=l dr=0 
soit 
y (z2 — 11)? + (go — y) =1 
alors 


dyj (z2 — £1)? + (y2 — y) = 0. 
Dans la première équation on lit directement que r est réglé à la valeur l. Il ne reste que 
le degrés de liberté décrits par £1, y1, 0. Dans la deuxième équation on lit qu’il existe un 
relation de dépendance entre quatre coordonnées (£1, Y1, £2, Y2). Algébriquement cela 
signifie que trois seulement des quatre coordonnés sont indépendantes. Ainsi donc un 
degré de liberté est décrit mathématiquement par une coordonnée indépendante. Cela si- 
gnifie que, physiquement, un degré de liberté correspond à une direction généralisée le 
long de laquelle le système peut se déplacer indépendamment des autres directions, i.e. 
en les gardant constantes. Clairement ici, si on varie z1, £2, et yı par exemple, alors y2 
n’est pas libre de prendre n’importe quelle valeur. y2 est contraint de prendre la valeur 
telle que 7 = l ci-dessus. Ce n’est pas un degré de liberté puisqu'il n’est pas indépen- 
dant des autres. Nous aurons à revenir sur la notion de degré de liberté. Notons ici que 
nous les comptons dans l’espace physique, en général l’espace à 3 dimensions dans le- 
quel se situe la mécanique classique (ou ses sous-espaces à 2 ou 1 dimensions). Il existe 
aujourd’hui des domaines d’études en physique, par exemple celui appelé systèmes dyna- 
miques, où on préfère travailler dans un espace de phase qui contient les vitesses en plus 
des coordonnées. Par exemple, l’espace de phase correspondant à notre espace physique 
habituel décrit, disons par les coordonnées x, y, et z, comprendra également les vitesses 
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Figure 1.9 
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t, y, et 2. C’est un espace à 6 dimensions et il est commun en système dynamique de 
compter coordonnées et vitesses comme étant des degrés de liberté. Comme nous le ver- 
rons la chose se justifie aisément mais nous garderons ici notre notion de degré de liberté 
défini dans l’espace physique seulement. Simple question de convention. 


———< }— 


sr 
2  FORMALISME DE LAGRANGE 


Le formalisme de Lagrange permet d’étudier une vaste gamme de problèmes en mé- 
canique. En ce sens il est équivalent au formalisme de Newton mais, il a sur ce dernier un 
certain nombre d’avantages. D'abord, il est fondé sur un principe théorique fondamental 
et élégant. Il utilise des quantités scalaires plutôt que vectorielles et, en ce sens, sa forme 
est indépendante des coordonnées utilisées. C’est également la porte d’entrée à une foule 
de méthode qui forment la base de la physique moderne en mécanique quantique et dans 
les théories de champs classiques et quantiques. 

Nous présenterons d’abord la méthode dans un cadre assez simple pour ensuite en 
souligner certaines limites d’application. L'intérêt et les avantages de ce formalisme de- 
viendront graduellement évident. 

Afin de souligner l’invariance de forme selon les types de coordonnées utilisées, nous 
les noterons q; et on les appelle souvent coordonnées généralisées. Elles sont absolument 
quelconques sauf pour les limitations que nous verrons dans la section sur les contraintes. 


2.1 Résultats d'expérience et principe de base 


Nous discutons ici d’une particule ponctuelle dont la position instantanée est donnée 
par les trois nombres notés {g;| i = 1,2,3} . Cette particule suit une trajectoire qui se 
développe avec le temps t et dont l’équation 
est le résultat recherché. Le long de la trajectoire, on définira les composantes de la vitesse 
généralisée {d;| i = 1,2, 3} définies par 

di = dat i = 1,2,3. (2.2) 
dt 

Notre expérience consiste en une source de particules (identiques) que nous nous si- 
tuons au point P, et en un détecteur que nous situons en Pz. A un temps noté {1 nous 
émettons une particule en P; (de coordonné q;(t1)). Nous ne nous intéressons qu’ aux 
particules détectées en P> à un temps tz tel que t2 — tı est le même pour toutes les expé- 
riences. Nous répétons l’expérience un bon nombre de fois. À priori il y a un nombre 
infini de trajectoires possibles pour les particules satisfaisant les paramètres de l’expé- 
rience: Co, C1, C2, C3, ... (voir figure 2.1). Pour les distinguer les unes des autres, uti- 
lisons un paramètre tel que la trajectoire C obéit aux équations 


a = (6) (2.3) 
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où, pour un À donné qi (t) £ qe? (t) pour a # a’ (deux trajectoires différentes). Ayant 
filmé l’ expérience, nous constatons que les particules ayant satisfait les paramètres de 
l'expérience ont toute utilisé la même trajectoire, disons Co. La nature semble donc 
préférer cette trajectoire et la choisit toujours. 


Figure 2.1 


La méthode de Lagrange compare les trajectoires possibles entre elles et nous donne 
un critère pour choisir la bonne. Pour ce faire nous calculerons (en principe) une quantité, 
notée S(a), qui caractérise la trajectoire 


tı 
S(a) = f L (aP 0, iH, t) dt 2.4) 
to 
où la fonction L, qui reste à déterminer, dépend des q;(t), des qi(t) et possiblement ex- 
plicitement de { lui-même. On aurait pu prévoir que L ait une dépendance en ÿ;— mais 
l expérience nous indique que ce n’est pas nécessaire. Ayant calculé (en principe) S(«) 
pour toutes les trajectoires nous décidons de la bonne en comparant les différentes va- 
leurs obtenues pour S. Pour pouvoir choisir un donné il faut que S prenne une valeur 
particulière en ce point (trop arbitraire) ou ait un comportement particulier. Le compor- 
tement le plus simple à identifier, c’est le point stationnaire, là où S est un extrémum. 
C’est le cas de a sur la figure 2.2, mais également de «1. Dans ce qui suit nous suppose- 
rons toujours qu’il s’agit de (bien que ce soit difficilement démontrable. Nous écririons 
donc que 
al, (2.5) 
da |z 
définit & et fixe ainsi la bonne trajectoire sur laquelle S prend la valeur extrême (mini- 
male) S(&). 

La quantité S s’appelle l’action et le principe énoncé ci-dessus est le principe de 
moindre action. C’est un principe variationnel, i.e. nous recherchons un point fixe de S 
tel que dS = 0. Aujourd’hui, on tend à baser toutes les lois de la physique sur de tels 
principes. 


2.2 Variation fonctionnelle et application du principe 
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Figure 2.2 
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Ci-dessus nous avons écrit 
dS = 0 (2.6) 
comme si la variation de S en était une au sens habituel, i.e. le long d’une trajectoire. 
Or, ce n’est pas le cas du tout, la variation est faite en comparant des trajectoires, i.e. en 
variant selon les fonction q (t) (voir figure 2.3). On notera de telles variations à l’aide 
du symbole plutôt que du symbole 6. La différence est très nette 


daf (E) = qi (t + dt) — gi) (#) 2D 
aa) = PO — ai A). 2.8) 
Sur une trajectoire donnée on connaît q;™ = A) et les vitesses a” (t) sont 


q; N 


t t+dt 


Figure 2.3 


fixées. Mais en comparant des trajectoires on constate sur la figure 2.4 qu'entre a et a’, 
ôq (t) est le même qu’en comparant a avec a”. Ceci n’est pas vrai des 6q(%)(#). Les 
variations des vitesses sont donc indépendantes des variations des coordonnées dans ce 
formalisme parce que nous comparons des trajectoires différentes. Ces variations étant à 
temps constant i.e. par exemple 


dP O — af E) = PO a O = qi (4) 2.9) 
mais , n 
TOR OE O a O (2.10) 
les variations en a et en temps t sont indépendantes et 
(9 = 5 (Lo) =L( ow) 
HO = 6 (F0) =$ (0). 1n 


Si nous calculons la différentielle ordinaire d’une fonction f(x, y), i.e. df, nous ob- 
tiendrons 
df df 


df (z, y) = dde + y” (2.12) 


si les variations dx et dy sont indépendantes. Le même type d’opération s’applique au 
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Figure 2.4 
calcul, par exemple de L 
. dL dL. 
SL dit) =) Ttt O gë (2.13) 


et il n’y a pas de terme ÎL ôt puisque ôt = 0, les variations étant à temps constant. 


Pour appliquer le principe de moindre action nous aurons à calculer 


to t2 
tı tı 


puisqu'on peut intervertir les variations en t et en a et puisque q; (t1) = 0 = q;(t2) étant 
donné que selon les paramètres de l’ expérience en t la particule est nécessairement en 
P» et en {2 elle est en P2. Ces deux points ne sont pas variés, toutes les trajectoires 
considérées devant les relier. 


Nous aurons donc 


t2 
0 = f ÔL(qi, di, t)dt 


tı 
t2 dL dL 
= — qit ð — | dt 
FEES 
i dL dL dqi 
+ | dt. 2.15 
Lae en 


Intégrons par parties le deuxième terme du crochet | ] 


t2 : t2 t2 
1 CRL -f n (à) jna (2.16) 
t tı tı dt \ di 


di dt Oi 
où le premier termde de droite est zéro puisque ôq; (t1) = ôqi(t2) = 0. Remplaçant nous 


II 


1 
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avons 


t2 dL d /dL 
an E cs ddt = 0. 2.17 
L Ela ala). Q17 


Pour aller plus loin nous devons faire l’ hypothèse que les q; sont indépendants les 
uns des autres. En termes physiques, cette indépendance des ôq; signifiera que les q; 
sont indépendants les uns des autres, i.e. qu’il n’existe aucune contrainte les reliant. Ils 
devront donc correspondre à des degrés de liberté physique du système. 

Posant donc que les q; sont indépendants et comme ils sont quelconques, la seule 


façon de satisfaire cette équation est que chaque terme dans | ] de (2.17) soit nul, i.e. 
dL d {dL 
ba —|—=0 2.18 
généralement écrit comme 
d (dL dL 
— |> |] -~ =0. 2.19 


Ce sont les fameuses équations d’ Euler-Lagrange. Nous posons qu’une fois solutionnées, 
elles définissent une trajectoire privilégiée 

qi = qi(t) (2.20) 
qui est identifiée à la trajectoire physique. 

Nous avons débuté en parlant d’une particule mais clairement, cela n’a eu aucun 
impact dans le développement de cette équation. Elle demeure valable pour un système 
à un nombre arbitraire, n, de degrés de liberté pourvu qu’ils ne soient pas contraints. 
Nous obtiendrons alors n équations pour à = 1, 2....n. De plus, rien n’a été dit sur les 
{qi} . Ils sont quelconques et mesurent des longueurs, des angles, des... La forme de 
l’équation n’est pas affectée par le choix des {qi} . 


Remarque 3 
On remarque ici qu’étant donné que les q; sont quelconques, ils n’ont pas nécessairement 
les mêmes dimensions. Là où dans l'équation de Newton 


F = ma (2.21) 
toutes les composantes de cette équation vectorielle ont une dimension de [M LT-?], il 
wen va pas de même des composantes de l’équation d’Euler-Lagrange. Elles n'auront 


dimension de forces que si qi a les dimensions de longueur L'approche Lagrangienne 
fait automatiquement la cuisine des dimensions. Elle est dimensionnellement homogène. 


2.3 La fonction L(qi, di, t) 


Il est évident, toute la validité de la méthode repose sur le choix ou la définition de L. 
Il devrait être également évident, étant donné que les équations d’Euler-Lagrange préten- 
dent résoudre le problème mécanique en ayant la trajectoire physique comme solution, 
que ces équations devraient correspondre aux équations de Newton. On peut de fait dé- 
montrer la forme de L à partir des équations de Newton. Nous en postulerons la forme 
et vérifierons le bien fondé de notre hypothèse. 
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Forces conservatrices 


On appelle une force conservatrice (sur une particule), une force F telle que V x F = 
0. Une telle force F(r) peut s’écrire alors 
F(r)=-VV(r) (2.22) 
où V (r) est appelé le potentiel ou l’énergie potentielle. 
On vérifie facilement alors qu’on peut écrire 


L=T-V (2.23) 
où T est l’énergie cinétique. 
Vérifions-le pour une particule soumise à une telle force et utilisons les coordonnées 
cartésiennes que nous noterons x; = (x, y, z). Alors 


1. 
T = 5 mij V = V(x;). (2.24) 
j 
Donc 1 
L=T-V= D gmi — V{x;). (2.25) 
j 
L équation ď’ Euler-Lagrange pour le degré de liberté x; (i fixé) demande que l’on calcule 
OL OV 
a 2.26 
OL d (ƏL 
— = mii —|—)|=mi . 2.27 
mb => (=) mä (2.27) 
L équation d’ Euler-Lagrange donne donc ici 
V 
miti + =0 (2.28) 
Ozi 
on Əv 
mti = —— = F;. (2.29) 
Oxi 


Il y a donc équavalence complète avec Newton. 

Dans l’ approche Lagrangienne, on apprend à raisonner à partir de concepts d’ énergie, 
potentielle et cinétique, au lieu de concepts de force. Les deux approches sont évidem- 
ment équivalentes physiquement, mais les énergies n’étant pas des quantités vectorielles, 
elles sont conceptuellement plus faciles à utiliser dans une vaste gamme de problèmes. 
En physique quantique par exemple, la notion de force n’a aucune signification mais les 
notions d’énergie demeurent valables. C’est une raison de plus pour se familiariser avec 
leur utilisation. De plus, la force au sens de Newton est une action instantanée à distance. 
En relativité, une telle chose est impossible. La notion de force est donc une création 
purement classique et macroscopique et contrairement à notre intuition, son intérêt est 
limité. 


Quelques exemples importants: 


Il est important de noter que nous n’avons identifié que quatre types d’interactions 
(forces) fondamentales dans la nature: gravitationnelle, électromagnétique, faible et 
forte. Les deux dernières étant purement quantiques, seules les deux premières se mani- 
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festent en physique classique. Or, la force gravitationnelle est du type en r7? et dérive 


donc d’un potentiel 
GMm 1 
Vaza = —— x ——. (2.30) 
r r 
Il en va de même de l’interaction coulombienne qui fait partie des interactions électro- 
magnétiques (nous reviendrons plus tard sur l’ensemble des forces électromagnétiques) 
a aussi une force en r7? et de ce fait dérive d’un potentiel 
€1e2 1 


— x ——. (2.31) 
4TEor r 


Vcoulomb = 


Un autre cas important est celui de la force harmonique (typiquement le ressort par- 
fait) qui est —kzx et dérive donc d’un potentiel 
kz? 


Viam. = =y: (2.32) 


Bien que n’étant pas une interaction fondamentale de la nature, elle joue fréquemment 
un rôle important dans les calculs. En effet dans des systèmes à géométrie un peu com- 
pliquée, l’énergie potentielle d’une particule peut prendre une allure assez quelconque 
comme sur la figure 2.5. Cependant au voisinage de x, correspondant à un extrémum 


Figure 2.5 


de V (x) on peut faire l’expansion 


V(x) & V(xo) + (x — zo) V’ (£o) + EZ Ed yr(so) +. (2.33) 


Le premier terme est une constante sans grand intérêt. Le deuxième terme est nul, V’ (xo). 
La première approximation non triviale est donc 
(x — xo)? V 


V(x) ~v ——_——— — 
2 dx? m 


k 


où 


: un nombre. (2.35) 


Ainsi pour bon nombre de matériaux, la première réponse à un (petit) déplacement hors 
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d’équilibre est harmonique. Ceci est d’une grande importance pratique. 


Forces non conservatrices 


Mathématiquement, peu de fonctions F dérivent d’un gradient 


F(r) £ -VV(r). (2.36) 

Il en est ainsi par exemple des forces de frictions que l’on écrit souvent empiriquement 
comme 

Fhict. (r) £ =k” (2.37) 


où typiquement n ~ 1 pour les basses vitesses (écoulement laminaire) et n ~ 2 pour 
des vitesses plus élevées (écoulement turbulent). La constante k dépend entre autre de la 
géométrie du problème et sa détermination est généralement empirique. 
Pour tenir compte de tels effets, il faut alors définir une force généralisée, de compo- 
santes Q; et notre équation ď’ Euler-Lagrange devient 
d (dL dL 
T (Z) — = Qi. (2.38) 
En général il faut être prudent dans la détermination de Q; puisque les composantes 
ne sont pas nécessairement cartésiennes et que les équations n’ont même pas toutes les 
mêmes dimensions! 
Il existe une exception notable et qui apparaît aujourd’hui comme extraordinairement 
importante. Nous en discuterons plus loin dans le cadre de l’invariance de jauge et nous 
verrons qu’elle correspond à l’ interaction électromagnétique complète. 


2.4 Coordonnées curvilignes 


On écrit communément 


1 
T= zov (2.39) 
où il est entendu que 
v=v.v=it. (2.40) 


Cette notation peut rapidement prêter à confusion. En effet, en coordonnées cartésiennes, 
il n’y a pas de problème 


r =(x,y, z) = x + y + kz (2.41) 
ii + y + ki (2.42) 
donc 
t= (à,9,2) (2.43) 
et 
vi? + + 22. (2.44) 


Cette simplicité vient du fait que x, y et z ont tous les trois des dimensions de longueur 
et que leurs axes sont fixes et orthogonaux. Qv’ arrive-t-il lorsqu’on passe à d’autres 
coordonnées? Prenons par exemple les coordonnées sphériques où 


r =(r,0,4) (2.45) 
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où 0 et y sont des angles (voir figure 2.6). Doit-on sans ambiguïtés définir 
i=(r,0, D)? (2.46) 
Deux problèmes surgissent ici, dus aux fait que (voir figure 2.7): 


i 


Figure 2.6 


1. r,0 et y n’ont pas les mêmes dimensions, 


2. leurs axes sont orthogonaux mais ne sont pas fixes. 


Les axes cartésiens (a) demeurent parallèles à eux-mêmes en différents points, ici 
P, et P2 alors qu’en (b) on voit que les axes du système sphérique sont en tous points 
perpendiculaires l’un à l’autre mais (P, ) n’est pas parallèle à (P: ), etc... 


En effet, si nous écrivons le rayon vecteur 


r —ir + jy + kz, (247) 
nos obtenons 
. dr + 7. +. 
nr = it + jy + kz (2.48) 


parce que i = j = k = 0. On a cette simplicité qu’en coordonnées cartésiennes. De fait, 
sachant que 


r sin 0 cos y 
= rsinĝsing (2.49) 


= rcosô 


8 
| 


on calcule 
à = tsinb cosy + rô cos 0 cosy — rý cos0 sin p 
ÿ = řsinĝsing + rô cosð siny — rÿsin 0 cos y (2.50) 
à = řcosô+rôsinð. 


De 
v=? Hy H. (2.51) 


P, 

y a 
La 

Figure 2.7 
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on obtient 3 
v=? +720 + r? sin? 0? (2.52) 
Les coordonnées étant r, 0 et p, clairement 
1 1 2 1 
T £ =mi? + =m + =mp?. (2.53) 
2 2 2 
mais bien _ ee d 
T = St + mi + A sin? 0%? (2.54) 
De façon générale on écrira alors 
T — 5 = gijdid; : coordonnées généralisées. (2.55) 
i,j 


Pour les coordonnées cartésiennes, on identifie, q; = (x, y, z) et 


1 0 0 
gy=ļ| 010), (2.56) 
0 0 1 
alors que pour les coordonnées sphériques, q; = (r, 0,4) 
1 0 0 
gj=| 0 7? 0 : (2.57) 


0 0 r?sin?8 
Ici, dans les deux cas, g;; est diagonal parce que les deux systèmes d’axes restent ortho- 
gonaux en tout point. Pour le cas sphériques, les g;; ne sont pas des constantes, mais des 
fonctions de la position qui tiennent compte simultanément du fait que 8 et p n’ont pas 
des dimensions de longueur et du fait que les axes F, 8 et @ varient en direction d’un point 
à l’autre de l’espace. gi; s’appelle la métrique (tenseur métrique) et il apparaît générale- 
ment dans la définition de l’élément de longueur souvent noté 


ds? = dr-dr. (2.58) 

En coordonnées cartésiennes 
ds? = (dx)? + (dy)? + (dz)? (2.59) 

alors qu’en coordonnées sphériques 

ds? = (dr)? + r? (d0)? + r? sin? 0 (dy) (2.60) 

et on écrit de façon générale 
ds? = DETTE (2.61) 

j 


C’est là la définition formelle de la métrique, gij, qui a une dépendance sur les coor- 
données (en général). Fondamentalement la métrique permet de définir la longueur dans 
un espace donné. On vérifie facilement ci-dessus que g;; est identique au g;; qui nous 

permet de définir sans ambiguïté l’ énergie cinétique T par 

m He 
T=} zusii. (2.62) 
i,j 

Tout ceci est très important pour obtenir les équations du mouvement, spécialement 


lorsque les coordonnées utilisées ne sont pas les coordonnées cartésiennes. En effet, dans 
le cas des coordonnées cartésiennes, l équation de Newton est F = ma avec les compo- 
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santes 


d 
là où x;—(x, y, z) pour i = 1, 2, 3. Du Lagrangien 
L= >. — V(xi), (2.64) 

les équations d’Euler-Lagrange nous donnent 

d OV 

— (mi) = -—. 2.65 

TA R T (2.65) 
Identifiant F = —VV, les deux équations sont identiques et Lagrange concorde avec 


Newton. 


En coordonnées sphériques, par contre, et l’équation de Newton pour 0 west pas 


d ; 
Z (mbi) = Fo. 2.66 
dt ( Í (apm 

Sachant que le Lagrangien sera 

L= z (r° + r20° + r? sin? 6?) _V(r,0,@) (2.67) 
l’équation d’Euler-Lagrange pour 0 nous donnera, avec 
E . 
2 = mr?0 (2.68) 
00 
d L - f 
. (5) = mr? + 2mri0 (2.69) 
et ƏL 
39 — mr? sin 0 cos b? (2.70) 
donc 
z E ak .2 OV 
mr 0 + 2mrr0 — mr” sin 0 cos 0 = —-—, (2.71) 
ce qui n’était pas à priori évident. On sait retrouver ce résultat à partir de l’équation de 


. . . 2 . . 
Newton si on fait attention dans le calcul de dz, Cependant la cuisine est relativement 
désagréable. La méthode Lagrangienne nous donne ” automatiquement” la bonne équa- 


tion. 
On remarquera qu’en divisant par mr?, l équation en 0 est 
1 OV 


T (2.72) 


MRE a 

Ë + “50 — sin 8 cos 4? = — 

r 

Le côté gauche est de la forme 
äi + X Tindig. (2.73) 

j,k 

C’est ce qui s’appelle la dérivée covariante par rapport au temps du vecteur vitesse de 
composante ġ;. Ici, si q; = (r, 0, p) pour i = i, 2, 3, équation en 0 correspond à i = 2 
et au lieu de Zii = Gi = 0, la bonne définition de la dérivée par rapport au temps, tenant 


compte des unités et du fait que les vecteurs unitaires varient d’un point à l’autre, donc 
dans le temps le long de la trajectoire nous avons des termes additionnels 


Ô + Dir + Tr + Dir + TA ôr + 12,00 
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+1 330p + Tai + T30 + Tab. (2.74) 
Tenant compte du fait que T° el j On identifie, pour les coordonnées sphériques 
Tî = 0, 
1 
r = I=- 
21 T , 
TÎis = IT} =0, (2.75) 
15: = 0, 
Ts = T2, 
Ti = sinl cosð, 


À partir de L on peut identifier les T1 jp et les r3 k de la même façon. Ce qui distingue 
les coordonnées cartésiennes, c’est que tous les I, = 0, c’est le seul système de coor- 
données pour lequel c’est vrai (et uniquement parce que l’espace considéré ici est plat, 
i.e. sa courbure est nulle). Ces facteurs géométriques, T° p, appelés symboles de Chris- 
toffel, jouent donc un rôle important. On peut les calculer par la formule 


gij | Ogre _ Og5r 

a j j 

2 2.76 
J -32 [6 [Se og dq Fo 
où g`! est la matrice inverse de g. On voit qu’ils sont entièrement déterminés par la 


métrique, g. Cette cuisine compliquée, la méthode Lagrangienne la fait automatiquement. 
Ce n’est pas le moindre de son intérêt! ! 


2.5 Les contraintes 


Il peut exister plusieurs types de contraintes, par exemple x = a signifie que le 
mouvement est gelé en x et qu’il est contraint de ne se faire que dans le plan yz passant 
par x = a. Il ne reste que deux degrés de liberté, y et z. On peut également avoir une 
contrainte du type 

ÿ = 4, (2.77) 
i.e. la vitesse selon y est contrainte d’avoir la valeur a. Cette équation s’intègre triviale- 
ment pour donner 


y = at + b. (2.78) 
Soit le Lagrangien (avant pa tenir compte des contraintes) 
L= T (a +4? +27) — V(x,y,2) (2.79) 


Si la contrainte est x = a e £ = 0, on devra écrire 
L= T p + 2) — V(x,y,2) (2.80) 


et nous n’aurons que deux i d’ Euler-Lagrange, une pour y et une pour z. 
Si la contrainte est ÿ = a cu y = at + b, on devra écrire 


L=— e +a? + 47) — V(x,at +b, 2) (2.81) 


1 Vector Analysis, M. Spiegel, Schaum. 
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et nous n’aurons que deux équations d’Euler-Lagrange, ici une pour x et une pour z. 
Notons que ces solutions seront paramétrisées par b s’il reste inconnu. 


De façon générale une contrainte s’écrit sans la forme 


> 0 


F(q, 02, dis doses) = 0 . (2.82) 
< 0 


Les cas d’inégalité correspondent à des contraintes non-holonomes. En fait on définit 
comme contraintes holonomes, les contraintes qui s’écrivent 


d 
fai, di, t) = grld =0 soit h(qi,t) =C (2.83) 


où h(qi, t) est une fonction quelconque des coordonnés (et du temps). On appelle non- 
holonomes celles qui n’obéissent pas à une telle relation, soit que 


2. ou Flai dit) < 0 ou I tds dt) > 0. 
Nous parlons de trajectoires, i.e. de l’existence de fonctions 


qq) — di=d(). (2.84) 
Par conséquent, pour une contrainte holonome 


d ðh ðh 
Flai dt) = gt) IË t e (2.85) 
De telles trajectoires satisfont 
h(qi,t) =C : une constante. (2.86) 


Dans les deux exemples vus précédemment les contraintes sont holonomes. Nous avions 
d’abord étudié x = a. Ici, nous aurons simplement h = x = a où C = a et 


ðh 
= — =ċł=0. 2.87 
f= 5? 2.87) 
Par contre, le deuxième cas étudié correspond à 
ÿ=a (2.88) 
et nous écrivons 
f=ġ-a=0 (2.89) 
ce qui même à 
h=y-at=C. (2.90) 


De façon générale, une contrainte holonome est intégrable au sens où on peut (même 
si c’est compliqué) l’écrire sous une forme permettant une substitution exacte dans le 
Lagrangien, faisant ainsi disparaître les degrés de liberté contraints. Physiquement on 
peut visualiser la contrainte comme étant due à une force extérieure telle que son effet 
impose au mouvement d’être contraint. Si cette force est indépendante des (i.e. la même 
pour) trajectoires possibles, alors la contrainte est holonome. Si cette force dépend de la 
trajectoire (raie d’une trajectoire à l’autre) alors la contrainte est non holonome. 
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Méthode des multiplicateurs de Lagrange 


Si un Lagrangien L dépend de degrés de liberté contraints, les équations d’Euler- 


Lagrange qu’on peut en déduire 

d (dL dL 

— | —]-—=0 (2.91) 
ne sont pas valides. Elles ne peuvent donc pas représenter nos équations de mouvement. 
Ce Lagrangien est inutile. Or, lorsque les contraintes sont non holonomes nous sommes 
en général incapable d’extraire exactement les degrés de liberté contraints du Lagran- 
gien. Même pour certaines contraintes holonomes, l’exercice peut être difficile. Il existe 
une méthode, dite des multiplicateurs de Lagrange, qui peut alors être utile. Nous la pré- 
sentons sans démonstration. 


Soit un Lagrangien, L(qi, ġi, t), à = 1, 2...n décrivant un système mécanique dont 
les trajectoires doivent obéir à une contrainte qu’on sait exprimer comme 


F(q;, dj, t) = 0. (2.92) 


On construit alors un Lagrangien auxiliaire, L’ 


L'=L+); (2.93) 
pour lequel on suppose que la contrainte est (temporairement) levée. Ceci étant, les n de- 
grés de liberté peuvent être considérés comme indépendants et les n équations de Euler- 


Lagrange 
d {dL’ _ dd 
— | — =0 :1=1,2,: 2.94 
sont valides. En principe on peut résoudre pour obtenir les équations de la trajectoire 
qi = q(t, À) (2.95) 
qui seront paramétrisées par À puisque L’ en dépend. On peut en calculer les 
d 
hi = —qi = di(t, À). 2.96 
di = gt = dit, À) (2.96) 
On remplace alors dans l’équation de contrainte 
F(Q(E, À), GE, à), t) = 0 (2.97) 


qui permet de calculer la valeur de À 
À (2.98) 


permettant à la contrainte d’être satisfaite. On remplace alors cette valeur de À = À dans 
les équations de la trajectoire pour obtenir les équations de la trajectoire contrainte 


qi = qilt, A) ; i= 1,2,...n. (2.99) 


Pour simple qu’elle soit en apparence, cette méthode n’est pas triviale d’application. 
En effet, on doit prévoir de 


falt, E i(t, à), t) = (2.100) 
que la solution dépende de t i.e. A = À = A(t) n généralement de t. Or, si 
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l’équation de contrainte dépend des ġ; alors 
dL’ dL of 
dåi då; dd 


d (dL'N\ d {dL d (Əf\ :Əf 
q (T) “T (T) B (5) TAi Rak 


et nous voyons apparaître non plus seulement mais aussi inconnu. C’est d’ailleurs tou- 
jours le cas pour les contraintes non-holonomes. 


(2.101) 


et 


Nous ne pousserons pas plus loin la présentation de cette méthode qui nous mènerait 
à des divergences considérables. Pour ceux qui sont intéressés on peut consulter les livres 
de Goldstein ou de Saletan et Cramer par exemple. 


2.6 Invariance de jauge 


On appelle transformation de jauge une transformation de L en L’ 
; à d 
L' (qi, di t) = Lai, dit) + E) (2.103) 
où F est une fonction des q; et de t (appelée génératrice de la transformation) et 


dF dF. OF 
I~ X ii + — (2.104) 


Remplaçant dans la définition de l’action, S devient S” 


t2 
S' = f L' (qi, &,t)dt 


ti 


t2 t2 d 


ti tı 


= S+F(2)-F(1) (2.105) 
où F(1) = F(qi(t1),t1) donc F (2) = ôF (1) = 0 puisque ôq; (t1) = ôqi(t2) = 0. 
Ainsi 
ôS’ = ôS. (2.106) 


Or comme la physique est déterminée par l’extrémisation de S, ou de S”, rien n’est 
changé ici. La physique sera inchangée, les trajectoires seront les mêmes. 
On constate cependant que le passage de L à L’ ne laisse pas la forme du Lagrangien 
inchangée. En effet, supposons que 
L=T—V(qi,t). (2.107) 


Alors une transformation de jauge générée par F(q;,t) nous donnera 


dF(qi,t). , OF(q:t) 
L'=T-V(g LG, + 2, 2.108 
Viat+2 de ËtT > (2.108) 
Puisque est une fonction de q; et de t, appelons 
OF(q;,t 
V'qist) = V(qt) — 2L (2.109) 
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Ainsi dF( 
=T- V'(qi,t EA (0) ds. (2.110) 


Le dernier terme fait que la forme de L n’est pas RA 
Opérons une deuxième transformation de jauge, générée par la fonction G(q;, t). 
Nous obtiendrions de L’ un nouveau Lagrangien, 5 


pot) D -F +G) (2.111) 


où 


V'(ait) = Vlant) — Š [F(q t) + Glai 0). eu? 


On voit donc que L’ est invariant de forme sous une transformation de jauge qui laisse 
la physique inchangée. Aujourd’hui on a admis le principe théorique qu’il s’agit de la 
forme la plus générale que peut prendre un Lagrangien, i.e. que les seules interactions 
possibles sont des interaction de jauge. C’est cette philosophie qui a permis l’unification 
de trois des quatre interactions fondamentales en théorie du champ. 

Le terme en di, i.e. ÿ;, Eii est la forme q-VF , i.e. le produit scalaire entre le 
vecteur et un champ vectoriel (local) que la transformation de jauge nous donne comme 
étant le gradient de F, VF. 

Supposons maintenant que notre Lagrangien L s’écrive 


où A (qi, t) est un vecteur quelconque, et non un gradient. Alors une transformation de 
jauge L z L’ donnera 


L' = L- V' (qi, t) +A A' (qi, t) (2.114 
où 
F(o; 
V'(q,t) = E (ant) (2.115) 
A'(qt) = Alq, t) + VE(q;t). (2.116) 


La forme du Lagrangien est clairement restée la même et nous savons que la physique (la 
trajectoire) n’est pas affectée par la transformation de jauge. En physique moderne, on 
adopte aujourd’hui une approche basée l’axiome suivant: la nature est telle qu’ observée, 
invariante de jauge (interaction électromagnétique). Nous devons donc développer un 
formalisme physique qui respecte cet aspect de la nature et qui soit invariant de jauge. 
En mécanique classique, cela signifie que le Lagrangien le plus général que l’on peut 
écrire à priori devra être invariant de forme sous une transformation, i.e. devra être de la 
forme 
L=T-—V(qt) + àA(q,t) (2.117) 
où V et A sont des champs locaux, scalaire et vectoriel respectivement. La conclusion 
qui s’impose est que les seules interactions permises par la nature sont celles décrites par 
ce Lagrangien. Il nous reste donc à vérifier quel type d’interaction existe dans la nature, 
au niveau classique, sur la base de cet axiome d’invariance de jauge. 
Typiquement donc un interaction invariante de jauge dépendra des vitesses puisque 
L ~ à'A(q,t), donc la force dépendra des vitesses. Clairement cette force n’est pas 
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conservatrice au sens vu dans le chapitre précédent, néanmoins de telles forces trouvent 
leur place dans le formalisme Lagrangien. Examinons le type de forces qui émergent de 


-= Ta? — V'(x,t) + x A(x.t) 


= SD -Vt Y yA t) (2.118) 
7 7 


qui restera invariant de forme lors d’une transformation de jauge même dans le cas géné- 
ral où 


A 4# VF. (2.119) 
L équation ď’ Euler-Lagrange pour la composante x; demande que l’on calcule 
ðL OV . Ô 
ne pe D (x,t) (2.120) 
= = md; + Á;(x,t) (2.121) 
d (ƏL o 
n (=)= mä; + 2 _ (ct) + Ait) (2.122) 
puisque, sur une trajectoire x; = s0. 
d = . ðf əf 
gret To D (2.123) 
Au total nous avons donc 
0À; 04: i 
i 7 2.124 
IR a ET Z p2 = GER 
donc 
SS ƏV (x,t)  OA;(x,t) . (0A; ƏAj 
më; = (ET +) + 2 an (2.125) 
C’est la composante x; de l’équation vectorielle 
. OA(X,t) | 
mx = — | VV (x,t) + rc +x x (V x A). (2.126) 


On sait qu’en électromagnétisme les champs électrique et magnétique peuvent être ob- 
tenus des potentiels scalaire et vecteur Väect €t Aélect 


E(x,t) = -VVasa (x,t) — Ps CE D (2.127) 


B(x,t) = V x Adea(x,t). (2.128) 


On sait également qu’une particule de charge e placée dans des champs E et B est sou- 
mise à la force de Lorentz 

FLorentz = e(E+x x B). (2.129) 
On peut donc aisément identifier V avec eVgect et A avec e A4. et conclure que l’in- 
variance de jauge du Lagrangien qui nous a permis de poser la forme la plus générale 
possible pour L nous mène directement à l’interaction électromagnétique. C’est un ré- 
sultat remarquable. 
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L'électromagnétisme possède également son invariance de jauge, i.e. que les champs 
physiques E et B sont invariants si on change simultanément 


0F 
Vétect = Vétect + ot (2.130) 
Acea > Aaea + VF. (2.131) 


Cette invariance de jauge électromagnétique est identique à l’invariance de jauge La- 
grangienne. 
Plus haut, nous avons identifié la partie de l’interaction de jauge qui dépend des vi- 
tesses comme étant de la forme 
X-A (2.132) 
par analogie avec le terme 


oF 
5 ii — =XVF. (2.133) 
Si A = VF i.e. si A est le gradient d’une fonction scalaire alors le Lagrangien 
T-V(xt)+xkAZ=T-VRi)+xkVF (2.134) 


décrit, par invariance de jauge, la même physique que le Lagrangien L qui apparaît au 
début et par conséquent il n’y a pas ici d’interaction nouvelle. Il n’y aura interaction 
nouvelle que si 

AZVPF, (2.135) 
i.e. il n’y aura interaction nouvelle ou de jauge que si A n’est pas le gradient d’une fonc- 
tion scalaire. De fait physiquement, en électromagnétisme, un potentiel vecteur qui n’est 
que le gradient d’une fonction scalaire ne génère las de champs. Il est évident que le cas 
particulier À = 0 est possible. Il permet de couvrir les interactions à potentiel habituel 
i.e. V (qi, t), ce qui permet les interactions électromagnétiques et gravitationnelles par 
exemple. 


2.7 Quelques caractéristiques, propriétés, limites... 


1. On ne saurait trop insister sur l’indépendance des coordonnées généralisées, les q;, 
qui décrivent les degrés de liberté physiquement indépendants. Si cette condition 
n’est pas satisfaite en écrivant le Lagrangien, celui-ci n’est pas valide et les équa- 
tions d’Euler-Lagrange qui en découlent non plus. Les trajectoires, solutions de ces 
équations n’ont rien de physique. 


2. En mécanique classique non relativiste, pour chaque vrai degré de liberté du système, 
qi, le Lagrangien contient un terme en ġ?. Le Lagrangien peut également dépendre 
linéairement de ġ;et sa dépendance en q; est quelconque. La dépendance en q? est 
nécessaire pour garantir que l’équation d’Euler-Lagrange sera en g;. Depuis Newton, 
on sait que la connaissance des deuxièmes taux de variation (d;) des q; est nécessaire 
et suffisante pour déterminer l’historique du système. La dépendance en g; apparaît 
avec les potentiels de jauge (potentiel vecteur) discuté à la section précédente. La 
dépendance en q; est quelconque. Elle dépend du système de coordonnées et des 
interactions. 
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3. Coordonnée cyclique: Une coordonnée q; (j fixé) est cyclique si elle n’apparaît pas 
dans le Lagrangien alors même que ce dernier dépend de ġ;. De l’équation d’Euler- 
Lagrange pour ce degré de liberté 


d (ƏL OL 
— |= |] -——=0 2.136 
di (3) Da: Mt 
il ne reste alors que 
d (ƏL OL 
— | — | = i ; 2.137 
k ( z) 0 puisque a (2.137) 
Formellement la solution est 
L 
al = T; = constante (2.138) 
Oi 


une telle équation est généralement plus simple à résoudre que l’équation d’Euler- 
Lagrange complète. 


4. Le Lagrangien L = T — V est structuré comme les énergies cinétique, T, et poten- 
tielle V. Il en partage plusieurs propriétés, en particulier l’additivité. Si Li et Lə sont 
des Lagrangiens de deux systèmes physiques indépendants, alors le Lagrangien du 
système physique constitué de l’union des deux précédents systèmes est 


L= Lı + Lə. (2.139) 
Si les deux systèmes interagissent alors le Lagrangien total sera 
L = Lı + Lə — V (1,2). (2.140) 


Peuvent jouer le rôle de sous-systèmes des particules différentes ou une même parti- 
cule à qui on octroie des degrés de liberté additionnels. 


Exemple 2.1 


Si Li = miii — Vı (rı) décrit le mouvement de la particule 1 et si La = maii — V2(r2) 
décrit le mouvement de la particule 2 alors le système physique constitué des deux particules sans 
interaction est 


L = L+L 


mu .2 
= = lita 
2 2 
Sienpl t aux d ticules d'interagir vi F; = -702 = l 
i en plus on permet aux deux particules d'interagir via une force Fij = or- = force sur la 
particule à due à la particule j, alors le Lagrangien est 


L = I1+Li—V(1,2) 
= L+Li—-V(r,r2) 
z a Vin 1) 
Va (re) — V2 (r1, r2). (2.142) 


M2. 2 


i2 — Vi (rı) — V2(r2) (2.141) 


2 


— r 
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3 APPLICATIONS ET PRO- 
PRIÉTÉS 


3.1 Cas simples en mécanique 


Particule dans un champ gravitationnel 


Une particule de masse m dans le champ gravitationnel près de la surface a une énergie 
potentielle V = mgz où z mesure sa hauteur et g est l’accélération due à la gravité (voir 
figure 3.1). Son énergie cinétique est 
\ 
x ou y 
Figure 3.1 
m,. ; í 
T= z +9 +27) (3.1) 
donc Ha 
L= ra +9? + 2°) — mgz. (3.2) 


Nous aurons trois équations d’Euler-Lagrange, celles pour x et y étant identiques. Voyons 
celle en x. On constate que 2L = 0. Dans un tel cas, on dit de x que c’est une variable 
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cyclique. L équation d’Euler-Lagrange pour x se limite donc à 


d (ƏL 

a Re) ET 3.3 

dt (5) a 
ou 2L =constante (d'intégration). De bL = më = C, on tire 

mz =ct+a 8.4 


où les constantes c et a sont déterminées par les conditions du problème. De la même 
façon 


my = ct+ a (3.5) 
Pour z nous avons 
OL 
Zo åA = 3.6 
5z mg (3.6) 
OL ; d {OL : 
alors 
mž+mg=0 ou ž= -g (3.8) 
donc 
gt? (14 1 


où c” et a” sont déterminées par les conditions du problème. 


Particule suspendue à un ressort 


© 


Une particule de masse m est suspendue à un ressort de constante k dans le champ gra- 
vitationnel près de la surface de la terre (voir figure 3.2). Son énergie potentielle 


m 
Figure 3.2 
est k 
V= G — z0)? + mgz (3.10) 
où 2ọ est la longueur au repos du ressort. Le mouvement n’étant que vertical, 
M 9 
T = —ż (3.11) 
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et 
k 
t= ur - (2 — 2) - mgz (3.12) 
donc 
E = —k(z— zo) — mg (3.13) 
d (ƏL 
—|=] = 2 .14 
E ae 
l’équation d’Euler-Lagrange (il n’y a qu’un seul degré de liberté) est 
mž + k(z — 20) + mg = 0 (3.15) 
m2 + kz = kzo — mg. (3.16) 
Posons z = 2} + Zp où 
mžn + kzn = 0 (3,17) 
alors (pour ms? + k = 0 a -E s = iw où w = È) la solution 
générale s’écrit , | 
zn(t) = A'e™* + B'e™* = Asin (wt + ô) (3.18) 


avec 


w= Je (3.19) 


m 
et À et ô qui sont des constantes d’intégration devant être déterminées par les conditions 
initiales. Le terme non homogène étant constant, on n’est pas surpris de trouver 


Zp = constante = C (3.20) 
z(t) = zat) + C (3.21) 
où Cest une constante. Remplaçant, avec 
2 = zp + 0 = —-w?A sin (wt + ô) (3.22) 
on obtient 
—mu? À sin (wt + ô) — kA sin (wt + ô) +kC = kzo — mg (3.23) 
a 
kO =kzxa-mg = C=%- (3.24) 
z(t) = Asin (wt + ô) + zo — Ts, (3.25) 
Ainsi, au repos ou À = 0 nous avons z = zọ — 4#, i.e. le champ gravitationnel cause 
un étirement du ressort (vers le bas) d’une longueur ma 


Particule suspendue au haut d’une tige rigide 


© 


Une particule est suspendue au haut d’une tige rigide sans masse et se déplace dans le 
plan xy (voir figure 3.3). Cependant la tige rigide constitue une contrainte telle que la 
particule ne se déplace que sur la courbe C qui n’a qu’une dimension. Le système ne 
possède donc qu’un seul degré de liberté. Dû au champ gravitationnel, l’énergie potentiel 
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Figure 3.3 
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est V = mgy. Géométriquement 


x =lsing, y= —lcos g (3.26) 
donc 
t =lþpcosy, y=lpsinys (3.27) 
donc yi m ni 
T = zË + ÿ) = zU o (sin p + cos? p) = a (3.28) 
et 
V = —mgl cos y (3.29) 
et ainsi 2 
Le Ty + mgl cos p. (3.30) 


C’est un Lagrangien pour un système à un seul degré de liberté, ici y, comme il se doit. Un 
choix intelligent de coordonnées généralisées, lorsqu’il y a contrainte, consiste à choisir 
les degrés de liberté contraints (ils ne sont plus des degrés de liberté alors) comme faisant 
partie des coordonnées généralisées. Ainsi dans cet exemple, le mouvement dans le plan 
xy peut être décrit en coordonnées polaires (r, p) où r = l est précisément l’équation 
de contrainte ici. C’est le choix que nous faisons, ce qui ne laisse que le degré de liberté 
décrit par y. Ce degré de liberté est un angle et n’a pas de dimension. Il n’y a qu’une 
seule équation ď’ Euler-Lagrange 


L 

—- = —mglsiny, (3.31) 
OL 
— = m°p (3.32) 
op 

d (ƏL 2 

— | =] = ml ọ . 
=> J (S=) mo (3.33) 

donc l’équation d’Euler-Lagrange se lit 

m + mglsinç = 0 (3.34) 


ou + Ẹsin y = 0, après division par ml?. La solution n’est pas triviale et fait appel à 
des intégrales elliptiques. Cependant si ọ reste petit, alors 


3 5 
npp- À +E + (3.35) 
et ne retenir que 
siny © p (3.36) 
et l’équation devient 
p+ = = (3.37) 
qui a comme solution 
p(t) = Asin(wt + ô) (3.38) 


avec w = V7 , À et ô étant des constantes déterminées par les conditions initiales du 
problèmes. C’est le fameux problème du pendule plan qui a longtemps servi de référence 
pour mesurer le temps. 
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Remarque 4 
Dans les deux derniers exemples, le mouvement est harmonique”. C’est le cas à chaque 
fois que l'équation du mouvement est du type 


ü + w?u = 0 ü=-wu ; w? >0 (3.39) 
qui a comme solution 
u(t) = Asin(wt + 6) = A' cos(wt + 0) (3.40) 
ou 
u(t) = Bsinwt + B'coswt (3.41) 


où w est la fréquence (angulaire) du mouvement. u(t) revient au même point à chaque 
fois que 

wt=2nr : nentier (3.42) 
onart = © = 1 où T est la période, v la fréquence du mouvement et w = 27v est la 


w 
fréquence angulaire. 


On notera également que dans les deux premiers exemples, la coordonnée utilisée 
pour décrire le degré de liberté a les dimensions de longueur, ainsi les équations d’Euler- 
Lagrange ont des dimensions de force, comme l’ équation de Newton. Il n’en va pas de 
même de le dernier exemple où la variable n’a pas de dimension (un angle). L équation 
d’ Euler-Lagrange est l’ équation du mouvement même si elle n’a pas des dimensions de 
force. 


Pendule plan suspendu par un ressort de masse nulle 


© 


Soit un pendule plan dans lequel la tige rigide est remplacée par un ressort de masse 
nulle, en fait négligeable (voir figure 3.4). Le mouvement étant dans le plan xy on attend 
deux degrés de liberté. Puisque la tige n’est pas rigide le mouvement n’est pas contraint 
à une trajectoire et on conserve deux degrés de liberté. On pourrait conserver x et 


Figure 3.4 


y pour les décrire mais ces coordonnées ne collent pas très bien avec la géométrie de 
l’objet. Clairement u et y collent mieux à cette géométrie où y est l’angle du pendule 
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avec la verticale et u sa longueur. On obtient facilement 


T—=USINYP, Y= —u cosy 
donc 
t = using +upcoswy, 
= —ùcosy + uy sin y 
et 
T = Ba +i) 


= (2 sin? p + 2uUp Sin p COS Y + ug? cos? y 


2 


+ù? cos? y — Quüÿ sin y cos y + uùġ? sin? p) 


m,. : 
= 7” + u2p?). 
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(3.43) 


(3.44) 
(3.45) 


(3.46) 


L'énergie cinétique a un terme en ù? et un en g?, ce qui est correct dans un Lagrangien 


destiné à décrire un système physique à deux degrés de liberté . 


L'énergie potentielle est 


V 


K 
mgy + = (u — uo)? 


K 2 
=  —mgu cos p + z“ — uo) 


et finalement 


K 
L= Ta +u? h?) + mgu cos  — z“ — uw)? 


(3.47) 


(3.48) 


où ug est la longueur au repos du ressort-tige. Ici nous avons une coordonnée qui a des 
dimensions de longueur (u) et une qui n’en a pas (4 ) puisque c’est un angle. Néanmoins, 


nos équations ď’ Euler-Lagrange seront de la même forme 


a (ƏL) ƏL a 
dt \ Où ðu 


dt \ 3p op 
Dans (3.49), 
ðL 
PT T ug? + mg cosy — K (u — uo) 
d fOL\ ə 
dt \au) 7 
et donc i 
mü? — SUP — mg cosy + K(u — uo) = 0 


qui a des dimensions [MLT?], i.e. de force. Pour (3.50) 


OL ; 

— = —mgusing 

dp 
d OL .: d 2e) | e 2e 
d (5) = (mu g) = 2muug + mu 


(3.49) 


(3.50) 


(3.51) 


(3.52) 


(3.53) 


(3.54) 


(3.55) 
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ce qui donne 
mu? + 2muü + mgu sin g = 0 (3.56) 
qui a des dimensions [MLT?], qui ne sont pas des dimensions de force 
La solution de ces deux équations couplées n’est pas triviale. Cependant leur obten- 
tion, par une méthode raisonnablement simple constitue déjà en soi un résultat intéres- 
sant. 


3.2 Exemples non mécaniques 


Principe de Fermat 


© 


On peut baser toute l’optique géométrique sur le principe de Fermat qui, remarquable- 
ment, est un principe variationnel. Les trajectoires des rayons lumineux à des équations 
d’ Euler-Lagrange. 
Énoncé: Entre deux points, 1 et 2, le rayon lumineux suit la trajectoire qui prend le 
moins de temps. 
Si n est l'indice de réfraction, la vitesse du rayon lumineux est c/n . Appelons ds 
T élément de longueur de la trajectoire, alors le temps requis pour parcourir ds est dt = 
ds/v. Entre les points 1 et 2 le temps requis sera T 
2 2 
T = G = 7 nds (3.57) 
1 v c Ji 
où 7 peut varier d’un point à l’autre comme dans une fibre optique par exemple. 
Pour simplifier limitons-nous à un système à deux dimensions, (x, y) donc 


ds = y dx? + dy?. (3.58) 


La trajectoire du rayon est une équation du type y = y(x). (on aurait pu choisir x = 
æ(y)). Écrivons donc 


du\? 
ds = dz4|1 + (4) = de V1 + (3.59) 
où y =% et ici x est considéré comme le paramètre et y une variable. Nous aurons donc 


1 2 2 
Te [nn ir P= | Loiad 3.60) 
1 1 
Le Lagrangien est ici 


L(y, ġ, 2) = Lyle, y)v 1+9. (3.61) 

Insistons sur le fait que le problème est semblable à un problème de mécanique à un 

degré de liberté, décrit par y. Ici x est le paramètre , ne décrit pas un degré de liberté 

et joue le rôle joué généralement par le temps en mécanique. Ainsi ce qui joue le rôle 

de la vitesse (un degré de liberté donc une vitesse) est ý =% i.e. la dérivée totale de la 
coordonnée y par rapport au paramètre x. 
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On cherche à minimiser T = Je Ldzx entre deux points fixes en comparant toutes 
les trajectoires y(x), qui les relient. On calcule donc 


ôT. (3.62) 
Le résultat est connu, c’est l’ équation d’ Euler-Lagrange pour le degré de liberté y (ici le 
seul avec le paramètre x) 
d (ƏL oL 
— | — | -——=0. 3.63 
dx (5) dy Fe 
Nous calculons donc, avec L défini ci-dessus 
L 1 
3 Z 2 iy 
ôy c 2/1 +9? 
ny 
— (3.64) 
cy/1 + 7%? 


homme 
dt \ dy cy/1+ÿ cÿ/1+3ÿ2|[0y0x ðr 
nyf 1 ar 

+8 ( :) (1+%°) 2yÿ 
c1+g  cy1+3ÿ2 LOy ox 

Eis 

L (3.65) 
c(1 +42)? 


et 
ƏL 1 an 


ðy cyi F g? ðy 


a a 
cy 1+9? cy1+ÿ Lôy” ôx 
nÿ°ÿ 1 ôn 


1 f) 
c(i+g)?  evT+ ôy 


(3.66) 


ce qui donne 


(3.67) 
et puisque c4/1 + ÿ? n’est jamais nul, on peut simplifier en 


a waon à à n Ô 
Cm + (+0) = 0 (3.68) 


3.3 Problème à deux corps 


C’est le système physique fermé le plus simple qui existe. Deux particules, de masses 
ma et m2, dont les positions instantanées sont r et ro interagissent via un potentiel 


V(r, r2) = V(rı E r2) (3.69) 
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pour respecter l’homogénéité de l’espace. Ainsi leur Lagrangien s’écrira (le Lagrangien 


est additif) 
Mi. M2. 
Le 42 V(r1 —r2). 


(3.70) 


Opérons un changement de coordonnés définissant la coordonnée relative r = rı — r2 


et la coordonnée du centre de masse (voir figure 3.5) 
MT + Mara 


R 
mı + Mo 
ainsi F 
rı= Rir 
mi Aus m2 
r= R-——r 
M1 + M2 
donc 
=R+ À 
mı pa m2 
i=- — à 
| mı + M2 
Remplaçant dans le Lagrangien nous obtenons 
À 
£ 
mi 
qxr 
TE a Sj 
# R Pa i 
/ Ts 
J f r, y - 
© 
x P 
P 
# 
Figure 3.5 
M: m. 
L= =R? + >r? — V(r) 
2 2 
où 
M = mi+ma2 = masse totale du système 
m ; 5 
= —— = masse réduite du système 
M1 + M2 


Le Lagrangien se décompose en deux éléments qui ne sont pas reliés 
L = Lem + La. 


Lou est simplement l’énergie cinétique globale du système et 


M 
Lem = TR 


(3.71) 


(3.72) 
(3.73) 


(3.74) 
(3.75) 


(3.76) 


(3.77) 
(3.78) 


(3.79) 


(3.80) 
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puisque 2L 0 2k C = constante 
MR=C. (3.81) 
Le C.M. se déplace à vitesse constante. La deuxième partie, relative, est 
Lea. = Ti? — V(r) (3.82) 


apparaît comme le Lagrangien d’une particule de masse m et de position r. Aucune des 
deux particules n’a la masse m et r ne donne la position d’aucune des deux particules. 
L décrit le mouvement relatif entre les ceux particules en le réduisant à un problème 
d’une seule particule (fictive), de masse m et de position r. Parce qu’il peut se réduire de 
cette façon à un problème à un corps, le problème à deux corps peut avoir une solution 
analytique. 

Le problème à N corps, ou N > 2, n’a pas de solution analytique. 


3.4 Le potentiel central 


Nous étudions ici un problème à un corps qui peut être le problème relatif d’un 
système à deux corps qui est aussi assimilable à celui d’une particule soumise à une 
force centrée à l’origine. Le Lagrangien est de la forme 


L= TE -V(). (3.83) 


Dans bon nombre de cas, l’interaction ne dépendra que de la distance, soit entre les (deux) 
corps, soit entre le corps étudié et le point d’origine de la force. On a alors V (r) = V (r) 
et la force est dans la direction r. C’est le potentiel central. 


Pa 


A 


w 
Figure 3.6 


Physiquement, si le problème de base est un problème à deux corps (voir figure 3.6) 
alors la force est purement dans la direction de la droite les reliant. C’est le cas de l’inter- 
action gravitationnelle entre deux corps massifs ainsi que le l’interaction coulombienne 
entre deux corps chargés. 

Puisqu’ici la force est purement radicale (aucune composante 0 et 4) le torque r x F 
s’exerçant sur la particule est identiquement nul et par conséquent le moment cinétique 
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(voir figure 3.7), 


l=mrxr (3.84) 
est une constante du mouvement. En effet 
l=mixi +mrxé=0+rxF (3.85) 
puisque r et F sont colinéaires. La conséquence physique est que le mouvement est 
À 
Zz 
o e 
“A, 
5 y 
Pa j z 
K n s ? 
P 
Figure 3.7 


dans un plan perpendiculaire à l puisque 1 = constant signifie constant en grandeur et en 
direction. Choisissons ce plan comme étant le plan xOy, i.e. le plan 0 = + donc 0 = 0. 
Le Lagrangien se réduira (avec sin 0 = sin 5 = 1) à 


L= Zo +r?) — V(r) (3.86) 


comme en coordonnées polaires. Immédiatement, on constate que y est une variable 
cyclique et donc 

ðL 

~i 

dp 
et donc a = mr? = l, une constante même si r = r(t) et ġ = ġ(t). 


(3.87) 


On vérifie trivialement que cette constante l est précisément la longueur du moment 
cinétique qui pointe ici selon l’ axe Oz. Comme est une constante du mouvement sa valeur 
est fixée par les conditions initiales. L'équation en r se calcule aussi: 


ðL 


OL d (ƏL 
F ~ m? = T (5) = mr (3.89) 
donc av 
mr — mro? + = (3.90) 
et toujours en @ : mr?ġ = l. 
De l’équation en ç on tire 
| l 
p = — (3.91) 
mr 
que l’on remplace dans l’équation en r 
È OV 
mr — =0 (3.92) 
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ou 
© &V Ê _ aValr) 
2 
2mr2° 


où 


Var(r) = V(r) + (3.94) 


Tout se passe donc en r comme dans une équation à la Newton pour un système à un 
degré de liberté 
: OVa(r) 


mi = x - Far). (3.95) 


Ce potentiel efficace Var(r) est constitué du potentiel original V (r) plus qui repré- 
sente une répulsion centrifuge: un corps qui tourne par rapport à l’origine O est effec- 
tivement repoussé de l’origine (il ne peut pas l’atteindre) et plus il tourne i.e. plus l est 
grand, plus il est repoussé. L exemple de la figure 3.8 est pour V = — K (gravitationnel 
ou électrostatique). On constate dans ce cas que 


vV 


Figure 3.8 
K P 
Ve = —— + — 3.96 
(r) ro 2mr? (3.20) 
a un extremum Vọ en ro qui obéit à 
OVe(r) K P 
ZE = — .97 
ðr |, +? mr? 3 Pon 
donc à 
T E RS (3.98) 
, 0 Km RE 
Puisque 
OV. 
Wah cp es pt 6 (3.99) 
Or | o 


ce qui correspond à r = ro donc À = 0 donc ë? = 0. C’est un point stationnaire qui 
correspond physiquement à une orbite circulaire (seul & varie) avec une vitesse angulaire 
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constante donnée par 


D = —; = constante (3.100) 
mr 
l 
=> t) = —;t | 3.101 
p(t) n2 + Po ( ) 


0 
Ici, les conditions initiales ont fixé r = ro, E = Vo, la valeur de l et de celle de Yg. Nous 
continuons d’avoir 0 = ? donc 0 = 0. 


Pour chaque intervalle de temps, 7, pour lequel y(t) augmente de 27 , nous complé- 
tons une orbite, donc 


K2 2 [3 3 : 
TT 2r = -o 7 T = période. (3.102) 
La fréquence est v = 1 = m et la fréquence angulaire de l’ orbite circulaire est 
K2 
D T. (3.103) 


Les conditions initiales pour que l’orbite soit précisément circulaire sont relative- 
ment peu probables. Etudions la situation lorsque l’orbite se dégage légèrement de ce 
cas particulier, i.e. E > V5 


OVar(r) + (r == ro)? d2Var(r) 


Varr(r) SVeælr) +(r — ro) = ——— (3.104) 
ps Or HS 2! Or2 " 
constante N— u 
=0 


Si r ne s’éloigne pas trop de ro, c’est le terme harmonique en (r — ro)? qui va gérer le 
mouvement par 


OV. 
=D on tuent 620 
Ici, avec le choix particulier V = -#8 que nous avons fait 
OV: 2K 3P 3 K4 
CS S OUR (3.106) 
ðr? fro rå mô IS 


on obtient pour Vr(r) 
Vie MR (3.107) 


et donc en terme des paramètres de V(r) 


(r — ro)? mè K4 


Vlr) * — r e (3.108) 
L équation de mouvement en r est mr = — aln), Nous calculons 
OV ate T m3 K4 
Tant) ~ (r— ro) J6 (3.109) 
donc l’équation de mouvement donne 
3K4 
mř x —(r — ro) Z (3.110) 
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Définissant u(t) = r(t) — ro, nous obtenons (F = ü) (voir figure 3.9) 


3K4 
mü = -u (8.111) 
16 
ou ia 
= -77 " = -Qu (3.112) 
l’équation harmonique qui a comme solution 
u(t) = r(t) — ro = Asin(Qt + ô) (3.113) 
où r(t) apparaît comme une constante, ro plus une fluctuation d’amplitude A et 
r(t) = ro + Asin(Qt + ô) (3.114) 


avec Q = mK? = fréquence angulaire de la fluctuation u(t). 


Figure 3.9 


Cette fréquence de fluctuation de r autour de r, se fait à la même fréquence que la 
rotation sur l'orbite circulaire puisque u(t) mesure la variation ou fluctuation de r(t) 
autour de ro. Ceci est caractéristique du choix particulier V (r) = -E que nous avons 
fait. On dit d’un tel potentiel qu’il génère des orbites stables. 

Cette fluctuation de r(t) autour de ro (orbite circulaire) cause un étirement de l’orbite 
à des extrémités opposées et un écrasement aux extrémités perpendiculaire (voir figure 
3.10). De plus comme Q = w, le mouvement de fluctuation est synchronisé avec la 
fluctuation et la particule revient au même point. C’est la première déformation du cercle 
vers l’ellipse que l’on sait être la trajectoire normale des planètes. On note ici que la 
proto-ellipse est encore centrée sur l’origine. Ce sont les termes asymétrique de Veg(r) , 


dont le premier est 
(r — ro)? Ver) 


qui seront responsables de ce déplacement. 


Tro 


3.5 Constantes du mouvement 


Nous avons vu quelques exemples de situation où le Lagrangien dépend d’une cer- 
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Figure 3.10 
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taine variable q mais ne dépend pas de q. On appelle q une variable cyclique et de l’é- 
quation d’Euler-Lagrange pour q 


d {OL OL 
—|[—]-——-0 3.116 
di ( a) ðq ns 
on tire, du fait de l’indépendance de L en q, que 
OL 
— =0 3.117 
Da ( ) 
et donc iat 
_{—|—0 3.118 
dt a di 
d’où nous concluons que 
ðL 
—— = constante. (3.119) 
å 


Cette constante s’appelle constante du mouvement. 
Pour un système à n degrés de liberté {q| i = 1,2....n} nous aurons n équations 


d’Euler-Lagrange 
d (ƏL ðL 
=S | — — — = 1,2...n. .12 
(5e) De = tehan (3.120) 


En mécanique, ces équations sont des équations différentielles du 2ième ordre, i.e. chaque 
équation est du type 

Pour fixer de façon unique la solution d’une équation du 2ième ordre nous avons besoin 
de deux conditions, qu’elles soient initiales, finales, limites... Techniquement cela signi- 
fie que l’intégration de chacune de ces équations requiert deux constantes d’intégration. 
Comme il y a n équations cela fait 2n constantes qui seront indépendantes puisque fixées 
arbitrairement dans le laboratoire. 


qi = qi(t, Cj, C5); 1,9 =1,2,..n (3.122) 
Ajoutant 

ġi = dilt, Cj, Ci) i,j=1,2,..n (3.123) 
nous avons 2n équations qui dépendent des 2n constantes (les n C} et les n C;). Un tel 
système peut en principe s’inverser pour obtenir 

Ci = Ci(t, q5,d;) 

Ci = Cie, dj, dj) 

Physiquement, ce sont les données d’un problème qui fixent ces constantes. Le but 

de l’exercice est d’arriver à exprimer les q; en fonction de ces constantes et du temps. 


} = 2n constantes, (3.124) 


Exemple 3.1 
Prenons l'exemple simple de l’oscillateur harmonique à une dimension 
L= TE — mur? (3.125) 


Ici, n = 1, nous n’aurons qu'une seule équation donc deux constantes d'intégration. L’'équation 
d’Euler-Lagrange est 
ï= wr (3.126) 


dont la solution peut s’écrire de plusieurs façons 


x(t) = Asin(wt+6); A,0 = const. d'intégration 
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z(t) = Bocos(wt+A); B,A = const. d'intégration 
c(t) = C'snwt+Dcosut; C,D = const. d'intégration 


Ici A, B,C et D sont des amplitudes et 8 et À des phases. Ces trois solutions sont absolument 
équivalentes. Pour les fins d'illustration prenons la dernière forme 


x(t) = C'sinwt + D coswt (3.127) 
tlt) = wCcoswt —-wDsinuwt. (3.128) 


Ces deux équations s’inversent assez facilement en 


il) 


C =  x(t)sinuwt + =o cos wt (3.129) 
D = x(t)cosut — #0 sinwt. (3.130) 


a) Conditions initiales 
Soit que, dans le problème étudié on sache qu’à un temps initial t = to la position est x(to) = to 
et la vitesse initiales &(to) = čo où to et to sont connus. On identifie facilement 


C = to sinwto + Z coswto (3.131) 
w 
to . 
D = zgzocoswto — — sin wto (3.132) 
w 
et la solution est 
3 to | 
c(t) = feo sin wto + — cos wto sin wt 
w 
to . 
+ feo cos wto — — sin usto] cos wt. (3.133) 
w 


b) Conditions limites 
Soit que dans le problème étudié on sait qu’à un temps t = to, la position est x(to) = £o et qu’à 
un autre temps t = t1, la position est x(tı) = %1 OÙ £o et xı sont connus (mesurés): 
Àt = to: zo=Csinwto + D coswto (3.134 
Àt = to: zı = Csinwts + D'cosut: (3.135) 


On peut inverser ces deux équations 
To cos wti — xı coswt 
C = RE (3.136) 
sin w(to — t1) 
zo sin wti — xı sin wt 
DS (3.137) 
sinw(to —t1) 
et la solution s'écrit 
To cos wti — xı COS Wt, 
ft) = | o 1 1 0 
zo Sinwti — zı sin wto 
RS — 
sinw(to — t1) 


| sin wt (3.138) 


sin w(to — t1) 
| coswt. (3.139) 


c) Conditions mixtes 
Toutes sortes de quantités peuvent être déterminées (expérimentalement) pour fixer la solution : 
position, vitesse, angle, énergie, moment cinétique (plus d’une dimension), etc... 
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4 LE FORMALISME CANONIQUE 


Le formalisme canonique n’introduit pas une nouvelle physique mais nous propose 
une nouvelle gamme d’outils pour étudier les phénomènes physiques. Son élément cent- 
ral, le Hamiltonien, joue un grand rôle en mécanique quantique. Comme dans le forma- 
lisme de Lagrange nous travaillerons avec des quantités comme l’énergie, T et V, plutôt 
qu'avec des quantités vectorielles comme la force F de Newton. Ici encore le formalisme 
sera invariant de forme. 


Dans le formalisme de Lagrange, la description d’un système mécanique à n degrés 
de liberté décrits par les coordonnées généralisés q;| i = 1,2, ...n indépendantes (non 
contraintes) nous mène à n équations d’Euler-Lagrange qui sont des équations différen- 
tielles du 2% ordre. 

Dans le formalisme canonique, ou de Hamilton, un système mécanique à n degrés de 
liberté toujours décrits par des q; indépendants nous mènera à 2n équations du premier 
ordre. 

Chez Lagrange on compare des trajectoires et par conséquent les q; et les g; sont 
tous indépendants (tant que nous n’avons pas résolu les équations d’Euler-Lagrange qui 
choisissent la trajectoire extremum). Chez Hamilton nous devrons d’abord apprendre à 
définir les moments généralisés, les p;, pour remplacer les ġ;, et qui eux aussi resteront 
indépendants entre eux et indépendants des q;. 


4.1 La transformation de Legendre 


Cette transformation est souvent utilisée en thermodynamique où elle permet de re- 
lier entre eux les différents potentiels thermodynamiques. En mécanique elle permet de 
définir le Hamiltonien à partir du Lagrangien. Nous en donnons une description simpli- 
fiée. 

Soit une fonction f (u, v) où u et v sont les deux variables indépendantes dont dépend 
f. Définissons 

ðf (u, 
w = raa) De w(u, v). (4.1) 
v 
La transformation de Legendre permet de définir une fonction g(u, w) qui peut remplacer 
f(u,v): 
g(u,v) =v:w— f. (4.2) 
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On vérifie facilement la chose. En effet 
Ô 
a = SE f du s pe f 


+ Su —du + wdv. (4.3) 


De la définition de g nous calculons 
dg =  wdv + vdw — df 


= wdv + vdw — de, — wdv 
Ou 


=  vdw — Î gu = g=glu,w) (4.4) 


ce qui confirme que g est bien fonction de u et de w. Pour opérationaliser cette transfor- 
mation et la disparition de v dans g on doit, à partir de la définition de w 


= 2700) = w(u,v) (4.5) 
pouvoir l’inverser en v = v(u, w) 

glu, w) = wo(u, w) — f(u, v(u, w)). 46) 
Puisque g est fonction de u et w 


_ Ôg 0g 
dg = au T pu Où =—>dw (4.7) 


et on identifie avec l’expression pour dg plus haut 
0g  Og _ of 


4.2 Le Hamiltonien 


Posons un Lagrangien L(q:, di, t) que nous traiterons comme la fonction f ci-dessus 
avec les q; jouant le rôle de u et les q; le rôle de v. A la place de w, nous définissons les 
moments généralisés 


OL i . 
Pi = dd: mn Pildi, ġj, t); (O — 1,2, en (4.9) 


un système de n équations que, comme pour v et w, nous devons pouvoir inverser pour 
obtenir les n relations 


Nous définissons donc, en analogie avec g, une fonction des q; et des p; que nous noterons 
H (qi Pis t) 


H(qi, pit -Yn — L(qi, dit) (4.11) 


dans laquelle expression g; est présumé être di (qi, Pi, t). 
De 


dL = 2 ED DE di + — dt 
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7 >. OL 
- $3 da + > pidi; + (4.12) 


on calcule dH à partir de sa définition 
nm r n : n ƏL nm | 
dH = > didpi + . - >» Ba it > nd -Ë 


> didpi — 2- a -Ea (4.13) 


ce qui vérifie que H est fonction ds di, des p; : de t). On peut donc écrire 


dH = >, da + DS di + at (4.14) 


et comme les q; et p; sont es on identifie, en comparant nos deux expressions 
pour dH 


oH ðL souati 
= ——; néquations 
qi qi 
oH 
T = i; n équations (4.15) 
OH OL 
ðt at 
On sait que la trajectoire physique obéit à l’ équation ď’ Euler-Lagrange 
ƏL d f0L d 
= — = —pi = Pi 4.16 
ðq d (ž)- dt À Sy 
et ainsi les 2n équations ci-dessus se liront 
` ƏH : 
di = p i=1,2,..n 
Di = ne i=1,2,..n a Fe 


ce sont nos 2n équations canoniques du mouvement. Ce sont des équations différen- 
tielles du premier ordre et on voit que les q; et les p; y sont traités de façon beaucoup 
plus symétrique que ne l'étaient les q; et les ġ; dans l’équation d’Euler-Lagrange. L ap- 
parition du signe moins (—) entre les équations pour les q; et celles pour leurs moments 
conjugués, s’appelle une symétrie symplectique. Les 2n équations canoniques rempla- 
cent les n équations d’Euler-Lagrange. 


De 
dH = >. da + DS di + at (4.18) 
ôðqi 


on peut, sur une trajectoire qui obéit aux E canoniques, calculer 


dH OH. oH. OH 
D = Dot) 
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ðH OL 
ôt ðt 

Ainsi, H est une constante du mouvement à moins de dépendre explicitement du 
temps, i.e. à moins qu’un agent extérieur n’agisse sur le système étudié et ce de façon 
non constante dans le temps. 


(4.19) 


4.3 Quelques exemples 


Particule soumise à une force en une dimension 


Soit une particule de masse m se déplaçant en une dimension (disons x) et soumise à une 


& force F = — y. Nous savons que son Lagrangien est 
L= Ti —V(a). (4.20) 
Nous n’aurons qu’un seul moment, noté p, conjugué à x et défini par 
OL 
= —, 4.21 
ðt ste. 
équation que nous pouvons (on doit pouvoir le faire) inverser 
=. (4.22) 
m 


On note qu’ici le moment p correspond à la composante x de la définition élémentaire 
p = mw. Ce ne sera pas toujours trivialement le cas. Selon la définition de H 


H = à(p}p— L(x,à(p)) 

= Lpa LE 

= m 2 A EVA) 
p? 

= 2— 4.2 
ma (4.23) 

que l’on écrit souvent 
H=T +V (4.24) 


où T est l’ énergie cinétique exprimée en fonction des moments. Ici, H est indépendant 
du temps et est égal à l’énergie totale, une constante du mouvement. 


Particule soumise à une force en trois dimensions 


Comme les énergies, H est additif. Ainsi pour une particule de masse m se déplaçant en 
& trois dimensions sous l’influence d’une force F = —V V (r) nous obtiendront, de 


L= TE +i +2) —V(x,u,2) (4.25) 
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que 


1 
H = (pe + pi + p2) + V(x, 2). (4.26) 
Ainsi les équations canoniques donneront (nous regardons celle en x seulement) 
H _ Prz 
er (4.27) 
pz m 
oH OV 
Dm = ——=-— 4.28 
p ox ox ER) 
Trivialement, redérivant par rapport au temps (4.27) 
se (4.29) 
m 
et utilisant alors (4.28) nous obtenons 
1 OV 1 OV 
É=-—— => ü=- = (VV): = Fr (4.30) 
mör m Ôx 


qui n’est autre que l’équation de Newton. On vérifie trivialement la même chose pour y 
et z. De plus, les moments pz, Py, pz soit ici les trois composantes de mv. 


Particule dans un champ central 


© 


La forme des équations canoniques ne dépend pas du choix qui a été fait des coordonnées 
généralisées, les q;, choix qui influencera la signification et même les dimensions des p;. 
Rappelons le cas étudié plus tôt d’une particule dans un champ central V(r) et dont le 
mouvement sera limité à un plan que nous choisissons être 0 = % ou le plan xOy. Le 
Lagrangien se réduit à 


L= Zo +r?ġ?) = V(r). (4.31) 
Nous avons vu que les équations ď’ Euler-Lagrange donnent 
mr?@ = l = constante, (4.32) 
puisque 4 est cyclique 
1? OV 
MP — —> — — =0 (4.33) 
| mrè ðr, 
(après remplacement de @). Nos moments généralisés seront 
OL 
Pr = = =mi (4.34) 
or 
ƏL 2 
= —=mr4 4.35 
Po F p (4.35) 
que l’on peut inverser en è = ® , © = 2 2- où on voit que pr et p, n’ont même pas les 


mêmes dimensions! 


On construit H selon la formule générale 


H=ŅX ġp:i-L (4.36) 
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qui donne ici 


_ Pr Po m (2y 1 | 
H = —p += po - — |(— — 
m? i mr 7 | m Tr mr? EVAN 
2 2 
pE Pr Po 
= mn t ama V(r) (4.37) 
Ici y est variable cyclique donc 
) = tant (4.38) 
— constante : 
Po 2 Po 


où on voit que dans ce formalisme une variable cyclique, en plus d’être (automatique- 
ment) éliminée, éliminera aussi son moment conjugué qui sera une constante. Deux des 
variables de H sont ainsi éliminées, ce qui n’est pas le cas dans le formalisme de La- 
grange. Les équations pour r et py sont 


T (4.39) 
pr m 
2 

is 047 07 (4.40) 


ðr mr ðr 
Si on veut comparer avec Euler-Lagrange, on dérive par rapport au temps (4.39) et on 
compare p, avec (4.40). L'identification est immédiate avec py = l. 

On voit que, comme dans le cas de Lagrange, le formalisme canonique est invariant 
de forme, i.e. il prend à son compte la cuisine algébrique qui entoure le choix de coor- 
données généralisées dont les propriétés géométriques et dimensionnelles peuvent être 
quelconques. 


Remarque 5 
Les exemples ci-dessus donnent tous 


H = T (qi, pi) + V (qi), (4.41) 


ie. le Hamiltonien est la somme de l'énergie cinétique plus l'énergie potentielle du 
système. Cette forme de H demeurera vrai tant et aussi longtemps que les interactions 
qui apparaissent dans le Lagrangien ne dépendent pas des vitesses comme dans le cas 
de l'interaction électromagnétique par exemple. 

Tant et aussi longtemps que H ne dépend pas explicitement du temps, c’est une constante 
du mouvement. Cependant H (qi, pi) west identifiable à l énergie physique que si les co- 
ordonnées généralisées, les {qi}, wont pas été obtenues de coordonnées inertielles par 
une transformation dépendant du temps. 


Exemple 4.1 


Le cas suivant en est un exemple (voir figure 4.1). Soit une bille contrainte de se déplacer sur 
une bouche circulaire (disons centrée à l’origine) et qui elle-même tourne autour de laxe Oz 
avec une fréquence angulaire w, entraînée par un moteur (extérieur). À priori, on peut écrire, en 
coordonnées sphériques 

T= D + 20° + r? sin? 0%?) (4.42) 
mais icir =a: rayon de la bouche, donc ôt = 0 et de plus & = wt donc & = w = constante. 
Donc 


N 


ma 
2 


T= (0° +r? sin? 0), (4.43) 
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Figure 4.1 
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© 


un seul degré de liberté, 0. En l’absence d'autre interaction L = T, i.e. 


ma? £ 


L=T (8° + r? sin? 6), (4.44) 
de | f 
po=mad = = 2, (4.45) 
ma 
on trouve à 2 
___Pg MA 2 . 2 
EF z% sin 0. (4.46) 


Ici H ne dépend pas du temps, c’est donc une constante du mouvement mais on ne peut pas l’iden- 
tifier à l'énergie physique de la particule parce que en posant 


p = wt (4.47) 


on fait l'équivalent d’une transformation de coordonnées dépendant du temps. En fait on se ret- 
rouve dans un repère non-inertiel puisqu'il tourne donc est accéléré par rapport au laboratoire 
(que nous considérons inertiel). 


Exercice 4.1 


Obtenez les équations du mouvement. 


4.4 Les crochets de Poisson 


Le crochet de Poisson {A, B},, est la façon standard de noter une certaine opéra- 
tion qui implique les quantités A(q;, p;) et B(q;, pi) ainsi que l’ensemble de variables 
canoniques (qi, Pi) 


” [040B  0A9B 
A,B = — ; 
(4 Bar D a Opi Opi qi aa 
De cette définition on déduit un certain nombre de propriétés 
{4, B} = {B, A} (4.49) 
{4,B +C} = {4A, B} +14,0} (4.50) 
{A,BC} = B{A,C} + {4, B} C (4.51) 
{4, {B,C}} + {C, {4, B}} + {B,{C,A}} =0 (4.52) 


oùcette dernière expression est l’ identité de Jacobi. 
Au delà d’une simple notation, leur calcul assez facile permet d’obtenir un certain 
nombre de résultats intéressants. D’autre part, ils sont intimement reliés aux commuta- 
teurs de la mécanique quantique. 
Considérons une fonction quelconque F(q;, pi, t). Sa dérivée totale par rapport au 
temps le long d’une trajectoire s’écrit 
dF ss E OF J OF 
= (3 


CE A 4. 
di aa E T Dp) T F3) 


i 
Si cette trajectoire est une trajectoire physique, elle obéit aux équations canoniques 
du Hamiltonien H du système 
. H , H 
aig ðp; Am qi 


(4.54) 
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et alors 


dE 5 0FOH 0FO0H| 0F 
dt + LOqi pi Op: Oqi ot 
ƏF 
= {F H}+ >. 4.55 
(PH}+ (4.55) 
En particulier, cette équation permet un calcul facile des constante du mouvement, 2e = 
0. En effet, le calcul de Æ est immédiat et le calcul de {F, H} est un simple exercice 


Ainsi donc, on calcule facilement . Par exemple, si F ne dépend pas explicitement 
du temps i.e. 2e = 0 alors F(q;, pi) est une constante du mouvement si son crochet de 
Poisson avec H est nul. Ceci permet d’identifier rapidement bon nombre de constantes 


du mouvement. 


Par exemple nous avons déjà vu que la conservation du moment angulaire 1 = r x p 
a comme conséquence que le mouvement est dans un plan. L inverse n’est pas vrai cepen- 
dant. Considérons un mouvement dans le plan rOy d’une particule obéissant au Hamil- 
tonien. F 
H = (pe + Py) + V (2,9). (4.56) 
Sous quelles conditions le moment angulaire 1 sera-t-il constant? Ici 1 n’a qu’une com- 
posante, soit l, où 


l; = XPy — YPr. (4.57) 
Or l, ne dépend pas explicitement du temps donc ol, =0 


iI = {L,H} 
ôl, 0H OL 0H l, 0H El, 0H 


= Á 4+4 -M NMM. 4. 
Ox Öp,  Öy Öp, Op: 0x Öp, Öy (4.58) 
On calcule 
olz əl, 
Dm = M = Per (4.59) 
olz dl, 
= —y, TS (4.60) 
OP A OPy 
et 
OH Px OH č y 
= Z Z =u 4.61 
pz m’ p m (ren 
OH OV ðH OV 
TR 4.62 
ðr ðr’ óy dy (10 
et on obtient av av 
D Pare a Pabu ee (4.63) 
m m ox oy 
l, sera donc une constante du mouvement ssi 
OV OV 
— = y—. 4.64 
£ dy Y (4.64) 


Mathématiquement cela n’est possible que si, indépendamment de z qui n’apparaît pas 
ici, V(x, y) ne dépend de x et de y que sous la forme (x? + y?)”. 


Si on s'intéresse aux trois dimensions on obtiendra des conditions sur ly — V ~ 
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(y? +22)" etsur ly — V ~ (x? + 2?)®. Pour avoir conservation de 1 on doit donc avoir 
V = (x? +y? + 22)" ~ r?”: un potentiel central. 

Il existe toute une famille de résultats intéressants du crochet de Poisson. Parmi les 
plus importants, calculons certains de ces crochets entre des variables canoniques : co- 
ordonnées et moments: {qk, qj} , {pr,p;} et {qr,p;} où k et j sont fixés 


7. faq ðq;  Əqk 0q; 
M ` = = 4. 
{ar 95} = E Op; Op; ðqi 9 FES 


puisque 
o o 
ti o, JE =0 (4.66) 
Opi pi 
parce que les variables canoniques sont indépendantes et 
{Pr p;} = 0 (4.67) 


pour la même raison, mais 


a S ôa 0p; 3dr ôp, 
{ar,p;} = > E Op; Op: qi 


i 


+= Ôk Op; 
~ Üqi pi 
= 5 Ôkiðji = Ôkj (4.68) 
où ôk; est le delta de Kronecker: | 
ë = { =. (4.69) 


Ces résultats sont très importants parce qu’on peut démontrer que leur inverse est vrai, 
i.e. si un ensemble de n q; et de n p; obéit aux relations ci-dessus, alors l’ensemble des q; 
et des p; constitue un ensemble de variables canoniques. Ceci est très important et trouve 
éventuellement des applications dans les théories quantiques du champ. 

Une autre utilisation intéressante des crochets de Poisson est qu’ils permettent de 
symétriser les équations canoniques. Puisqu’il est vrai que pour une fonction quelconque 
des variables (p;, qi), soit F(p;, qi), sa dérivée par rapport au temps est donnée par 


F={F H} (4.70) 


la chose est certainement vrai pour les q; et les p; eux-mêmes et les équations canoniques 
peuvent s’écrire 


di = {&,H} (4.71) 
D = {p H} (4.72) 


A cause de cette symétrie on est parfois amené à parler des 2n variables canoniques. 
Une telle symétrie n’existe pas dans le formalisme Lagrangien entre les q; et les di. 
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4.5 Les moments généralisés 


On parle des p; comme étant des moments généralisés de la même façon que les q; 
sont des coordonnées généralisées. Comme les dimensions des q; peuvent être à peu près 
n’importe quoi, il en va de même des p;. 

Dans les exemples que nous avons vu les p; étaient les composantes de p 


p = mv (4.73) 


ce qui est particulièrement évident en coordonnées cartésiennes. Mais même en coor- 
données cartésiennes cette définition n’est pas toujours vraie. En effet, lorsque l’interac- 
tion dépend des vitesses p Æ mv. L'exemple le plus important est sans doute celui des 
interaction de jauge. En effet, nous avons vu que le Lagrangien d’une particule de masse 
m et de charge e dans un champ électromagnétique est, en coordonnées cartésiennes, 
Ti = (x, y, 2) 

L= oi te) iAV (4.74) 


i i 
où V et A sont les potentiels scalaire et vectoriel du champ électromagnétique et dépen- 
dant généralement des x; et de t. 
Définissant les moments généralisés p; 


ðL 
Pi = — = Mi; + eÅi (4.75) 
Oti 
on constate que p n’est plus mv mais 
p= mv + eA. (4.76) 
Ces équations s’inversent en 
i — eÅ; 
p= = (4.77) 
Avec la définition de H = $; pit; — L, on obtient 
1 2 
H = D (Pi -eA) + eV. (4.78) 


K3 


Exercice 4.2 


Vérifiez ce résultat et vérifiez que les équations canoniques redonnent les équations du mouvement 
d’une particule soumise à une force de Lorentz. 


Ce résultat est important puisqu'il nous dit comment écrire le Hamiltonien pour une 
particule soumise à une interaction de jauge. Aujourd’hui on cherche à écrire toutes les 
interactions qui apparaissent dans la nature comme des interactions de jauge. 


4.6 Les transformations canoniques (T.C.) 


On dit que des q; et des p; que ce sont des variables canoniques généralisées. Ce n’est 
pas un euphémisme puisqu'il n’y a pratiquement aucune limite à ce qu’elles peuvent 
représenter physiquement. Nous en venons d’ailleurs quelques exemples. Puisque tel est 
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le cas il doit exister des transformations entre ces différents choix. Nous noterons Q; et 
P; les nouvelles variables canoniques obtenues suite à une telle transformation. 


On n’est pas surprit par contre de constater que ces transformations sont soumises à 
des conditions assez sévères. En effet les q; et p; sont généralisés et obéissent à 


{ax qj} =0, {Pk Pj} =0, {qp;} = di (4.79) 
et les équations canoniques 


OH OH 
= ES 4.80 
Dp, Di Ög; (4.80) 
sont invariantes de forme. Ainsi, à la suite d’une transformation des q; et p; vers les Q; 
et P; et définissant un nouvel Hamiltonien que nous noterons K (Q;, P;) nous devrons 


avoir 


di 


{Qk Qijap =0, {Pk Pijan =0, {Qr Pi}ap = di (4.81) 
et les équations canoniques 
oK . oK 


aap" Ro (4.82) 


Qi 
Strictement les équations de transformation peuvent s’ écrire 


Qi = Qildi, pj, t) AE 

i,j=1,2,..,n 4.83 
P; = P;(q;, pj, t) >] ETES ( ) 
et doivent pouvoir s’inverser puisque la physique reste indépendante des variables qu’on 
emploie pour la décrire, donc on doit pouvoir écrire les transformations inverses 


q=œ(Q;, Pit) .. 
à, =1,2,..,n. 4.84 
pi = pi(Q;, P;,t) i te 


Les q;, Pi, Q; et P; forment 4n variables mais il est évident que seules 2n d’entre elles 
sont indépendantes. D’ autre part s’il est généralement possible d’écrire par exemples les 
n équations de transformation 


Qi = Qi (dj Pj, t) (4.85) 
de façon assez arbitraire il est généralement impossible d’écrire les n autres équations 
P; = P;(qj, pj, t) (4.86) 


de façon aussi arbitraire puisqu’un tel choix ne satisfera pas en général les conditions 
énoncées plus haut. 


Il faut donc apprendre à faire correctement ces transformations. La façon standard de 
le faire est de considérer que pour les fins de la transformation n des anciennes variables 
et n des nouvelles sont linéairement indépendantes, par exemple les q; et les P; alors que 
n anciennes et n nouvelles restantes sont linéairement dépendantes, ici les p; et les Q;. 
Dans ce cas précis nous écririons donc les équations de transformation 


pi = pilq, Pj t) (4.87) 
Qi = Qig;, P;t) (4.88) 

Pour obtenir la forme habituelle on inverse les n premières en 
Pj = Pj (qi, pi,t) (4.89) 
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que l’on remplace dans les n dernières 


Qi = Qi(a5, P (de, Pk, t),t) = Qildi, pj, t). (4.90) 


Pour générer ces transformations, nous retournerons au principe variationnel lui- 
même. Nous savons que 


2 
ôS = af Ldt = 0. (4.91) 
Par ailleurs de 2 i 
H= pġ- L où pi= a (4.92) 
7 ôdi 
on peut obtenir 
L= n -H où Gi = 2L (4.93) 
7 pi 


ce qui permet d’écrire 


2 n 2 
ôS —6 | D Didi — J dt =6 | Ldt = 0. (4.94) 
1 3 1 


Or si L correspond à H, L’ correspondra à K et nous exigeons d’avoir également 
2 n 
ô | L'dt=0 où V= 5 PQ; -K. (4.95) 
1 i 


Pour que L’ et L décrivent la même physique nous avons déjà vue que L et L’ ne peuvent 
différer l’un de l’autre que par la dérivée totale d’une fonction F, i.e. 


dF 
L=L'+ ET (4.96) 


Nous poserons donc 


2 n 2 
‘| nanas f 
1 7 1 


Cette fonction F, que l’on appelle le générateur de la T.C. sera choisie comme ne dépen- 
dant que des variables indépendantes de la transformation. On identifie généralement 


Variables Variables Générateurs 
ea permes 


dt (4.97) 


Zira dF 
Dino 


quatre cas 


Fi (qi, Qù,t) 
Fola, Pit) 
Fs(pi, Qi, t) 
Fa(pi, Pit) 


Étudions un peu plus en détails les cas F} (q, Qi, t), Fo (qi, t). Dans le cas F; (qi, Qi, t) 


DE or Igt D ar 10 . (4.98) 
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Ici, il est suffisant de comparer les fonctions à intégrer dans ô f Ldt = 0 de qui donne 


> pidi — H(t pit) = $ Pà- K(Qi, Pi, t) D 


2 Q+ a IL (4.99) 


Cette équation est satisfaite si on identifie, q;, Q; et donc ġ; et Q; comme étant indépen- 
dantes. Les facteurs de ces variables indépendantes, doivent donc être identiques 


o 
Pi = gg (0 Rit) = Pildi Qt) 
Ô 
P; = | 2Q; Fı (qj, Qj, t)= P;(g;,Q;, t) (4.100) 
E OF, 


Clairement ces lois de transformation nous permettent d’écrire les 2n variables dépen- 
dantes (ici les p; et P;) en fonction des 2n variables indépendantes (ici les q; et Q;). Ces 
2n équations peuvent se mettre sous la forme plus habituelle 


di = q(Q;,P;,t): n équations (4.101) 
pi = pi(Qj, Pj,t): n équations (4.102) 
ce qui permet de calculer K 
K(Q:,P;,t) ma H(qi(Q;,P;,t),p:(Q;,P;,t),t) 
o 


En étudiant les équations de transformation obtenues ci-dessus, on constate que 
Op: ə? F 1 OP. j ©? F 1 


= $ = — (4.104) 
oQ; qQ; di 0Q;0q: 
ainsi donc, un test du caractère canonique de la transformation 
pi = pi(g;, Pt): n équations (4.105) 
P; = P(g;,P;,t): n équations (4.106) 
est qu’elle doit satisfaire 
Op; oP; 
es (4.107) 
0Q; qi 
ce qui est équivalent à 
{Qi Qj} =0 = {P;, Pj} et {Qi, Pi} = ôiz. (4.108) 


Dans le cas F3 (qi, P;, t) ce sont les q; et les P; qui sont considérés indépendants. On 


calcule 


SE D ed DT + a (4.109) 


Ici la comparaison des RS à oA n "est pas suffisante et nous devons récrire au 
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complet 


sf DD na d = 
1 i 


Ar E he (4.110) 


Le problème vient de ce que les Q; ne sont pas considérés indépendants ici et donc les 
Qi ne le sont pas. Intégrons par partie le premier terme à droite 


2 Nn n 2 2n 
‘| > P.Q;dt = D PQ; | ÿ P.Qidt. (4.111) 
1; i 1 1 ił 


=0 pcq points fixes 


Maintenant nous pouvons comparer les fonctions à intégrer 


i OF. 
Di —A = $ Ra- K +) t 


+ 2o a +72 (4.112) 
et identifier les facteurs des variables indépendantes, d; et È; 
pi = Pala Pyat) = pilan Ppt) 
Q = -3b Palan Ppt) = Pilay Ppt) (4.113) 
OF: 
K = H+ a 


Ici encore les variables dépendantes apparaissent exprimées en fonction des variables 
indépendantes. Comparant les expressions pour p; et Q;, le test du caractère d’une trans- 
formation de type F2 est 


Opi — Qj 
a a (4.114) 
{1Q:,Q;}=0 = {P,P;} et {Q:,P;} = di. (4.115) 


Remarque 6 

Les fonctions Fı (qi, Qi, t), Fo (qi, Pi, t) etc, ne génèrent pas des transformations diffé- 
rentes mais sont simplement des façons différentes de générer une transformation donnée. 
Évidemment, deux fonctions F, différentes par exemple vont en général générer des 
transformations différentes. 
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Quelques exemples 


a) Fréquemment on cherche à effectuer une transformation de coordonnées, i.e. on 
connaît les fonctions 
toujours inversibles en 


qi= qŅ(Q;; t), i,j =1,2,.. n. (4.117) 
En général il n’est alors pas possible de poser à priori les équations de transformation des 
moments, on doit s’assurer que ces derniers seront des moments canoniques généralisés. 


Une façon ”simple” de procéder est par le biais d’une transformation de type F2(q;, P;, t) 
définie comme 


Fo(qi, Pat) = X Qi(qt)Pi. (4.118) 
Dans ce cas les équations canoniques de transformation seront 
OF. 
Qi = = Qi(qj,t): tel que désiré (4.119) 
ðP; 
et 
0F> +=, 0 
= = 2 Fr 25 (art) = pi(P, at). (4.120) 
Ces n dernières équations peuvent s’inverser en 
Pi = P;(q;,p5;t) (4.121) 


complétant ainsi l opération et garantissant que les P; ainsi définis seront canoniques. 


Exemple 4.2 

Voyons un exemple simple (voir figure 4.2), celui du passage aux coordonnées cartésiennes en deux 
dimensions, x et y aux coordonnées polaires, r, p. Ici qi = (x, y) et Qi = (r, p) avec i = 1,2. 
De plus, pi = (Pa, Py) et P; = (P,,P,). Nous savons que 


Figure 4.2 


qı = T = r cos Y 


EE } du type qi = qi(Q;) et { (4.122) 


P2 = Py 


® ® 
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donc 
r = (+ Qi = (+) (4.123) 
p = tan™ (2) =S Qta (B), (4.124) 
T qı 
Nous écrivons la fonction Fz(qi, P;) = F2(x, y, Pr, Po) 
Pal) = (È+ ÈP. tta (B) p, 
1 
= F(x,y, Pr, Ps) = (x? + ÿ)?P, + tan”! (2) Pa 
Les lois canoniques d’une transformation F2 sont 
__ 0P 
P£ = PA 
_ 7 m y 
(£2 + y2)2 (z? +y?) ° 
= cospP, - EP, (4.125) 
_ ôR 
Py = y 
2 y g 
(x? + y2)? CES 
= snpP.+ Er, (4.126) 
d’où on obtient facilement 
P, = er (4.127) 
(x? +y?) 2 
P; =  XPy — YPx = (r Xp): (4.128) 
où (r x p), = composante z du moment angulaire. Supposons de plus que nous ayons 
l,a 2 
T = — 4.12 
+p?) (4.129) 
alors on calcule facilement 
T= (P? + À) (4.130) 
o m" r j 


que l’on sait déjà être le bon résultat. 


Exercice 4.3 


Calculez les {Q:, Pj} pour vérifier que les nouvelles variables sont canoniques. 


Exemple 4.3 


Quelques exemples illustrateurs sur l’oscillateur harmonique dont le Hamiltonien (1 dimension) 
est 


H= SE + D. (4.131) 
Tentons de passer des variables canoniques x et py à de nouvelles, notées q et p par 
na JE es emo (4.132) 
mw 
1 
P: = VMWP = p= (4.133) 


Jma 


Copyright © 1997 P Amiot, L. Marleau 


74 


Chapitre 4 LE FORMALISME CANONIQUE 


On vérifie que 


ôq p ôq p 


torn} = > D an (4.134) 
1 


de qui est correct et de toute évidence {q, q} = {p, p} = 0, donc la transformation est canonique 
et H devient K : 


K=H+ TH (4.136) 
K(q,p) = Sp? +g’). (4.137) 
La solution est triviale IK JK 
ÿ=— = wpetp = wq (4.138) 
Op ôq 
alors 
à = wp = -w°q (4.139) 
d’où 
q(t) = Asin(wt + ô) (4.140) 
donc 
x(t) = VmwA sin(wt + 6). (4.141) 
Exemple 4.4 
Au lieu de cette transformation, essayons plutôt de passer de (£, Px) à (q, p) définis par 
z= 2 et Pr = V2mp (4.142) 
mw 
ou 
mw? Pz 
= t p= 7 4.143 
q y TARR ( ) 
Nous passons alors de H à K défini par (Æ = 0 => K = H) 
K(q,p) =  H(x(q),p(p)) 
2 
=. + 2 mM 2 2 
o t 2 mul 
= p+? (4.144) 
La solution est triviale ok ok 
ÿ=———=2p, p — 2q (4.145) 
p ôq 
donc 
ÿ = 2p = —4q (4.146) 
donc 
q(t) = Asin(wt + ô). (4.147) 


On s'attend à ce que q(t) = N sin(wt + ô) où w # 2 en général, donc ce résultat est faux. La 
raison est que la transformation faite ici, même si elle semble très simple, n’est pas canonique. En 
effet on vérifie que 


ôq Op ôq Op 


(PP Dion on 0 
muw? 1 w 
= == 002i 4.148 
= = 17 ( ) 


sauf pour le cas particulier w = 2. 
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Reprenons l’Hamiltonien K (q, p) correctement obtenu précédemment et supposons que w = 1 


1 
= Klap) = 3 +) (4.149) 
dont la solution sera trivialement 
q = Asin(t +6). (4.150) 
Faisons une T.C. additionnelle définie par 
1 ; 
Q= 0t ip) (4.151) 
et i 
i 
P = —(q — i 4.152 
34 ip) (4.152) 
qui s’inverse facilement en 
Q—iP i(Q +iP) 
= ; = — << 4.153 
e (4.153) 
Même si la transformation est complexe on vérifie facilement que 
_ QƏP aQƏP 
(P) = Go dp dq 
Re © 
V2 V2 V2 V2 
1 1 
= 21. = 1; (4.154) 
Trivialement le nouveau K, noté ici Kı est 
K1(Q,P) = —iQP (4.155) 
Les équations du mouvement seront 
> __ OK1 _ . 5 _ _ OK _  . 
Q = = iQ, P= 5Q = +iP (4.156) 
donc | 
Q(t) = Ae”, P(t) = Be”. (4.157) 
Ainsi m m 
Ae “ —iBe 
t) = ——————. 4.158 
q(t) 7 (4.158) 
Par les définitions même de Q et P, P = iQ*, donc 
A*e* = Be” > B=iA* (4.159) 
et . i 
Ae™™ + A*e' 
t) = —— =. 4.160 
q(t) 7 (4.160) 
Mais À est une constante (complexe) que l’on peut écrire 
A=|AjeS = A*=|Aje (4.161) 
ce qui donne 
q(t) = |A| cos(t — A) (4.162) 


qui est une bonne solution. Cette façon de faire peut sembler étrange mais en plus de mélanger 
coordonnées et moments, elle trouve une application en mécanique quantique. 


4.7 Une transformation canonique très spéciale: La 
méthode de Hamilton-Jacobi 
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L'objectif 


Les transformations canoniques ont pour but de simplifier les problèmes. L'une d’entre 
elles est tellement systématique qu’elle porte un nom, la méthode de Hamilton-Jacobi. 

Soit un Hamiltonien (H (qi, pi) dépendant de 2n variables canoniques, les (qi, pi), 
i = 1,2...n. Nous savons que nous pouvons passer à un nouvel ensemble de variables 
canoniques, les (Q;, P;) également au nombre de 2n et dont sera fonction un nouvel 
Hamiltonien K(Q;, P;). Par ailleurs nous savons que le système compte 2n constantes 
du mouvement. Le but de la méthode est d’opérer une T.C. telle que les (Q;, P;) soient 
précisément 2n constantes du mouvement. Si tel est le cas, alors 


Qi E 0 Qi = b; = constantes (4.163) 
È, = -2 0 > P; = Qi constantes (4.164) 
d’où 
di = qi(Q;, Pj j t) (4.165) 
2 arand (4.166) 
pi = pi(Q;,P;,t) (4.167) 
= pilbj aj t) (4.168) 


ce qui trivialise au maximum les équations du mouvement dans la cadre Q;, P; et K (Q;, P;). 
Pour y arriver nous chercherons la fonction génératrice, ici choisie de type F2(q;, P;, t) 
donc du type F2(qi, œi, t), que nous noterons de façon standard S(q;, &;, t), telle que 


K(Qi Pnt) = Hlp t) + D = 0. (4.169) 


La méthode 


La T.C. est de type F2 et donc 
OS 


ðq 
Le but recherché est H + 25 = K = 0, ce sera notre équation fondamentale après 
remplacement des p; dans H 


Pi = (4.170) 


ðs aS (qi, ai, t) 
A e A 
Gr x 


c’est l’équation de Hamilton-Jacobi, une équation différentielle pour S. Une fois solu- 


=0 (4.171) 
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tionnée i.e. une fois que l’on connaît S il ne reste qu’à opérer les T.C. 


OS(q;,a;,t 
Pi — ʻi 4 ) = pi(qj, aj, t) (4.172) 
qi 
et 
OS  OS(qj, aj, t 
Q, = p= 2S 2Sa t) 
OP; Oai 
= Qil, aj, t) = bi (4.173) 
Ces n équations peuvent s’inverser en 
qi = qi (aj, Bj, t) (4.174) 
ce qui est la solution! Si on veut les p;, on remplace dans les résultats de la T.C. pour 
Pi = pilti, ai, t) = pi(g;(ar, Br t), Qi, t) 
= pilan Bit). (4.175) 


Une simplification importante apparaît lorsque où = 0. Dans ce cas, par séparation de 
variables, on peut écrire (puisqu’alors H est une constante donc de H + 2 = 0 => 


Sa~t) 

S(qi, ai, t) = W (qi, ai) — aıt (4.176) 
où on identifie l’un des a, soit ici «1 comme la valeur numérique (constante) de H : 
H= Q1. 

Comme as aw 
i p= 4.177 
Di in Me ( ) 
T équation de Hamilton-Jacobi devient simplement 
ow 
H (qi, —) — ai = 0. 4.178 
(q Üa )— a ( ) 


C’est l’équation caractéristique de Hamilton-Jacobi pour la fonction W (q;, a;). La sim- 
plification peut aller plus loin. En effet, toujours par séparation de variables on constate 
que si une coordonnée, disons qk pour k fixé, est cyclique, elle n’apparaît pas dans H et 
pk est alors une constante qui peut être utilisé comme px. Dans ce cas on peut écrire la 
dépendance en W sur qk simplement comme 


OS oW 
W ~ azqr Pk = =— = =— = Q : constante (4.179) 
0qk  Əqk 
et de façon générale W s’écrira 
cycliques 
W(qirai)= XO arge + W(q;,a;) (4.180) 
k 
où les qp cycliques n’apparaissent pas dans W”. 
Exemple 4.5 
&& Illustrons la méthode par un exemple simple, soit un problème physique décrit par 
ad 2 2 2 mw? à 2 
H = (ps + py + pz) + (£? +y’). (4.181) 


2m 2 
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Ici, H ne dépend pas du temps, donc on peut écrire 


Si, ai, t) = S(x,y, Z, ai, t) = —ait + W (zx, y, z, ai). (4.182) 
De plus, z est variable cyclique et nous pouvons écrire 
W (x,y, z, oi) = azz + W'(x, y, as). (4.183) 
Nous aurons donc Si w 
Pe = -py PS py P5 (4.184) 
L'équation caractéristique de H.-J. sera, de 
oW 
H(qi, —) — a1 = 4.1 
(q 2 )— 1 =0 (4.185) 
1 fOW'\?, 1 fOW'\? oœ ma? a o 
=) +) tant az e +#)-a=0 (4.186) 


C’est une équation différentielle (non linéaire). Cependant la forme relativement simple de l’équa- 
tion permet d'espérer qu'une séparation de variables 


i ðw’ _ 3Wz(2) 
W = Wz (x) + W, (y) = aw E ou) (4.187) 
ðy — ðy 
donnera des résultats. On obtient en effet alors, regroupant 
1 oWz LE 1 0Wy ? _ mu? 2 4 
2m \ 0x 2 2m \ y z” 

CE 

+ -a = 0. (4.188) 
2m 


Les deux premiers termes contiennent toute et seulement la dépendance en x, leur somme doit donc 
être égale à une constante que nous appellerons a3. Ceci nous laisse 

2 2 
1 /ƏW;, mw 
— Z) + r = 0 (4.189) 
2m \ Ox 2 


2 2 2 
— (z=) + y + 08 + E -o =0. (4.190) 
Ainsi donc les deux termes en y sont aussi égaux à une constante mais ici il n’est pas nécessaire d’en 
introduire une nouvelle. Nous aurons donc trois constantes œi = (&1, &2, &3) qui représentent les 
trois nouveaux moments, P;, ce qui est correct puisque nous avons trois degrés de liberté. Isolant 


les dérivés ci-dessus nous obtenons 


sie nor (4.191) 


ma = 2maı — 2ma — a? — m?w?y? (4.192) 

ou 
W: = J mog — M?w?r? dr (4.193) 
W, = | 2maı — 2ma? — a? — m?w?y?dy. (4.194) 


Souvent il west pas nécessaire de faire ces intégrales puisque nous wavons pas besoin de W en 
soi. Ici, S s'écrira donc 


S = -ort toaz + | Vamai- murs 
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+ f 4/ 2Mmaı — 2ma? — a? — m?w?y? dy (4.195) 


Souvent, il reste à appliquer les règles de transformation pour une T.C. de type F2 i.e. où nous 
savons que 


Qi = 7 = B, : constante (4.196) 
ce que se lira ici, avec @i et B; constantes í 
pae. (4.197) 
ai 


Explicitement nous aurons donc 


6, = ðS A Í dx 
dB : V/2ma2 — mwg? 
—2maz T —__  Ọ_ 
2maı — Zma? — a? — m?w?y? 
(4.198) 
d 
h= mra | a (4.199) 
Oaz 2maı — 2ma? — a? — m?w?y? 
d 
B= tm |- (4.200) 
ðaı 2ma1 — 2mai — a2 — m?w?y? 
Si on intègre les équations pour B; et B>, nous obtiendrons deux expressions du type 
f(æ,y,a1,02,a3,83) = 0 (4.201) 
g(y,2,@1,02,a3,83) = 0. (4.202) 


Comme nous avons 3 dimensions, ces deux équations satisfaites simultanément nous laissent un 
espace à une dimension: la trajectoire, exprimée en termes de x, y et z, sans la dépendance en 
temps. En d’autres termes f = 0 et g = 0 laissent une des coordonnées indépendante, disons y et 
inversant f et g on peut en principe écrire 


x = x(y, ai, b2, ba) (4.203) 
z = Z(Ņ,@i, b2, ß3). (4.204) 


La dernière équation, celle en 6, donne y = y(t, œi, B1). 

C’est de là qu’on obtient le développement dans le temps de la trajectoire. Dans un certain nombre 
de cas cette dépendance en t west pas le but recherché et on peut alors se limiter aux deux premières 
qui nous donnent la trajectoire uniquement en fonction des coordonnées. 

Voyons voir ce que cela donne ici. 


2a3 . —1 2 z] 2a3 . _1 y 
= —— SIn — ——— sin -p |, 4.205 
PT as v2 w m ue a (En) 
mu? mw? mw 
2@2 . 1 y 
Ba = zZz — —— SIn aa aaea (4.206) 
Ww y 2a1 2a2 a 
mw? mw? mw 
et 
1 = 
Bacaan ea (4.207) 
(02) 2@œ1 2a2 a 
mw? mw? mw? 
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Des deux premières expressions nous tirons z(y) et z(y) après un peu d’algèbre élémentaire 


2 
z(y) = A sin vps +sin | 2 U (4.208) 
mm 2a3 Jai 2a3 a2 
@2 . —i y 
2) = Ba + mo y Joy  2a2 a2 € ) 
MW = MW = Mme 


qui sont les expressions donnant la trajectoire sous la forme x = z(y) et z = z(y). L'équation en 


B; donne 
241  2a ai . 
Den un us sin(wt + wb) = y(t) (4.210) 
Dans le problème étudié, il est évident que le mouvement en x et y est harmonique de fréquence et 
le mouvement en z est libre. Clairement la solution y(t) est de la bonne forme 


y(t) = Yo sin(wt + 6). (4.211) 
Remplaçant ce résultat dans les solutions z(y) et z(y), nous obtenons 
PEN A fiot + wc, + £) (4.212) 
m w 2a3 
aussi de la forme 
x(t) = Xosin(wt + A) (4.213) 
et 
a 
z = 2+ UHN 
= iy (5, m Pi) (4.214) 
m m 
de la bonne forme 
z = Viet + zo (4.215) 


où vz est une vitesse constante. On voit directement ici comment relier les constantes @; et b; aux 
conditions initiales du problème. 


L'exemple ci-dessus est simple mais il illustre de façon claire que pour la première 
fois nous obtenons la trajectoire sans passer par une équation du mouvement se ramenant 
à F = ma. La méthode est apprécié pour son intérêt théorique, sa relation avec l’ optique! 
et lorsqu’on cherche la trajectoire sous la forme 


x= z(y), z=2(y) (4.216) 


plutôt que sous la forme 


= y(t) (4.217) 


4.8 T (qi, pi) en coordonnées généralisées 


Nous avons vu dans le cadre Lagrangien que l’ énergie cinétique s’écrit de façon géné- 
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rale en fonction des vitesses q; : 
m a: 
= 7 D tij. (4.218) 
2 
Ainsi les moments généralisés a d'interaction dépendant de v) sont 


= ikli 4.219 
Pk = de = m > gni ( ) 


en utilisant la symétrie gik = gpi (espace de ou On peut écrire cette équation de 
façon matricielle 


p = mgg (4.220) 
où 
pı gı 912 t Jin di 
p2 g21 92 ‘‘* on d2 
ou RTE ; E O E a ll & (4.221) 
Pn ni In2 °° Onn On 
et m est la masse, un simple nombre. Multipliant de la gauche par s où g7! est la 
matrice inverse de g, i.e. 
g lg = gg ! = I = la matrice identité (4.222) 
on obtient i 
9 D à (4.223) 
ou de façon explicite 
1 -1 
&=— > (g he Pj- (4.224) 
j 
Utilisant la notation matricielle on peut écrire 
T= pa (4.225) 
et donc en formalisme de Hamilton où T = T (p) nous aurons 
a Ce Ka D 
Dr 
F =p T (g7 1) g gg 
1 11T 
= 4.226 
a T (gt) p (4.226) 
ou explicitement 
1 LT 
T= RD mg Ps 
ij 
1 = 
= nd Pig Ps (4.227) 
ij 


Lorsque g est diagonal, ces opérations sont encore plus simplifiées puisqu’alors g;; = 
Jiiðij et que les opérations de transposition sont sans effet. Par exemple nous avons vu 
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qu’en coordonnés sphériques, qi = (r, - A ), la rar 


. (4.228) 
r? sin? 0 


Trivialement 
1 0 
0 + (4.229) 
0 0 zim 
donc ici 
T = 1- p? + Li + ni (4.230) 
TU p20 l p2gin 0 ? i 
ce qui est le bon résultat et F la forme 
1 —1 
T = > Pig ÔijPj 
ij 
1 = 
=a >, Cr lp? (4.231) 
ij 


qui est générale pour les cas où la métrique g est diagonale. 


4.9 La fonction S (ou comment refermer la boucle) 


La méthode Hamilton-Jacobi voit apparaître une fonction génératrice de transforma- 
tion canonique et noté S. Nous avons déjà utilisé ce symbole pour désigner l’action. Ici, 
S est défini par 

OS 
H + D — = 0 (4.232) 
où S = S(qi, ai, t).Nous avons donc, calculant la dérivée total de S par rapport au 
paramètre t, 


aS | yo, , SaS 
d ` ag Oo" Dt 
OS OS 
— à; + — 4.2 
D Bat + z (4.233) 
puisque &; = 0, et 
OS OS 
= pi — = -H 4.234 
Ba: pi et St (4.234) 
dS . 
ou ’ 
dS = Ldt (4.236) 
ou , 
S= Ldt. (4.237) 


tı 
La fonction génératrice de H.-J. est donc simplement l’action. Formellement intéressant, 
ce résultat est cependant pratiquement inutile parce qu’il faut avoir complété la solution 
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du problème pour la vérifier. 
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5 THÉORIE DES PERTURBA- 
TIONS 


5.1 Buts de la méthode 


Il s’agit d’une méthode approximative pour obtenir une solution analytique à un pro- 
blème de mécanique qui n’a pas de solution analytique exacte ou pour lequel cette solu- 
tion est trop difficile à obtenir. En fait il n’existe que très peu de problèmes de mécanique 
qui ont une solution analytique exacte. Le problème à trois corps par exemple n’a pas de 
telle solution. Les ordinateurs d’aujourd’hui permettent de résoudre numériquement ces 
problèmes avec pratiquement la précision désirée mais à chaque fois pour un ensemble 
donné de conditions initiales. Pour avoir une vision générale du type de trajectoire, il faut 
faire plusieurs fois les calculs et ceci peut être onéreux. 


5.2 l’idée de base: la variation des constantes 


Soit un système décrit par un Hamiltonien H (q;, pi), i = 1,2...n. On sait que la 
solution du problème dépend de 2n constantes d’intégration, appelons-les les a; et les 
bi, i = 1,2...n, qui sont évidemment des constantes du mouvement. La solution devrait 
alors s’écrire 

qi = qi(t, aj, bj) 
Pi = pi(t, Qj, b;) 

Supposons que nous soyons incapables d’obtenir ces solutions analytiques mais que 
pour des raisons du type énumérées en ci-haut, nous désirons obtenir une solution ap- 
proximative et analytique. La méthode des perturbations peut permettre d’obtenir cette 
solution approximative. Elle requiert que nous soyons capables d’écrire 


H(qi.pi) = Ho(qi-pi) + Hı (qipi) (5.2) 


ij =1,2,. n (5.1) 


de façon telle que 


1. il soit possible d’obtenir une solution analytique pour Ho et, 


2. que H soit petit devant Ho. 
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Cette dernière condition est souvent difficile à vérifier à priori. Elle requiert une cer- 
taine stabilité du mouvement face aux changements dans les conditions initiales et de fait 
la méthode n’est pas appropriée au traitement des mouvements chaotiques par exemple 
qui sont caractérisés par une très grande sensibilité aux conditions initiales. 

L'idée de base est relativement simple et elle compte les étapes suivantes: 

i) On résout analytiquement pour Ho et on obtient les solutions 


de (0) b0 ) 


di =q 
( (5.3) 
Pi = pa ro pO), 
ii) On inverse ces 2n équations pour obtenir les 
AE a(t, qe), p®) … 
bi = 00t, q9, PO, (a 
qui vérifient évidemment (ce sont des constantes) 
— [a0 = a® 
& = {a; ", Ho} + S (5.5) 


bi = {00 Ho} + 4 HO 
iii) On pose que la solution complète pour H peut prendre la même forme que celle en 
i) et ii) mais avec des aP et p® remplacés par des a; et b; qui ne sont plus des constantes 
du mouvement i.e. pour lesquels à; Æ 0 Æ b;. 
iv) On calcule les a; et les b; par leur équation du ”mouvement impliquant H au 
complet et dans lesquelles les a; et b; sont présumés avoir la même dépendance dans les 
qi et p; qu’en ii). Ainsi 


Ô Ta 
ài = {a H} + 5 ui = à = ai, Ho} + {ai H1} + = (5.6) 
Mais selon la seconde pe 
Ô 
{ai, Ho} + ai = 0 (5.7) 
et il ne reste que 
ài = {ai, Hı} (5.8) 
et . 
bi = {b;, Hi}. (5.9) 


v) Il reste à intégrer ces équations pour obtenir a;(t) et b;(t) et à les replacer dans les 
équations (i) en lieu et place des aP et p® pour obtenir la solution désirée 


a = G, a, O) 


pi = pa al 


mêmes fonctions que pour la solution non-perturbée mais ici a; = a;(t) et b; = b;(t). 


(5.10) 


5.3 Les approximations 


À ce point-ci, il n’y a aucune approximation de faite. Elles apparaissent dans l’inté- 
gration des équations pour a; et b; souvent elles-mêmes trop difficiles pour être résolues 
exactement. 
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On présente souvent la méthode perturbative comme l’approximation d’une expan- 
sion en série de puissance d’un paramètre qui caractérise H1. C’est d’ailleurs généra- 
lement le cas en mécanique quantique. Ce n’est pas cependant la seule approximation 
possible. Mentionnons la méthode itérative et celle de la moyenne, cette dernière étant 
utile lorsque les trajectoires de Ho sont des orbites fermées. 


Méthode par série 


Sous une forme simplifiée, on peut la présenter de la façon suivante. On identifie 
d’abord un paramètre À, idéalement sans dimension (et petit) tel que 


Hi = Xh(qi, pi) (5.11) 
et on pose que l’on peut écrire les a; (et les b;) en séries de puissance 
ai = a + Aah Aa H (5.12) 
Remplaçant dans l’ équation pour å; on obtient 
O Aah Aa oe = Ma + Aa Aa Hhh. (5.13) 
Égalant les termes en même puissance en on obtient 
a0 = 0= a® = = constante 
aD = {ah} 
aP = {al,h} (5.14) 
aP = af, 


et de même pour les b;. C’est la philosophie qu’on retrouve dans la théorie des perturba- 
tion de la mécanique quantique par exemple. 


Méthode itérative 


La méthode suppose que la séquence suivante converge. À partir de 


à; = {ai, Hı} (5.15) 
b = {b Hi}, (5.16) 
on calcule d’abord la première itération 
à = {a,H}lo,o (5.17) 
a®) = {ai, Hi} pO (5.18) 
du + Jha (0) (0) signifie que le résultat du calcul du crochet est évalué en af ) DA 


i.e. que les a® et aO apparissa à droite de l’équation après le calcul du crochet sont 
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remplacés par les déjà connus. Suite à l’intégration de Ho. La seconde approximation 
suit 


à = {a ÆH}loo (5.19) 


aP = {a,H},oso (5.20) 


etc... 


Méthode de la moyenne 


Elle est utilisable lorsque la solution non perturbée est une orbite cyclique de période 
. On peut alors calculer l’effet net moyen de la perturbation sur une orbite par 


1 F 
di = — | {ai, Hı }dt (5.21) 
T Jo 


1 T 
bi = J {b;, Hı }dt (5.22) 
T Jo 


Cette méthode est compatible avec la méthode itérative par exemple. Les changements 
orbitaux des satellites et des planètes dus à certaines excentricités ou aux autres planètes, 
sont généralement calculés de cette façon, comme le déplacement (rotation) de l’ orbite 
(elliptique) de Mercure par exemple. 


Remarque 7 

Les crochets de Poisson qui apparaissent ici sont présumés calculés en utilisant les va- 
£D riables canoniques q; et pi. Cela implique d’avoir constamment recours aux équations (i) 

et (ii). On verra en plus bas une façon systématique de choisir ces constantes en optant 

pour les a; et B; de Hamilton-Jacobi qui ont l’avantage considérable d’être variables 

canoniques. 


5.4 Exemple 


Exemple 5.1 
Voyons d’abord un cas très simple, soit celui d’une particule soumise à une force constante en une 
dimension donc A 
Pp 
H(p) = Z+ (5.23) 
où V = àq donc la force s'oppose au mouvement vers les q croissant si À > 0 et l'inverse si 
À < 0. On sait résoudre exactement 
: P i 
ġ4={4 H} =>, p= {p, H} = -A (5.24) 
et donc 


„à q(t) = qo + dot — À 
Te = { p(t) = må = mûo ÅN (aaa 
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Afin de tester la méthode perturbative décomposons H en Ho + Hı où 


H=, Hi= (5.26) 
m 
et reprenons les étapes i) - v). 


i) La solution analytique pour H est triviale, nous avons une particule libre et donc 


q = a°)++bp0) (5.27) 
p = ma). (5.28) 
ii) L’inverse de ces équations est 
an = 2 (5.29) 
m 
po = q- t (5.30) 
m 


et même si b®) dépend explicitement du temps on vérifie (c’est inutile en fait) que 


O) 2 140) 0, = 
a = {a , H} + = ET ={Ż2, Pit 0 (5.31) 
po) = w, H} po 2yo 
= D} P? 
= u- Ta 
— P 
= E 
2 
= Peter Pets) (5.32) 
2m m m m 
iii) On pose que la solution pour H sera 
q=at+b a= À 
sam R T (5.33) 
iv) On calcule les équations d'évolution de a et b 
: P 
= H}=1{— 
à = {aH}={2,x) 
À À 
= 2ipa} =À. (1 
T (5.34) 
m 
pe Her 
b = {b, Hi} = {4 m AD 
À At 
= zina- z Pa} 
At 
eu) Pa (5.35) 
m m 


v) On intègre trivialement en posant que la perturbation a été allumée à t = 0 et que pour t < 0, 
seul Ho jouait un rôle. Ceci nous donne les conditions initiales donc les constantes d'intégration 
pour les équations pour a et b. Donc ici 


af) = -tta = a(0) = a (5.36) 
b(t) àp Hb% — b(0) = b., (5.37) 
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Remplaçant dans les équations en iii) on obtient 


2 
= («° E x) tq A pO (5.38) 
m 2m 
2 
= BO pa AE (5.39) 
2m 
ainsi que 
p = ma) — At. (5.40) 


Ces résultats sont exacts avec qo = bO eta® = do- 


5.5 Méthode canonique de perturbations 


L'idée est la même mais les manipulations sont sensiblement simplifiées du fait qu’on 
choisit les constantes de la méthode de Hamilton-Jacobi, les œ; et 5; au lieu de laisser ce 
choix au hasard. Ces constantes ne seront vraiment constantes que pour Ho, l’introduc- 
tion de H; fera qu’elles ne seront plus constantes du mouvement. L'avantage vient du 
fait que les à; et 8; étant variables canoniques (respectivement moments et coordonnées 
généralisés) on peut les utiliser pour calcules les crochets de Poisson. On évite ainsi cet 
incessant va-et-vient entre les (a;,b;) et les (p;, qi) qui ressort dans les exemples de la 
section ci-dessus. Une fois que Ho a été résolu par H.-J. on obtient 


qi = dt,a;,B;) (5.41) 
Pi = pilt, aj, B3) (5.42) 
que l’on inverse en 
Bi = Bi(t,g;,p;) (5.44) 
C’est toutefois la première forme qui est utile puisqu’elle permet d’ écrire 
Hı (qi, pi) = Hi(qi(t, aj, B4), pilt, aj, B;)) (5.45) 
= Ki(aj, bjt) (5.46) 
Par la suite, tous les calculs des crochets de Poisson se feront par 
ðA B ðA ðB 
{4, B} = FE — E (5.47) 
2 0B;0a; aj ðb; 
ainsi 
: Z Oai 0H: Oai 0H: 
Qi = {ai, Hı} = | a A (5.48) 
2 bj ðaj ða; ðb; 
— ôK, OK 
= a Ôi; = (5.49) 
2595, = 08, 
et 3K 
To (5.50) 
Oai 


C’est sous cette forme que la théorie des perturbations est généralement présentée 
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dans la littérature. 
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© 


Exemple 5.2 
Voyons un exemple assez ”classique” parfois appelé l’oscillateur quantique décrit par 
2 2 
p mw 2 mk 4 
H = — — q". 5.51 
ET (5.51) 
Les équations canoniques du mouvement sont 
H 
ds Met (5.52) 
Op m 
H 
D = = mu? q — mkg? (5.53) 
ôq 
ou, en les combinant 
ÿ = —w°q — kg. (5.54) 
Intégrer cette équation n’est pas trivial. Choisissons de décomposer H en Ho + Hı où 
2 2 
Ho = Po m e : oscillateur harmonique (5.55) 
2m 2 
Hı, = ur, : perturbation. (5.56) 


Nous allons d’abord résoudre pour Ho par la méthode de H.-J. Puisqu'on peut écrire, Ho étant 
indépendant du temps, 


S(q,a,t) = —at + W (q,a) (5.57) 
et sachant que 
ðs os 
Ho(q, —)+— = 0. 5.58 
o(q, ôq ) + ðt ( ) 
On voit immédiatement que 
oW 
Ho(q, ET =0. (5.59) 
Ainsi notre nouveau moment canonique, à, sera une constante égale à l'énergie du système! Ex- 
plicitant l équation ci-dessus 
1 dWN\? ma? 2 
2m (=) + 2 q — (5.60) 
on obtient aw 
— = y 2ma — m2w2q? (5.61) 
dq 
W = faw 2ma — M?w?g?. (5.62) 
En plus de p — nos équations de transformation canonique comptent 
OS dq 
s = Zr |- (5.63) 
q 2ma — M?w?g? 
=) pag | net | (5.64) 
w 2ma 


qui s’inverse en 


q= 1/2 sino(t + 9). (5.65) 
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Nous pouvons récrire Hı 


RS 5.66 
1571 re w(t + 8). (5.66) 


a et 5 nos nouveaux moment et coordonnée généralisés sont des constantes sous Ho mais ne le sont 
plus lorsqu'on introduit la perturbation (disons à t = 0). Leur équation d'évolution est canonique 


0H: _ 4ko? 


à = — DE — us Sin w(t + 6) cosw(t + 8) (5.67) 
b = ʻi - -2a sint w(t + 8) (5.68) 


Ces équations ne sont pas triviales à résoudre non plus mais elles se prêtent à une d’approximation, 
que ce soit par développement en série, par itération ou par moyenne. 


Développement en série 


On constate que k est le paramètre qui caractérise H4. Ce n’est par un très bon choix 
puisqv’il est lui-même dimensionné, néanmoins nous allons tenter une expansion en série 
du type 


a = ao +ka + k’as+--. (5.69) 
B = Bot kpi +k’ Bat. (5.70) 
Cependant on constate ici que le côté droit des équations pour à et b dépend de fonc- 
tions trigonométriques d’argument w(t + 8). Comme les fonctions trigonométriques sont 


hautement non-linéaires dans leur argument, il n’est pas trivial d'identifier leur degré de 
dépendance en 8. Par exemple 


B ~ sint w(t + 8) (5.71) 
de 
sinw(t + 8) = sin wt cos wp + sin w8 cos wt (5.72) 
on aura 
coswð œ~ cos(wbs +wkb:) 
cos wbo cos wkp — sin wbo sin wk p1 (5.73) 
sinw ~ sin(wbo + wk) 
~ sinwbycoswkb + cos wbo sin wk bı (5.74) 
et prenant 
sinwkB, ~ wkßı+ø0(k?) (5.75) 
coswkB, ~ 1+O(k?) (5.76) 
alors 
sinw(t + 6) = sinw(t + Bo) + wk: cosw(t + Bo). (5.77) 


Si k est considéré petit alors 
sinw (t+ 8) & sinwt(t+ bo) 
+4wkb: sin? w(t + bo) cos w(t + Bo) + O(k?) 
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(5.78) 


et faire les remplacements appropriés pour identifier les termes d’une puissance donnée 
de k. L'exercice est assez lourd ici et nous ne le complétons pas. 


Solution itérative. 


À partir des équations pour å et b qui sont de la forme 


à = f(a,B,t) (5.79) 
B = g(a,B,t) (5.80) 
elle consiste à dire qu’on peut tendre vers a et 8 par une série d’itérations 
f oH 
å = fond (5.81) 
a=aQn_1,/08=8 4 1 
: 0H 
b = gant) = 5 (5.82) 
a &=Qn—1 =bn 
au sens où 
lim an > a (5.83) 
n— 00 
lim 8, > 8. (5.84) 
n— oo 
La première de ces itérations est 
; 0H 
àr = flao Bot) =- 7 (5.85) 
a=00,8=69 
: OH 
bi = go bot) = . (5.86) 
Q |,,— a 
a=a0,8=f0 
Nous nous limiterons ici à cette première étape et résoudre donc 
i 4ka? . à 
å = =z sin w(t + Bo) cosw(t + Bo) (5.87) 
. 2k 
b = -2 sint w(t + po). (5.88) 
Comme ao et g ne dépendent pas de t, l’intégration est triviale et donne 
k 2 
ar = = sint w(t + Bo) + C (5.89) 
mw 
2kao [3w sin 2w(t + Bo) 
= — (t PU? 
Fi mwt { 8 (E+ Bo) 4 w 
in 4w(t 
et) Lo (5.90) 


Les constantes d’intégration sont ajustées par les conditions initiales touchant la pertur- 
bation. Si par exemple, on dit que la perturbation est allumée à t = 0, alors pour t < 0, 
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la solution non perturbée prévaut et a(t < 0) = @o, O(t < 0) = Bo. Donc à t = 0, 
&1(0) = ao, 81(0) = Bo, ce qui fixe 


ka2 
C = ao — À sint wbo (5.91) 
mw 
2kao | 3wbo sin2wf, sin 4wpo 
1 
Sp SM UD y SU 92 
7 A ur 1 a a 
et alors 
2 
aÿ r. ; 
ai = ao- a [sinf w(t + 86) + sinf wbo] (5.93) 
B = B 2kao f3wt  sin2w(t + bo) — sin 2wbo 
1 0 mwt | 8 4 
sin 4w(t + Po) — sin 4wbo } | (5.94) 


Ci-dessous fixons 8g pour alléger les expressions. La solution perturbée est obtenue de 
celle non perturbée en y remplaçant a: et 84. Nous obtenons 


2 


2 l 
2 kaĝ . 4 
—; [00 — 7 Sin wt| x 
mu MW 


2kao E _ sin2wt sin =) 


q(t) 


8 4 F7 32 


sin (a + 5 
mw 
(5.95) 


On voit que l’amplitude, qui était (22%)? dans le cas non perturbé, a été modifié, œı 
étant remplacé par 


ka . 
feo — —© sin st (5.96) 
Mw 
ceci nous donne un test de petitesse de k puisque la quantité est à la puissance ł et que 
q(t) doit demeurer réel, donc 
mw 
k< —. (5.97) 
ao 
Pour k > 0, l’amplitude décrofît. De plus, le comportement qui était harmonique en w, 
i.e. sin wt a été modifié en 


| 3kao 2kao sin2wt  sin4wt 
1+ = jz = 5.98 
sin (w ( +) Fan | 4 Li 32 |) TP 
où, si on veut encore parler d’une fréquence 9, on doit définir 
3k 
nie. (5.99) 
Amuwÿ 


On voit que pour k > 0, la fréquence augmente. Il faut aussi préciser que s’ajoute une 
modulation en sin 2wt et sin Awt. Strictement le mouvement n’est plus harmonique. 


Méthode de la moyenne 


Elle ne donne pas des résultats aussi détaillés que cette en série mais c’est parfois 
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suffisant. Sachant que le mouvement non perturbé est ici cyclique de période (T = 2, 
ici) nous remplaçons à et b par & et B moyennés sur une période 
= 4k 1 f 
ài = -5 / a? sin? w(t + B) cosw(t + B)dt (5.100) 
MW? T 0 
= 2k 1 [f 
b = ns | asint w(t + Bdt. (5.101) 
MOT 0 


Il est trivial de voir que, ne connaissant pas a(t) et B(t) (c’est ce que nous cherchons), 
il est difficile sinon impossible de faire les intégrales. C’est pourquoi cette méthode est 
souvent augmentée de l’approximation itérative, remplaçant a et 5 dans les intégrales 


par a et Bo pour faire un premier calcul de & et b, le résultat duquel pourra être remplacé 
dans les intégrales...etc. Au premier ordre nous aurons ici 


27 
= 4k 1 Tu 
à x rc | sin? w(t + B) cosw(t + B)dt 
0 
4k w ; T=2r 
z= — sn 20 sinf w(t + 8)|, ” =0 (5.102) 
donc & ~ 0 — a = constante = &o 
27 
= 2k 1 TTo 
B x agen f sinf w(t + B)dt 
— 2kao w [3w(t+6bo) sin2w(t+ bo) + sin 4w(t + Bo) =i 
~ m2r 8 4 32 5 
— 2r 
k 3w 2 TTo 83k 
a $z + o] g (5.103) 
TMW 8 w o Mw 
et donc m 
a 
HOR i! + bo: (5.104) 
Ici, la solution perturbée se lira 
2ao . 3ka 
qt) 7x n sin (a + Bo + rent) 
20 g 3k@o 
1+ —> |t ; 1 
mz Sin (2 | + pres + bo (5.105) 
Le seul effet de la perturbation ici est une modification de la fréquence qui de w passe à 
3k 
ul) (5.106) 
AmuwS 


donc qui augmente si k > 0, et qui est d’ailleurs la fréquence Q déjà obtenue. 


Remarque 8 

À une énergie E donnée (voir figure 5.1), le mouvement va de —x0 à +0 avec une 
fréquence w. Introduisant le terme quartique diminuera l’amplitude entre —x1 à +æ: 
tout en affectant le fréquence de w > Q.. 
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Figure 5.1 


5.6 Autre exemple 97 
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pm memes, 
6 MOUVEMENT DU SOLIDE 


6.1 Degrés de liberté du solide 


Jusqu'ici nous n’avons considéré que des particules ponctuelles en nombre relative- 
ment petit. Nous allons maintenant permettre aux corps physiques d’avoir de véritables 
dimensions physiques, telles des longueurs, largeurs et épaisseurs. Nous allons cepen- 
dant nous limiter aux corps indéformables, ce qui est une approximation de la réalité 
physique, mais une approximation souvent très valables. Nous considérerons donc que 
chaque point du corps solide demeure à distance constante de tout autre point du même 
corps solide. Il n’y a pas de déformation. 

En mécanique classique la structure microscopique du corps solide est sans intérêt. 
On peut donc le considéré comme constitué d’un grand nombre de petites particules ou 
comme un ensemble continu de matière. Par exemple, la masse d’un tel corps s’écrira 


M = Dm = L p(x)d?x (6.1) 


où m; serait la masse de la particule à, constituante du corps solide, et p(x) serait une 
densité continue de masse (les unités de masse par volume). À l’occasion nous utiliserons 
donc l’une ou l’autre notation. La notation discrète (5), ) est parfois plus pédagogique 
puisqu'elle fait essentiellement la somme sur un grand nombre de particules ponctuelles, 
concept avec lequel nous sommes maintenant familiers. 

Dans l’espace physique à trois dimensions, une particule ponctuelle a trois degrés de 
liberté. Dans le cas du corps rigide on se convainc rapidement que l’état du solide peut se 
décrire par la position d’un des points du solide (3 degrés de liberté) et l’orientation du 
corps rigide (solide) par rapport à un système d’axes fixées en ce point, ce qui implique 
trois autres degrés de liberté. Au total donc un solide a 6 degrés de liberté. (Il est très 
avantageux, pour assurer la simplicité des expressions qui vont suivre, de choisir le point 
dont nous suivons le déplacement et par rapport auquel nous mesurons l’orientation du 
solide, comme étant le centre de masse du solide). Si on visualise le solide comme étant 
constitué de N particules, on devrait avoir à priori 3N degrés de liberté. Le fait qu’il n’en 
reste que 6, résulte de l’ensemble de contraintes qui font que chacune de ces particules 
est toujours à égale distance de chacune des autres. 

Pour décrire ces 6 degrés de liberté nous procéderons de la façon suivante. Imaginons 
un premier système de référence, noté XY Z, fixé dans le laboratoire et présumé inertiel. 
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Fixons ensuite rigidement au corps en O un nouveau système de référence £1 £2£3 (voir 
figure 6.1). Ce système se déplace et tourne avec le corps rigide. Ce n’est donc généra- 
lement pas un système inertiel mais comme il est solidaire du solide il apparaît comme 
immobile à un observateur se trouvant sur le solide. Pour faciliter le travail à venir nous 
centrerons souvent le système x1x2æ%2 sur le centre de masse du solide, auquel car R 
est la position du centre de masse du solide. r est la position d’un point P de ce solide, 
mesurée dans le système inertiel XY Z. Notons par x la position de ce même point P 
mesurée dans le système x122t3. Comme ce dernier est fixé au corps, x est constant 


pr 
N 


X 


Figure 6.1 


(mesuré dans z1x2£3). Lorsque le solide se déplace, le point P se déplace. Pour mesu- 
rer un déplacement dans le système inertiel, dr, nous le décomposons en déplacement 
du point O, dR, plus un changement possible d’orientation du solide, dy , mesuré par 
rapport à l’axe de rotation instantanée. Alors 


dr=dR+dpxx où dy =ñdọ (6.2) 
ce qui nous donne 
._ d__ č dR dẹ _ 
verca ug y VTO (6.3) 
ou 
v=V+AOxx. (6.4) 


Si R mesure la position du C.M., i.e. si O est positionné sur le C.M., alors V est la 
vitesse de C.M. et correspond à une translation du solide comme un tout. Q est la vitesse 
angulaire du solide et sa direction, comme celle de d@ coïncide avec l’axe de rotation 
du solide. Notons qu’elle n’est pas constante en général. Comme le système z1%2%3 est 
fixé dans le solide, Q est également la vitesse angulaire de la rotation de ce système. Ce 
résultat ne dépend en aucune façon du fait que nous ayons centré le système z1%2%3 en 
O, le C.M. du solide. Nous aurions pu choisir ici un autre centre O’ déplacé de O par 
une longueur a au sens où la position x’ du même point P est reliée à x par 


x=x +a. (6.5) 
Nous aurions alors au lieu de (6.4) 
v=V+Qx(x +a)=V+9Q9xx+Q9xa. (6.6) 
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D'autre part, à partir de la forme de (6.4) nous pouvons écrire 


v=V'+0 xx (6.7) 
ce qui nous force à identifier, v étant identique à lui-même et a arbitraire 
V'=V+AOxx, Q =Q. (6.8) 


La deuxième de ces équations est importante puisqv’elle nous indique que la vitesse an- 
gulaire de rotation est totalement indépendante du système £1 £2%3 (fixé dans le solide) 
choisi. À un moment donné, tous ces systèmes tournent donc autour d’axes parallèles les 
uns aux autres (de direction donnée par celle de Q) avec une même vitesse angulaire Q. 
Cette propriété d’absolu dans la rotation du solide ne se retrouve pas dans la translation 
du solide puisque V’ # V. Notons qu’en général, lorsque le solide se déplace, Q n’est 
constant ni en direction ni en longueur. Il est parfois intéressant de choisir une origine 
O' telle que V” = 0. Instantanément le mouvement apparaîtra comme une rotation pure 
autour de l’axe défini par passant par O’ évidemment. On appelle cet axe, l’axe de rota- 
tion instantané du corps. Cependant à partir de maintenant nous choisirons l’origine O 
du système zı z2£3 comme étant le centre de masse du solide, à moins que la chose ne 
soit clairement spécifiée. 


6.2 L'énergie cinétique et le tenseur d'inertie 


Considérons le solide comme étant constitué de points matériels discrets et calculons 
l’énergie cinétique du solide, évidemment mesurée dans le système inertiel XY Z 


1 
T=} > mv? (6.9) 
part. 
la somme porte sur tous les points du solide A T v2 mais pour simplifier l’ écriture 


nous laissons tomber les indices identifiant ces points. Nous nous sommes évidemment 
placés dans le référentiel inertiel et v? = v? où v est défini par l'équation (6.4), ce qui 


donne 
part. 
= SX mV? +Y mV (Q xx)+5 ð m(Q x x) 
part. part. part. 


(6.10) 


Les vitesses V et Q sont les mêmes pour tous les points et peuvent donc sortir des 
sommes. Ainsi le premier terme devient 


lm É Dm = (6.11) 


part. part. 


où M = masse totale du solide. 
Le deuxième terme devient, par la propriété des produits triples 


Xmv. (D XX) = Yms (V x Q) 


part. part. 
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= (V xQ)-ÿ mx=0 (6.12) 
part. 
parce que la somme est nulle, ayant choisi l’origine du référentiel x1æ2x3 au centre de 
masse du solide. Pour le troisième terme nous développons le carré du produit vectoriel 
1 1 
Dm xx) =A m|- a.x] . (6.13) 


part. part. 


Au total (6.10) devient donc, lorsque (£1 £23) est centré sur le C.M., 


2 
PES [2x - (Q. x)’ , (6.14) 
part. 


Le premier terme est l’énergie cinétique de translation du solide. Ce serait le seul terme 
si toute la masse du solide était concentrée au C.M. Le deuxième terme est l’énergie 
cinétique de rotation, donc 

T = Tem + Trot. (6.15) 
Étudions la forme de To en décomposant les vecteurs selon les axes x117273 puisque 
Tot est une énergie cinétique de rotation impliquant certains concepts d’inertie qui serait 
propres i.e. intrinsèques au solide. Les composantes de Q et de x seront donc selon les 
axes du référentiel O£x1£2x3 (notons que T demeure mesuré dans le système inertiel) 


1 
Tiot -= 5 >" (Omar — Dit] 
1 
part. 


Les coordonnées? x; dépendent de la particule sur laquelle on fait la somme mais Q; 9y 
n’en dépend pas et peut sortir de la somme sur les particules, ce qui donne 


1 
Tot = 50 Sn LOiktiti — tit] 
part. 
1 1 iag 


Cette expression a la forme usuelle d’une énergie cinétique mais ce qui le rôle d’inertie 
est plus compliqué que dans le cas des translations. Ici on l’appelle le tenseur d’inertie, 
on identifie ses éléments 
Iik = Xom [£i£iðik — vite] = Iri. (6.18) 
part. 
On dit qu’il est un tenseur parce qu’il a deux indices. On peut dès lors lui donner une 
représentation matricielle 


hi h2 hs 
I = lı Iz L3 ; (6.19) 
Isı I32 J33 


2 Dans ce chapitre, nous utilisons la notation d’Einstein où un indice répété dans un terme est automatique- 


ment sommé à moins d’avis contraire, ainsi Q; Qiziz = J; ji DQrORkRTITI 
ija 
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Si on écrit aussi 


Qı 
Q=| % | = QT = (Q, R2, Q3) (6.20) 
Q3 
alors on peut écrire 
1 
Trot = 52710 (6.21) 
À partir de (6.18) on calcule directement 
Emr? +r) -P mr — Ý mrits 
I= -Ð mrt maxi +r) -P mrr3 (6.22) 
— D masi -Ð mzz Y` m(x? + x) 


où ici, on a alléger la notation par n — J. On voit aussi qu’on peut passer à une 
rotation et un calcul continues où la matière est réputée être distribuée de façon continue 
dans le solide selon une densité p(x) = p(x1, £2, £3). On écrit alors 


Lx = 1 p(X) [£i£iĝik — titre] dridr2dt3 (6.23) 
v 


Il est possible de choisir le référentiel Oxı £2£3en l’orientant de telle sorte que T est 
diagonal 


Ir 0 0 
I=| 0 w 0 (6.24) 
0 0 Iy 


Les éléments Jı; Iy et I% sont en fait les valeurs propres de la matrice (6.22). Ici nous 
avons noté prime (/) le référentiel (voir figure 6.2) qui garantit que le tenseur d’inertie 
est diagonal: Oxz’ £y £3. On appelle ces trois directions, Oz1/, Ox» et Ox3 les 


\ OON 


D 


Figure 6.2 


axes principaux (d’inertie) du solide et Jı; Ixy et T3 les moments principaux d’inertie. 
Si on choisit le différentiel Oxı£2£3 pour coïncider avec les axes principaux alors Trot 
devient 


1 
Tiot = AUS (6.25) 
Pour alléger la notation il sera entendu dans ce qui suit que, lorsque les éléments du 


tenseur d’inertie apparaissent avec un seul indice, c’est que nous aurons choisi de faire 
coïncider le référentiel fixé au corps et les axes principaux. Nous laisserons tomber les 
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primes (”) pour écrire simplement 


Tio = ERM. (6.26) 


6.3 Parenthèse sur les axes principaux et le tenseur 
J 
d'inertie 
| 
Il n’est pas nécessaire de diviser à priori la direction des axes principaux. Lorsque le 
solide a certaines symétries la direction de ces axes est parfois évidente. Il est toujours 
possible de choisir arbitrairement un système Ox1x2x3 et de déterminer par rapport à ce 
dernier la direction des axes principaux. 
On calcule d’abord les éléments de I par rapport au système d’axes choisi, ce qui 
nous donne J;k qu’on retrouve en (6.18), (6.22) ou (6.23). On remarque d’abord que la 
matrice (6.22) est symétrique 1.e. 
lik = Iki; lik réel (6.27) 
Les valeurs propres de cette matrice seront tout simplement les éléments J4, To et 13 de 
I sous sa forme diagonale, Ip. 
De façon plus technique nous dirons qu’il existe une matrice U, avec son inverse 
U- 1, telle que UT! = Ut : unitaire 
UIU! = Ip: diagonale. (6.28) 


Reprenant l’ expression pour Tot 
1 1 
To = z929 = SAUT UIU UQ 


= SOU UN = op (6.29) 


ùQ = UR. 
Strictement ceci termine l’ opération puisqu’en (6.29) nous avons Tot écrit en utilisant 
Ip. On voit qu’ici 


Qi Ui u2 u3 Qı 
Q =UQ — A = U21 U22 U23 Qə (6.30) 
3 U31 U32 U33 Q3 
ou encore 
QI = un + uN + u13Q3 
(ea =  U21Q1 + U2202 + U23983 (6.31) 
Q% = us Qı + u3z2NQ2 + u33 Q3. 


Le système prime qui définit les axes principaux du solide est obtenu du système original 
par une rotation du système original à condition que l’origine O ait été choisie comme 
le C.M. du solide (voir figure 6.3). Autrement il faudra effectuer d’abord une translation 
vers le C.M. Évidemment le vecteur ne bouge pas lors de cette rotation qui n’est en fait 
qu’un simple réalignement des axes du référentiel fixé dans le solide. Il n’a rien à voir 
avec le mouvement de rotation du solide. 
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Figure 6.3 
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Techniquement tout repose sur la matrice U, ce qui est relativement simple puisque 
les colonnes de U! sont tout simplement les vecteurs propres de T. Ceci se vérifie 
immédiatement en rappelant (6.28) 


UIU™ = Ip (6.32) 
qu’on multiplie par la gauche par UT! ce qui donne 
IU! = U™!Ip. (6.33) 


Il s’agit d’une égalité entre 2 matrices 3 x 3. Écrivant vij pour les éléments U7! nous 
explicitons (6.33) mais en ne spécifiant ici que la première colonne des deux matrices 
produites ce qui donne 


Davy + hova + lgti +++ ++: lvi . 
aivi + D22022 + I23v32 +++ ++ | = | vas ++. +. (6.34) 
J31v13 + 32023 + 33033 ++ oee Liv33 f 


Si deux matrices sont égales c’est que tous leurs éléments sont égaux et par extension 
les éléments d’une colonne de l’une sont égaux aux éléments d’une colonne de l’autre. 
Considérons la première colonne de chacune des deux matrices ci-dessus et égalons l’une 
à l’autre. On constate immédiatement qu’il est possible d’écrire l’égalité entre ces deux 
colonnes sous la forme 


hi h2 hs v11 v11 
lı Iz z vai | = 11 | va |. (6.35) 
Isı I32 T33 v31 v31 


Si nous écrivons VA pour la 1ère colonne de UT! i.e. (ne pas confondre VA avec une 
composante du vecteur vitesse) 


v11 
V= | va (6.36) 
U31 
cette équation s’écrit 
IV =V; JT = un nombre (6.37) 


où J; est la 1ère valeur propre de la matrice I alors que Vi est un vecteur colonne. C’est 
l’équation type du problème aux valeurs propres. Pour résoudre on obtient d’abord les 
valeurs propres de T et ensuite on obtient les éléments (non normalisés) de VA à l’aide 
de (6.35) ou (6.37). Nous aurons ici 3 équations de ce type à partir de (6.34), une pour 


chaque valeur propre 1; avec son vecteur propre V; qui constitue la 4 colonne de U-!. 


La séquence d’opérations est donc la suivante. 


1. On choisit un référentiel centré sur le C.M., Ox1x2x3 par rapport auquel on calcule 
les L;;; i, j = 1, 2, 3, éléments de T. 


2. On calcule les valeurs propre de cette matrice, ce qui nous donne Ip qui n’a que les 
éléments diagonaux I4, L et I3. 


3. On calcule les vecteurs propres V4 de I, chacun correspondant à une des 3 valeurs 
propres, J1, I> et I3. 


4. Les Vp sont les colonnes de la matrice U~! que nous pouvons inverser pour avoir la 
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matrice U, une matrice qui représente une rotation. 


5. Pour obtenir les axes principaux Ox1x2%3 il suffit de soumettre le référentiel Ox1172t3 
à la rotation représentée par U. 


Remarque 9 

La rotation représentée par U est généralement par rapport à un axe qui n’est pas un des 
axes de Oxit2t3 ni de Oxx tx £z. Il est cependant toujours possible d'opérer cette 
rotation à l’aide de 3 rotations successives faites autour d’axes choisis. La façon la plus 
courante est celle des angles d’Euler 


Laissant tomber les prime (”) et supposant que les axes du système Ox1x2%3 coïnci- 
dent avec les axes principaux la définition générale des éléments de J;; en (6.18,6.22,6.23) 
est toujours valide sauf que seuls les termes j = i ne seront pas nuls, nous les avons notés 
avec un seul indice 

I = In = X man — z?) (6.38) 
part. 
où il y a une somme sur l mais pas sur i. Ainsi trivialement 


Il = Ÿ_m(x3 + 23) 
part. 

b = X m? +z?) (6.39) 
part. 

B = $ maita) 
part. 


et on note qu'aucun des ces J; n’est plus grand que la somme des deux autres; tout au 
plus est-il égal à cette somme. Lorsque Jı = I> = 13 on dit que nous avons une toupie 
sphérique. Si deux seulement des moments d’inertie sont égaux, on parle d’une toupie 
symétrique et si les trois sont différents, d’une toupie asymétrique. 


Une remarque importante s’impose ici. Nous avons réussi en (6.14) et (6.17) à écrire 
1 1 
T = 5 MV° + SO (6.40) 


expression qui n’est valide que si le référentiel intrinsèque est centré sur le centre de 
masse, O. Cependant, pour calculer les J;x, il peut s’avérer utile d’utiliser d’abord un 
autre référentiel, également intrinsèque mais centré sur une autre origine O’ et dont les 
axes sont parallèles au premier. Il s’agit donc ici d’une translation du référentiel et non 
d’une rotation (voir figure 6.4). Appelons a le déplacement OO’, de telle sorte que 


x= X +a =T; = T; + ai. (6.41) 
On sait que 
Li = 5 mM(z£i£iĝik — tik). (6.42) 
part. 


Par rapport au système prime nous aurons 


f 1 R À EE À 
ik 7 ` MT ir — Li£p) 


part. 
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= Xom [Car + ar) (£1 + ai) ik — (zi + ai) (£k + ax)] 


part. 
= 5 m [£i£iðik — tite] + 2a1fik 5 ma 
part. part. 
+ai X maer — ak X mr + (aiaiĝik — aiak) Xom 
part. part. part. 
= Lg + Maud — a;ax) (6.43) 
OÙ D part. MTi = 0 et Dan, M = M. C’est le fameux théorème des axes parallèles. 
å 
4 3 
%3 
A 
Z 
x 7} x' 
ae i 
Aa dé 
A a 
at ‘ 
Figure 6.4 


Nous savons donc écrire l’énergie cinétique du solide en (6.40). Si le référentiel int- 
rinsèque correspond aux axes principaux alors cette expression se réduit à 


T= I MV? + ZRO. (6.44) 


Deux types de forces peuvent être présentes dans le système; des forces de cohésion 
particule-particule dans le solide dont le résultat global sur le solide est nul à cause du 
principe d’action-réaction comme nous l’avons vu au chapitre I. Il reste les forces ex- 
ternes et si ces forces sont dérivables d’un potentiel U alors on peut écrire le Lagrangien 


L= imv? + K — U. (6.45) 


Nous y reviendrons plus tard. 


6.4 Le moment cinétique/angulaire du solide 


Le moment cinétique dépend du point par rapport auquel il est défini. Dans l’ étude du 
mouvement du solide il apparaît raisonnable de choisir ce point à l’origine du référentiel 
intrinsèque qu’en (6.45) nous avons choisi pour coincider avec le C.M. du solide (voir 
figure 6.5). Nous avons noté x la position d’un point P mesurée à partir de O. Le moment 
cinétique de ce point matériel est 


lp = mx x v. (6.46) 
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Ici, la vitesse v est uniquement celle due à la rotation du solide, rotation que se fait 


Figure 6.5 


à la vitesse de rotation instantanée et par conséquent v est ici v = Q x x qui ne permet 
pas de changer la longueur de x comme il se doit puisque nous avons un solide rigide et 
donc, sommant sur tous les points matériels nous aurons pour le moment cinétique 


l= X mxx (Q xx). (6.47) 
part. 


Explicitant chaque composante du triple produit vectoriel nous avons 


li = XO m (21210 — mire Q) 


part. 
= Xm (£1£1ikNk = rt) 
part. 
= Qk 5 m (£ı£ıðik — TiTk) (6.48) 
part. 


et donc 

L= ipp (6.49) 
Evidemment si nous avions choisi de faire coincider le référentiel intrinsèque avec les 
axes propres du solide nous aurions simplement 


L=, b=bLQ, l3 = BN. (6.50) 


On voit donc qu’en général (sauf pour une toupie sphérique) la direction de 1 ne cor- 
respond pas à celle de Q. C’est là une différence dramatique avec le mouvement d’une 
particule et ce seul fait introduit déjà des différences notables entre le mouvement de la 
particule et celui du solide. Étudions brièvement la situation qui prévaux dans trois cas 
simples sans forces extérieures, i.e. le Lagrangien se limite à l’ énergie cinétique. Dans un 
tel cas on peut faire coincider les origines des référentiels inertiel et intrinsèque puisque 
le C.M. ne sera soumis à aucune accélération: V = 0. 


Exemple 6.1 


La toupie sphérique (essentiellement une sphère) a ses trois moments égaux Iı = 12 = 13 = I et 
par conséquent 


1 =IQ. (6.51) 
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En l'absence de force/torque extérieur l est une constante. Il en va de même de Q. La toupie 
sphérique libre tourne tout simplement par rapport à un axe fixe défini par l ou Q (le même axe) 
à vitesse constante . 


Exemple 6.2 


Un autre cas simple est le rotateur où la masse du solide est essentiellement répartie sur une droite. 
Plaçant les axes principaux comme sur la figure 6.6 on voit que I3 = 0, Iı = 12 = 1. De plus 


Figure 6.6 


Q3 = 0 puisque la rotation d’une droite sur son axe est comme la rotation d’un point et wa pas 
de sens (du moins classiquement). La vitesse de rotation n'a donc que des composantes Qı et Q2 
et Q se trouve dans le plan xı Oxz comme d'ailleurs l et Q sont colinéaires comme dans le cas de 
la toupie sphérique. 


Exemple 6.3 


La toupie symétrique (voir figure 6.7) est un véritable objet en trois dimensions (un beigne, un 
ballon de football sont des toupies symétrie) caractérisé par un axe de symétrie que nous choisis- 
sons comme Ox3. Ainsi Iı = 12 # Is avec Li, Iz et Iz Æ 0. En l'absence de torque i.e. en 
rotation libre l = constante. Pour décrire qualitativement le mouvement nous figeons le temps au 
moment où le plan 1Ox3 est perpendiculaire à l'axe Ox2 i.e. correspond au plan x1Ox3. À ce mo- 
ment lə = 0 mais puisque l2 = 12Q2 nous avons Q2 = 0. Donc le vecteur Q est alors également 
dans le plan 1Ox3.La propriété que 1, Q et Ox3 sont dans le même plan a été obtenue facilement 
à la suite d’un choix particulier d'orientation de l'axe Oxz mais une propriété indépendante de 
ce choix et reste vraie pour tout le mouvement. Ainsi tout point sur l’axe de symétrie Ox3, iden- 
tifié par Xa a une vitesse donnée par ® X X, qui sera nécessairement perpendiculaire au plan 
10 Oz3 et puisque l est constant, en longueur et en direction le mouvement sera forcément tel que 
laxe de symétrie tournera autour de la direction donnée par 1, laxe Ox3, ce dernier dessinant 
un cône autour de la direction constante, 1. Ce mouvement est appelé précession naturelle de la 
toupie symétrique. Le mouvement de la toupie se décompose donc en rotation de la toupie autour 
de son axe Oxs plus la précession autour de l, on décompose Q, qui est dans le plan 10%: ,selon 
Oxs et selon 1 qui ne sont en général pas orthogonaux. Pour faire le calcul, on se replace au 
moment où Ox: est dans le plan IQ. Selon la figure 6.8, clairement 


— Qı e vi Le 
cos x a, cos(> 0) = sin@ (6.52) 


donc o 
Q = Q,sin0 = Q,, — 2 (6.53) 
sin 0 
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Figure 6.7 


Figure 6.8 
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Nous avons également (voir figure 6.9) 


hQ 
Mb a 


7 7 (6.54) 


donc r 
Qp = UT (6.55) 


Figure 6.9 


6.5 Approche vectorielle et les équations d’Euler 
[E 


En pratique on cherche toujours l’approche la plus simple. Par exemple si le solide 
étudié n’est pas attaché en un point alors on utilisera comme point d’ancrage du système 
intrinsèque celui qui décompose l’énergie cinétique en un terme de translation plus un 
terme de rotation. Évidemment si le solide a un point fixé dans un référentiel inertiel, un 
tel point devient naturel pour définir une origine. 

Il n’en demeure pas moins que pour établir des équations de mouvement nous devons 
nous référer à un référentiel inertiel, quitte à traduire ensuite en termes de quantités me- 
surées dans un référentiel non inertiel si on choisit de le faire. Notre solide a 6 degrés de 
liberté que nous avons noté en page 99 par X, Y, Z, £1, £2, £3. Donc en principe 6 équa- 
tions de mouvement. Nous les établirons ici vectoriellement, à Newton en considérant P 
et 1. Rappelons que P est le moment linéaire total 


P= 5p =MV = Pom. (6.56) 


part. 
Nous avons déjà vu au chapitre 1 que 


È + =F: la force extérieure. (6.57) 

On solutionne ce problème exactement comme dans la mécanique d’une particule. Il n’y 

a rien de neuf ici. Considérons maintenant Í, i.e. la variation dans le temps du moment 

cinétique/angulaire. Dans un référentiel inertiel ou du laboratoire que nous notons par 
l'indice L 

dl 


= —| =N ; 
z es 


L 
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où N est le torque extérieur. Nous référant toujours à la figure 6.1 de la page 99 nous 
pouvons écrire 

ll = RXMV+Ilsr (6.59) 
où l|çp est mesuré par rapport au point O du solide qui sert à y ancrer le référentiel 
intrinsèque mais le F de l’indice SF signifie que les composantes de sont prises par rap- 
port à un référentiel fixe (F), i.e. qui ne tourne pas avec le solide. À ce point-ci cette 
nuance n’est pas significative puisque a une existence physique indépendante du réfé- 
rentiel par rapport auquel nous en mesurons les composantes. Par contre, dans ce qui 
suit nous allons le dériver par rapport au temps et là ça deviendra significatif, parce que 
dans un premier temps nous voulons éviter les dépendances dans le temps provenant de 
la rotation des axes. Nous avons donc: 


il = RxMV +RxMV + i| =N 
L ~a SF 
=0 pcq R=V 
z RxMV+ i| =N. (6.60) 
SF 


Ici nous allons nous limites aux cas où V = 0, ce qui veut dire que V = constante 
ou V =0, ce dernier cas étant utile lorsque le solide a un point fixe dans un référentiel 
inertiel. Nous gardons donc 
Ten, (6.61) 
SF 
Nous avons cependant pris l’habitude de décomposer ls selon un système d’axes qui 
tourne avec le solide avec la vitesse instantanée Q. Il s’agit toujours du même vecteur, 
qui a toujours la même réalité physique mais tout simplement dans le calcul de la dérivée 
par rapport au temps de ses composantes, nous devons maintenant tenir compte du fait 
que ces axes tournent. Ils sont donc accélérés et de ce fait le référentiel n’est par inertiel. 
Il a été vu en Mécanique Classique I que nous avons alors 


i| Š i| +92 xls =N. (6.62) 
SF S 


Gardant en mémoire que nous mesurerons toujours les composantes de ls selon les axes 
du référentiel intrinsèque qui tourne avec le solide, nous laissons tomber l’indice S et 
nous écrivons simplement 


d 
ren (6.63) 
Si les axes du système inertiel coïncident avec les axes principaux nous avons 
li = liQ; (pas de somme sur i). (6.64) 
Ceci permet d’écrire (6.63) 
Qi 
LUE + Eijk NkIk = N; (pas de somme sur i) (6.65) 


où 6;,x est le ”’tenseur antisymétrique”. En termes plus explicites nous avons ici 3 équa- 
tions 


LM - MDI) = M 
DO — Q3 Qı (I3 = L) = Mo (6.66) 
B93 — QOQ- h) = Nz. 


Ce sont les équations ď’ Euler 
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Évidemment s’il n’y a aucun torque extérieur, i.e. N = 0 alors (6.66) se réduit à 


D -R2 -I) = 0 
D- UBU -hA) = 0 (6.67) 
BÊ- M —-b) = 0. 


À l’aide de ces équations voyons si nous pouvons refaire le problème de la toupie 
symétrique de la page 110. Rappelons que nous avons Jı = I2 # 13 # 0. Avec Ji = D, 
(6.67) nous donne 


IÊ; = 0 — Q; = constante. (6.68) 
De plus les deux premières équations de (6.67) deviennent 
Ô = wh (6.69) 
Ô = wM (6.70) 
où w = Q3 G, et 
QG = wô (6.71) 
= wQ (6.72) 
et donc 
Q(t) = Acos(wt+ 8) (6.73) 
Q(t) = Asin (wt+ 8) (6.74) 


Donc 4/02 + Q2 = A une constante qui est la longueur de la projection de la vitesse 
angulaire dans le plan x10x2 (voir figure 6.10). D’autre part Q3 (projection de Q sur 
l’axe de la toupie) est une constante donc c’est l’ensemble de Q qui tourne autour de l’ axe 
de la toupie à vitesse angulaire w. A première vue ce résultat ne semble pas correspondre 
au résultat de la page 110 mais on s’intéressait alors à la vitesse de précession de l’axe 
Ox3 par rapport à un axe fixe donné par la direction de 1. Ici nous avons tout exprimé 
(tous les vecteurs) selon leur composantes mesurées sur des axes tournants. Il faudrait 
savoir faire le bien entre les deux. C’est ce que la méthode des angles d’Euler va nous 
apprendre à faire. 


Figure 6.10 


En première approximation on peut appliquer ce même genre de raisonnement au 
mouvement de la terre si on considère que (le torque dû à) la force soleil-terre reste faible, 
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donc V ~0 est une approximation raisonnable. Aucun torque extérieur ne s’applique 
dans ce cas et on peut utiliser (6.67). C’est que la terre n’est pas tout-à-fait sphérique (ni 
rigide) de telle sorte que 
LR < Ís (6.75) 

du fait de l’aplatissement de la terre si on prend l’axe Ox3 comme son axe de rotation. 
En fait E. 

A ~ —0.003 (6.76) 

2 

et par conséquent, appliquant les résultats du problème précédent nous avons une fré- 
quence 


& —0.00393 (6.77) 


or eus 7 1 jour et par conséquent on devrait avoir 
3 


2 
= ~ 333 jours. (6.78) 


De fait, un tel mouvement est observé mais son amplitude est très faible, l’ amplitude du 
déplacement du pôle étant de l’ordre de 5m. D’ autre part, la période est réellement ~427 
jours mais cette différence est imputable au fait que la terre n’est ni rigide ni uniforme. 


6.6 Angles d'Euler et approche Lagrangienne 


Il nous manque encore un outil, celui qui nous permettrait de faire systématiquement 
le passage entre référentiels tournants et immobiles. Ce même outil faciliterait aussi l’é- 
criture d’un Lagrangien. Ici nous ne nous intéressons qu’aux rotations. Nous faisons donc 
coïncider les origines de OXY Z et de Ox1x2æ3. Ceci semble indiquer que V = 0 mais 
en fait s’applique tant que O peut être l’origine d’un référentiel Ox1x2xz3 de OXY Z 
c’est une rotation. Cette rotation s’effectue instantanément par rapport à un axe. Malheu- 
reusement cet axe peut varier en direction avec le temps. Les angles d’Euler permettent 
de représenter cette rotation (en fait toute rotation) comme une séquence de trois rota- 
tions successives mais par rapport à des axes dont il nous est possible de garder la trace. 


Sur la figure 6.11 nous avons indiqué le repère fixe OXY Z et le repère tournant 
Ox1x2t3. On y remarque de plus l’axe ou la droite ON qui est la droite de contact 
entre les plans XOY et x10x2. On l'appelle la ligne nodale. Langle y est l’angle entre 
laxe OX et cette ligne nodale suite à une rotation dans le plan XOY, i.e. par rapport à 
l’axe OZ. Langle v est l’angle entre cette même ligne nodale et l’axe Ox:, mesuré dans 
le plan x,O%2, i.e. par rapport à l’axe Ox3. Quant à l’angle 0, c’est simplement l’ angle 
entre laxe OZ et l’axe Ox3, il correspond à une rotation par rapport à l’axe ON. Ce sont 
ce axes par rapport auxquels sont effectuées les rotations que nous avons représentés par 
p, xet b. Décomposant ces trois vecteurs vitesse selon les axes mobiles de Ox1£2£3 


nous avons , f | | | 
01 = 0 cosp; 02 —060sinv; 03 —<0 (6.79) 

pı = psinf siny; p= psinð cosy; 1 — pcosb (6.80) 

yı=0; y2=0; y=% (6.81) 


Les composantes Q1, Q2, et Q3 de Q sont simplement les sommes des composantes 
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Figure 6.11 
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respectives. Par exemple | | 
Qı = 01 +1 + Y1, (6.82) 
c’est-à-dire 


Qı = ÛcosŸ + ġsinð siny 
Q = sin + ġsinð cosy (6.83) 
Qs = ġcosð +y. 


Nous avons maintenant complété l élaboration des outils qui sont nécessaires pour at- 
taquer plusieurs problèmes impliquant le mouvement du solide. En choisissant les axes 
Ox1x2T%3 comme les axes propres du solides, nous pouvons spécialiser ces expressions 
pour écrire Tiot en fonction des angles d’Euler. Dans ce cas nous avons 


Tiot = 0% + 203 + Ta (6.84) 
Par exemple, pour la toupie symétrique où J} = I2 = 13 = I nous avons 
Ti = L (o + GE? + Ÿ + 20 cos 0) | (6.85) 
Pour la toupie symétrique où Jı = 12 Æ Iz nous avons 
Tha 4 (è +g? sin? 0) "E z ( + @cos o) | (6.86) 


6.7 Exemple 


© 


Comme exemple d’application retournons au cas de la toupie symétrique libre I4 = I> # 
T3. Choisissons l’axe OZ du référentiel pour qu’il coïncide avec la direction de 1 donc 


1—7/Z: une constante. (6.87) 
Trivialement 
l3 = l cos 0 = 03. (6.88) 
Le rôle du Lagrangien sers joué par Tot de la toupie symétrie ci-dessus 
l /:2 LB /: 2 
L= Tu = + (è + ọsin? 0) +2 ( J &cos0) | (6.89) 
On constate que Ÿ et & sont cycliques et par conséquent Py et Ps sont des constantes 
OL . 
Py = o% =Í; (% + cos o)- IQ = constante (6.90) 
— ee” 
Q3 
puisque pour des variables cyclique 2L = 0, 4 (2 = Py = 0. Par conséquent 


lcosð = 1303 est une constante, donc cos 0 est une constante, donc 0 est constant. 
Calculant 


OL | 
Pe = gp nesin + (ù + 50050) cosÿ 
= Jigsin? 0 + I3Q3 cos 0 = constante (6.91) 


où J4, 0 et 1:Q3 sont des constantes donc & = constante. Langle y indique une rotation 
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du référentiel intrinsèque, donc du solide, par rapport à un axe OZ fixe choisi dans la 
direction de 1. C’est donc une précession du solide ou si on préfère de l’axe de symétrie 
Ox3 du solide par rapport à cet axe fixe. De toute évidence il s’agit de la même précession 
que celle étudié en page 110. Elle est toutefois différente de celle étudié en page 114 où 
on étudiait une précession de laxe de rotation instantanée par rapport à l’axe Oxg i.e. 
par rapport à l’axe de symétrie du solide. Dans ce dernier cas nous avons vu que l’axe de 
rotation de la terre, qui ne correspond pas à l’axe de symétrie de la terre (axe des pôles), 
tourne autour de cet axe de symétrie avec une période de rotation de 1 an. Ceci n’interdit 
pas à l’axe de symétrie de la terre de précesser par rapport à un axe approximativement 
fixe qui serait celui du moment angulaire. Une précession libre de ce type se ferait avec 
une vitesse angulaire %. Il convient donc d’étudier celle-ci d’un peu plus près. En fait 
nous connaissons la réponse 


l 
p= T’ let Iı = constantes. (6.92) 
1 


Le vecteur l est constant et on peut le décomposer selon le système d’axes que l’on veut. 
Choisissant le référentiel intrinsèque nous savons que 


1 = ARI + BLR + IQ (6.93) 
où Q; = Q;X;. Or ces trois axes sont orthogonaux et ici I2 = J4, donc 
P=1.1= (Q? +923) + BA = constante (6.94) 


que l’on peut récrire en fonction des angles d’ Euler 
; 2 
P=1.1= 7? (i + %’° sin? 0) + (v + cos 0) = constante (6.95) 


où Å = 0 (8 = constante) et nous savons déjà par la définition de Py que la deuxième 
terme est simplement égal à 1303 = 13, donc 
P= p’ sin? 0 + 2. (6.96) 
Sachant que l3 = l cos et isolant &? nous obtenons 
P(l—cos 0) P 
2 
RE e 6.97 
j I? sin? 0 I (ep 


ou encore 
l b è kb% 


2 = — = —— = = —. 6.98 
P IL cos I; cos ( ) 
Par exemple, lorsque Jı œ~ 13 et que la vitesse de rotation est essentiellement w ~ Q3 


(le cas de la terre) alors 


I 
bax TO xw. (6.99) 


1 
Pour la terre w & ent Ainsi dans le cas libre il y a une précession de l’axe de symétrie par 
rapport à un axe dont la direction est donnée par le vecteur moment angulaire constant 
l, et cette précession a une vitesse angulaire donnée par %. Dans le cas de la terre cette 
précession fait que l’axe de symétrie de la terre tourne autour de 1 ~ une fois par jour. 
Dans le cas de la terre, il y a une autre précession causée par un torque cette fois, résultant 
de l’action de la lune et du soleil sur une terre non sphérique. Cette précession, dite des 
équinoxes est de plus grande amplitude mais a une période d’environ 26,000 ans. Elle est 


donc faible et il était raisonnable en première approximation de la négliger et de parler 
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ce cas libre. 


La situation n’est pas encore lumineuse. En effet en page 114 nous avons étudié la 
précession libre à l’aide des équations d’Euler et obtenu une période de précession de 
l’ordre d’une année pour la terre. C’est suffisamment différent de la période d’environ 
une journée obtenue ci-dessus, apparemment dans les mêmes conditions, pour se poser 
la question de la cohérence entre ces deux résultats. Physiquement, la situation est effec- 
tivement différente. La période d’environ 1 journée ci-dessus est celle de la précession 
de l’axe de symétrie de la terre par rapport à l’axe défini par la direction de 1. En page 
114 nous avons étudié la précession du vecteur Q| , tel que 


QD =M+® où Q =Q, +Q; (6.100) 


par rapport à l’axe Oxz, lui-même mobile. Cette précession apparaît dans la dépendance 
dans le temps de Qı et Q2 


Qı = Asin(Qht+ô) (6.101) 
M = Acos(NQprt+ ô). (6.102) 
Étudions donc ici ce que nous obtenons pour 94 (t). De (6.83) 
Qı = cosy + ÿsin 0 sin Y. (6.103) 
Ici 0 = constante = 09, donc = 0 et ġ = + et donc 
l 
Qı = ġsin 0 sin Y = 7, $20 sin Ÿ. (6.104) 
1 
Recherche de y(t) : 
Rappelant la définition de py déjà obtenue (6.90) nous avons 
Py = Is (ù + pcos0) 
= BR = l = lcosð. (6.105) 
On peut y isoler Ÿ en remplaçant p = T 
. 1 l 
p = — [| lcosð — Is cos0— (6.106) 
I3 IEI 
lcos0 (Å — I l (L — I 
BE A E EE) (6.107) 
Iz l L h 
avec l3 = 1303 nous avons 
. L-I L-I 
as, — y ge es) — Dito (6.108) 
1 1 


Nous avons donc pour 94 (t) 


l h-I: 
Q(t) = — sin 0o sin o, ds); + bol. (6.109) 
lh l 
Nous identifions donc, comparant (6.109) à (6.101) et (6.102) 
lsin 0 h-I 
|A] = 0l; =o [ol = a, , (6.110) 


Nous avons donc le même résultat pour Q2. Ainsi, le petit vecteur Q| dans le plan (non 
fixe) 4102 tourne autour de l’axe de symétrie de la terre avec la fréquence Qpr dont la 
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période est de l’ordre d’un an. 


6.8 Mouvement d'une toupie symétrique pesante à un 
point fixe 


C’est l’exemple typique de tout manuel de mécanique . Ayant complété l’étude du 
mouvement du solide libre, nous nous attaquons au mouvement du solide soumis à un 
torque. Pour ce faire nous choisissons l’exemple le plus simple d’une toupie symétrique 
dont la pointe est fixe et placée dans un champ gravitationnel uniforme. Comme la pointe 
est fixe, nous allons l’utiliser comme l’origine à la fois pour le système inertiel OXY Z et 
pour le système intrinsèque Ox1x2x3 qui lui, tourne avec la toupie (voir la figure 6.12). 
À priori ceci semble poser un problème puis que nous avions réussi en (6.40) à séparer 
T à condition que l’origine du référentiel intrinsèque coïncide avec le C.M. du solide, 
ce qui n’est pas le cas ici. Par contre ici, les deux origines coïncident et donc R = 0, 
R= V =0 et il ne reste que Tror. De plus, le corps étant symétrique, l’énergie cinétique 
de rotation peut s’écrire 


I I 
Ta = z (93 +0) + F% (6.111) 


avec I} = I> mais ici Iı n’est pas égal au moment d’inertie calculé par rapport à l’ axe 
principal 1: 1, qui lui, passe par le C.M. Cependant, par le théorème des axes 


Figure 6.12 


parallèles on calcule trivialement que le J} qui apparaît ici est simplement 
L = Li + Mh? (6.112) 
où M est la messe de la toupie et h la distance séparant la pointe du C.M.. Auparavant 


nous n’avions pas l’habitude de spécifier ” pr ” pour alléger l’écriture nous retenons le 
symbole même si ce n’est pas une moment d'inertie par rapport à un axe principal. 
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En terme des angles d’Euler (6.111) est identique à (6.86) puisque Z2 = Ji 
l /: Le fe 2 
rs (9° + g? sin? 0) +3 (ù + pcos) (6.113) 
tenant compte que J; est défini en (6.112). Comme le champ gravitationnel est constant 
on peut représenter son effet comme une force M g appliquée au C.M.. Puisque la pointe 
est fixe, cette force génère un torque ainsi l n’est plus une constante du mouvement. 


Pour écrire le Lagrangien nous n’avons besoin que de l’énergie potentielle résultant de 
la présence de ce champ de force, c’est simplement 


V = Mgh cos 0 (6.114) 


et alors 
; ; 2 
L= A (0° + g? sin? 0) +3 (ù +6 0050) — Mgh cos 6. (6.115) 


Ici, comme d’ailleurs dans le cas libre, y et Ÿ sont cycliques et par conséquent py et 
Py sont des constantes du mouvement. D’ autre part de façon générale = N (le torque) 
où N =r x F. Or ici la force est parallèle à OZ et donc l, = constante et de plus cette 
force est attachée à un point se trouvant sur l’axe Ox3 et donc l3 = constante aussi. Nous 
avons déjà remarqué (dans le cas libre) en (6.90) que 


OL 
Py = =I C + cos 6) 
Y EN % 3 (Y +o 
= IQ = l3 = constante (6.116) 
identifiant p à l3. D’autre part, il est également évident que 
L 3 
Po = - = I ýsin? 0 + Is (% + & cos6) cos 0 
p 
= l, = constante. (6.117) 
Combinant ce deux expressions nous obtenons 
l, — lz cos 0 
D = ——— 6.118 
I sin? 0 ( ) 


et combinant (6.116) et (6.118) nous isolons ÿ 
. la cos 0 
= — — —— (lL — l3cos0). 6.119 
Y Fe L ae 0 ( 3 ) ( ) 
Si on connaît O(t), on peut en principe intégrer (6.118) et (6.119) pour obtenir y(t) et 
p(t), ce qui donnerait la solution complète du problème. Pour obtenir 0 (t) on peut utiliser 


l’équation de Lagrange 
d (ƏL ðL 
— |=] -= =0 6.120 
(Se) -% (6.120) 
dont on élimine les G et E: à Te de (6.118) et (6.119) pour avoir une équation différen- 
tielle uniquement en 8, Ÿ, 6. Isolant Ê dans une telle équation nous donne 


(l; — la cos 0) (L; cos 0 — l3) o 
=; — Mghcos = —— Verl0 6.121 
I sin? 0 I 00 er(0) ( ) 
ce qui nous permettrait d'identifier (par définition) un potentiel efficace pour le mou- 
vement en 0. Est-il besoin de dire qu’intégrer une telle équation différentielle est tech- 


1,0 = 
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niquement assez difficile. Il existe cependant une autre approche qui nous donne VW 
plus facilement. Rappelons qu’en présence d’un torque, l n’est plus que une constante 
de mouvement mais l’énergie E continue à être une constante du mouvement là où 


E=T +V. (6.122) 
Utilisant (6.113, 6.114, et 6.116) nous avons 


21; 


= const. 


l 2 l .2 . 2 12 
E = z? + z sin 0+ = +- Mghcosð. (6.123) 
SZ 


Définissant E’ = E — oa et à l’aide de (6.118) nous avons 


L P % (l; — l3 cos 0)? 
2 21, sin? 0 

= To+V (6.124) 
définissant ainsi Vo, un potentiel efficace pour l’étude du mouvement en 0. Ici aussi 
l'intégration mène à des intégrales elliptiques et on perd les propriétés du mouvement 
dans les méandres techniques. Heureusement il est possible de déterminer qualitative- 
ment les propriétés intéressantes de ce mouvement. Avant de procéder cependant remar- 
quons que si Ø = 0, Vọ semble exploser à cause du facteur sin? 9 au dénominateur. En 
fait, à 0 = 0 les axes Ox3 et OZ coïncident, et l, — l3 cos 0 = 0. C’est donc une dé- 
termination. On peut vérifier, par la règle de l Hôpital que le terme litigieux de V5 — 0 
lorsque 0 — 0. 


E = + Mgh cos 0 


Faisons la transformation de variable 


u=cosð = ù= -—bsinð (6.125) 
remplaçons dans (6.124) et isolons ù? pour obtenir 
2E'  2Mghu (l; — lu)? 
2 _ A hp 2 NE 
U = ( T, l ) (1—u*) T (6.126) 
qui est de la forme 
ù? = (a — bu) (1 — u?) — (b — au)? = f(u). (6.127) 


La fonction f(u) est un polynômes cubique en u dont le coefficient de u3, 3 > 0. Donc 
f(—20) — —o et f(+oo) — +00 avec deux extrema entre ces deux limites. Puisque 
u = cos 6, seul le problème de u compris u = —1 et u = +1 nous intéresse. D’autre 
part en (6.127) 

fu) = & > 0 (6.128) 
et nous sommes donc limités au domaine entre u et u2. Pour qu’une situation physique 
existe, il faut que ces deux conditions soient remplies. Si tel est le cas et puisque 0 marque 
l’angle entre la verticale et l’axe de symétrie de la toupie, cet axe de symétrie aura, par 
rapport à la verticale, un angle qui oscillera entre les angles 0 et 02 où 


cos 01 = u1, COS 02 = u2. (6.129) 
C’est ce qu’on appelle une nutation. Rappelons qu’en (6.118) nous avons obtenu pour 
. EL —l3cos0 


= 227 6.130 
£ I sin? 0 ( ) 
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Selon les conditions initiales qui déterminent l, et l3 et selon le domaine de variation 
permis pour 0, donc pour cos 0 on peut identifier trois scénarios différents pour 


1.1, — l3cos0 > 0 pour tout le domaine de variation de 0. Alors @ ne change pas de 
signe et la précession est continue bien que de vitesse variable. (La même chose est 
valide si l; — l3 cos 0 < 0 dans tout le domaine, on change simplement le signe de 
p). 

2. l, — l3cos0 > 0 change de signe entre 0, et 02. Parce que la fonction cos 0 est 
monotoniquement croissante entre —1 et +1, alors ġ(01) aura le signe inverse de 


(01). La précession continuera de se produire mais avec des mouvements de va-et- 
vient. 


3.1, —l3 cos 0 > 0 ne change pas de signe dans le domaine 04 < 0 < 02 mais s’annule 
soit à 0. soit à 02. Dans ce cas la précession est toujours dans la même direction mais 
marque un temps d’arrêt lorsque la nutation atteint une de ses valeurs limites (soit 01 
soit 02) là où  s’annule. 


Il est habituel de représenter ces trois situations à l’aide de figures simples. On dessine 
une sphère qui est celle que l’extrémité libre de la toupie peut générer (puisqu'elle a une 
pointe fixe) et sur la surface de cette sphère on trace la trajectoire que la pointe libre y 
dessinerait. On a alors les trois sphères de la figure 6.13 respectivement. 


Figure 6.13 


Tous ceux qui se sont amusés avec une toupie ont pu constater la chose suivante. Si 
on démarre la toupie avec une vitesse élevée et une faible inclinaison par rapport à la 
verticale, alors elle dort, son axe demeurant pratiquement vertical. La friction aidant sa 
vitesse diminue jusqu’à un pointe où la toupie devient presque brutalement instable. Pour 
étudier ce phénomène rappelons la définition de V2 en (6.124) 

2 
V; = CONTE (6.131) 
21, sin“ 0 
On constate d’abord que Va (0 = 0) = Mgh et que le deuxième terme de V9 est répulsif 
parce que cos 0 est maximum à 0 = 0. Si la position de la toupie est stable à 0 ~ 0, 
c’est que Vọ doit avoir un minimum à 0 ~ 0. Pour étudier ce phénomène faisons une 
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expansion de Vọ valable aux petits angles (Taylor) 


OV» 0? V 
Vo = Va (0) + 0 — — 6.132 
0 o( ) 00 j 2 302 PE ( ) 
Se rappelant qu’à 0 = 0, OZ et Ox3 coïncident, donc l3 = lz 
Vo(0) = Mgh = constante sans intérêt (6.133) 
l-l l3 — l; cos 0 ; 
D, a EL e 
00 loo I sin? 0 9=0 
Alors utilisant la règle de 1’ Hôpital, 
Hu 2 =0—0=0: extremum à 0 = 0. (6.135) 
9-0 00 
Par ailleurs 
3? Vy 12 (1 — cos 0) > 
TA = La + cos” 0 — 3 cos 0) F — Mghcosô|,_, (6.136) 


Le premier terme donne, utilisant une fois la règle de l’ Hôpital, 
2 (1- 
m6 (1 — cos 8) 
0-0 11 sin*0 
une indétermination. Une double application de la règle de 1’ Hôpital donne cependant 


(2 + cos? 0 — 3 cos 0) = i (6.137) 


Vo 21 
= >- — Mgh. 6.138 
8 |, La g ( ) 
Au total donc 
E TEA $ — Mgh) + (6.139) 
9 — g 2 AT; g É 
2 
L'extremum à 0 = 0 sera un minimum si 4 — Mgh > 0 donc sil? = 120% > 4l, Mgh 
ou encore An 
g> LT (6.140) 


Tant que la vitesse de rotation de la toupie i dort, essentiellement (3, satisfait cette 
condition, la position verticale de la toupie est stable. Lorsque la friction fait tomber la 
vitesse sous cette limite l’extremum de Va à 0 ~ 0 devient instable et le mouvement 
devient rapidement désordonné. Le tout est en fait le résultat d’une compétition entre le 
terme de Vo, M gh cos 0, qui tend à faire tomber la toupie, et un terme qui provient de 
la rotation de la toupie et tend à la garder verticale. Plus la vitesse de rotation augmente, 
plus la stabilité est grande et moins l’effet de M gh cos 0 (donc du torque extérieur) est 
important. En fait, on peut dire qu’à grande vitesse la situation ressemble au cas libre. 


6.9 La toupie asymétrique libre: problème de stabilité 


Nous avons ici 11 Æ I> Æ T3. Posons ici 1 < 13 < 13. Dans la cas libre nous 
avons essentiellement conservation de 1 donc de 1? et de E. Développant selon les axes 
principaux nous aurons 


P =? = RO? + 1502 + I = constante (6.141) 
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1 1 
E = ino + ze% + zB% = constante. (6.142) 
On peut aussi écrire 
l =Ê ++ (6.143) 
HE 
= + + — .144 
211 2b 21 i 


Traçant les trois axes orthogonaux de coordonnées, l1, l2 et l3, on constate que la pre- 
mière équation définit la surface d’une sphère de rayon l alors que la deuxième définit 
la surface d’une ellipsoïde de demi-axes de longueurs V2E11, V2EI et V2EI3. Les 
deux équations doivent être satisfaites simultanément. Ainsi l’ extrémité du vecteur l ne 
pourra suivre que les courbes d’intersection de ces deux surfaces. Il est également clair 
que 

JE < < 2E (6.145) 
Si? = 2E1; (ou = 2EJ3), l est minimum (maximum) et 1 est selon l’axe 1 (axe 3). 
Toutes les valeurs intermédiaires sont permises. On voit sur la figure 6.14 un série de 
trajectoires tracées par la pointe de 1 pour différentes valeurs de l allant croissant de la 
trajectoire 1 où l est près de la valeur minimum jusqu’à la trajectoire 6 où l est près de 
sa valeur maximale. Dans le cas de la courbel, l est près de la valeur minimale donc 


Figure 6.14 


et l est surtout selon l’axe 1 ce qui indique une rotation autour de l’axe 1. On voit que la 
trajectoire est fermée et que 1 dérive peu de la direction 1. La rotation autour,ou presque, 
de l’axe 1 est stable. La même chose s’applique lorsque l est près de la valeur maximale 
alors que l est près de la direction de l’axe 3, indiquant une rotation autour de l’axe 3. Ici 
encore l trace une trajectoire fermée autour de l’axe 3. Mais tel n’est pas le cas lorsque la 
rotation se fait autour de l’axe 2 parce que les trajectoires tracées par l et qui passent près 
de ou par la direction de l’axe 2 (rotation autour de cet axe) ne sont pas fermées autour 
de l’axe 2 mais se promènent tout autour de l’ellipsoïde, passant même par les parties 
négatives de l2. Nous en concluons qu’une trajectoire initiée autour de l’axe 2, celui dont 
le moment d’inertie a la valeur intermédiaire, entre J} et I3, sera instable. C’est ce qu’on 
constate expérimentalement lorsqu'on fait tourner une raquette ou un livre (gardé fermé 
par un élastique) par exemple. 
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Cette explication est clairement plus qualitative que quantitative mais elle nous donne 
néanmoins une image raisonnable du phénomène. 


———< }— 


qq 
Annexe A: Notations, conventions... 


A.1 Notations et conventions 


Dans cet ouvrage, une certain nombre de conventions ont été adoptées pour faciliter 
la lecture. Les vecteurs sont notés par des caractères gras 


Xx,r,V,F,.. (A.1) 
L'alphabet grec est utilisé fréquemment: 


Majuscule Minuscule Prononciation 
alpha 
bêta 
gamma 
delta 
epsilon 
zeta 
eta 
theta 
iota 
kappa 
lambda 
mu 


ax © © 


D 
m 


omicron 
pi 
rho 
sigma 
tau 
upsilon 
phi 
psi 
chi 
omega 


D 
e a aD IJ OMNE OZ e e 


F 
6 


DHÉRRAMTOOMSÉR>-RS ONNEDL ERR 
x & 


E 
9 
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A.2 Systèmes de coordonnées 


Coordonnées cartésiennes 


az Le DOS 
ee | a | 
__ | P&yZ) ns êy g 
€ | | 
io y 
| y 
<L A o 
s e 
x plan a 
Xx=X] | L 
=y1 
Figure 6.1 


Les vecteurs unitaires d’un système de coordonnées cartésiennes â,, ây, â, ont les 
propriétés suivantes 


âr xây = À 
âx à = ê (A.2) 
â, X âr = ây. 


Un vecteur A dans ce système de coordonnées s’exprime souvent sous la forme de ses 
composantes A = (Az, Ay, Az) ce qui représente la somme vectorielle 


A = â, À, + ây Ay + à: 4. (A.3) 


Les éléments de longueur, dl = (dx, dy, dz), de surface, (ds, dsy, ds.) , et de volume, 
dv, sont respectivement 


dl = â,dx + &ydy + &.dz (A4) 
dsx = dydz 
dsy = dxdz (A5) 


ds, = dedy 


A2 Systèmes de coordonnées 129 


dv = drdydz. (A.6) 


Remarque 10 Dans ces notes, nous allégeons la notation en prenant 
âz, ây, â; =i,j,k (A.7) 


mais dans la littérature, les vecteurs unitaires Âz, ây, à; s'écrivent aussi souvent sous 
les formes variées 


l 


32 


XI 
» © 


, 


a 


X,y,Z 


êz, êy, êz. 


Coordonnées cylindriques 


cylindre 
r=r 
1 


Figure 6.2 


Les vecteurs unitaires d’un système de coordonnées cylindriques â,, â¢, â, ont les 
propriétés suivantes 


â xâ = À 
âp XÂ = ê (A.8) 
â xâ, = êy. 


Un vecteur A dans ce système de coordonnées s’exprime souvent sous la forme de ses 
composantes A = (4,, Ag, Az) ce qui représente la somme vectorielle 


A = â, A, + âp Ag + â, Az. (A.9) 
Les éléments de longueur, dl = (dr, d@, dz), de surface, (ds,, dsġ,dsz), et de volume, 
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dv, sont respectivement 


dl = â;dr + â,rdp + â.dz (A.10) 
dsr = rdpdz 

dss = drdz (A.11) 
ds, = rdrd® 

dv = rdrdpdz. (A.12) 


Les relations de transformations de coordonnées cylindriques à coordonnées cartésiennes 
sont les suivantes: 


= Trcoso® 
= rsin® (A.13) 
2 = :ÿ 
et inversement 
r = VER 


® = arctan Ÿ (A.14) 


Coordonnées sphériques 


sphère 
R=R; 


Figure 6.3 


Les vecteurs unitaires d’un système de coordonnées sphériques âp, âp, â4 ont les 
propriétés suivantes 


A 


âr X âọ = az 
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âg X gs = ÂR (A.15) 
âp X ÂrR = À. 
Un vecteur A dans ce système de coordonnées s’exprime souvent sous la forme de ses 
composantes A = (Ar, 49, À;) ce qui représente la somme vectorielle 
A = ârRAR + âp Ao + âg Ag. (A.16) 
Les éléments de longueur, dl = (dR, dọ, dz), de surface, (dsr, dso, ds4) , et de volume, 


dv, sont respectivement 
dl = â,dr + âg Rd + âg R sin 0d@ (A.17) 


dsr = R?sin0dOdé 
dsọ = RsinO0dRdp (A.18) 
dse = RdRdO 
dv = R? sin 9d Rdbdé. 
Les relations de transformations de coordonnées sphériques à coordonnées cartésiennes 
sont les suivantes: 


x = Rsinôcos® 
Rsin 0 sin @ (A.19) 
z = Rcosô 


et inversement 
r = / £2 + y? + 22 
2 2 
ve ty (A.20) 


zZ 


0 = arctan 


© 
I 


arctan CA 
£ 
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A.3 Aide-mémoire 


Mécanique lagrangienne 


L équation d’Euler-Lagrange pour un Lagrangien L(qi, qi, t) : 


d (ƏL) 0, 
dt \ di ðq 


(À compléter par létudiant.) 


Corps solide 


Moments d'inertie I par rapport à l’axe de symétrie: 
Tige mince p/r extrémité iM R? 
Tige mince p/r centre SM R? 

Sphère pleine: 2M R? 

Sphère creuse ou coquille mince: 3m R? 

Disque ou cylindre plein: 5M R? 

Cylindre creux ou anneau mince: MR? 

Anneau épais: + M(R?, + R2,) 


Dynamique 
T=L-la 
L = Iw 
Tot = Liw? 
2 


Condition de roulement sans glissement 
v=wR 
Théorème des axes parallèles: 
I=Icm +M- œ 
Théorème des plaques minces: 
I, = Iz + Iy 
Constantes usuelles: 


Accélération gravitationnelle: g = 9.8 m - s7? 
Rayon terrestre: R = 6378 km 


Vitesse angulaire terrestre: w = 7.27 x 1075 rad: s71 


A.0 133 


A.4 Références 


Les notes couvrent une partie de ce qui est traité dans les volumes suivants et ceux-ci 
peuvent être utilisés à titre complémentaire. 
1. Classical Mechanics, H. Goldstein, 2° édition, Addison-Wesley (1980). 


2. Mécanique., L. Laudau et E. Lifchitz, 4e édition, Éditions MIR. 
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